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Dear  Mr.  Ulsh: 


This  is  the  final  report  for  the  U.S.  Army  Research  Office  Grant  DAAL  03-86-k-0090. 
According  to  the  reporting  instructions  (ARO  Form  18,  page  S)  technical  material  which  was 
already  reported  will  not  be  included  here. 

The  thrust  of  the  research  under  this  grant  has  been  to  develop  the  methodology  for 
the  modeling  and  design  antenna  arrays  and  microstrip  discontinuities  for  microwave  circuit 
applications.  I  believe  that  certain  outstanding  results  have  been  obtained  during  this  period 
and  these  are  highlighted  in  what  follows: 

H-Y.  Yang’s  Ph.D.  thesis  involves  several  key  contributions  in  the  subject  areas  of 
modeling  microstrip  discontinuities,  microstrip  transitions,  and  the  synthesis  of  microstrip  di¬ 
pole  arrays.  This  woik  has  generated  several  journal  publications  within  which  one  finds  a 
fundamental  contribution  to  the  understanding  of  the  above  mentioned  subjects.  The  analysis 
and  the  generated  computer  programs  serve  as  practical  tools  for  the  design  of  microstrip  cir¬ 
cuits  and  microstrip  dipole  arrays.  The  methodology  was  substantiated  in  each  case  with  an 
experiment.  The  journal  publication  "Design  of  Transversely  Fed  EMC  Microstrip  Dipole  Ar¬ 
rays  Including  Mutual  Coupling"  IEEE  Transactions  on  Antermas  and  Propagation,  Vol.  38, 
No.2,  February  1990  is  only  one  of  many  examples  from  Yang’s  excellent  Ph.D.  thesis  work. 
This  is  a  seminal  paper  since  it  is  the  fint  publication  where  a  full  synthesis  procedure  is 
developed  foe  the  design  of  transversely  fed  electromagnetically  coupled  dipoles.  The  solu¬ 
tion  is  complete  as  it  accounts  for  all  substrate  effiects  and  mutu^  coupling.  It  demonstrates 
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that  mutual  coupling  contributes  25%  of  the  total  value  of  the  array  input  active  admittance 
and  it  therefore  can  not  be  neglected  The  design  procedure  evolved  into  the  fabrication  of  a 
microstrip  dipole  array.  The  experimental  results  show  very  good  agreement  with  theory. 


An  evolution  of  Yang’s  work  has  led  to  the  understanding  of  microstrip  discontinuity 
effects.  In  particular  the  theory  has  been  extended  to  model  microstrip  bends,  T-junctions, 
four  ports,  etc.  In  each  case  the  algorithms  account  for  a  precise  description  of  energy  loss  at 
discontinuities  due  to  radiation  loss  and  surface  wave  loss.  These  models  have  also  been  sub¬ 
stantiated  with  experiment  The  models  are  now  extended  to  provide  precise  designs  for  mi- 
crostrip  corporate  feeds.  This  will  lead  to  the  design  of  two  dimensional  transversely  fed 
electromagnetically  coupled  dipole  arrays. 


N.G.  Alexopoulos 
Professor 
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Design  of  Transversely  Fed  EMC  Microstrip 
Dipole  Arrays  Including  Mutual  Coupling 

HUNG-YU  YANG,  member,  ieee,  NICOLAOS  G.  ALEXOPOULOS,  fellow,  ieee,  PHILIPPE  M.  LEPELTIER,  and 

GEORGE  J.  STERN,  senior  member,  ieee 


Abstract— Dts^  IcdmiqMS  aad  proccdarcs  for  oucroslrip  dipole 
arrays  Iraasvcrscly  fed  by  proximity  coopted  nucrostrip  Koes  are  pre¬ 
seated.  Two  desiga  eqaatioas,  which  iadode  the  effects  of  mataal 
coapKag,  an  developed  aad  the  concspoadiog  desigo  carves  arc  obtaiaed 
by  a  rigoroos  iatcgral  eqaatioo  solatioa.  A  sevca-dcaMot  staadiag  wave 
liacar  anay  is  dcsigacd  to  illostratc  the  developed  dcsiga  proccdarcs.  The 
dcsiga  data  an  checked  by  a  compldc  iatcgral  cqaatioa  solatioa  of  the 
array  with  good  agrecawnt.  The  mcasareaicots  of  radiatioo  pattcra  aad 
faipat  inpc^Bce  an  fouad  ia  good  agrecaieot  with  the  dcsiga  goal. 

I.  Introduction 

OLTMAN  INTRODUCED  a  class  of  electromagnetically 
coupled  (EMC)  dipole  antennas  [1],  [2].  The  advantages 
of  EMC  dipoles  are  greater  bandwidth,  higher  efficiency  and 
an  easier  match  to  the  feed  lines,  when  compared  to  classically 
fed  printed  antennas.  Based  on  the  transmission  line  circuit 
model,  Oltman  and  Huebner  [2]  built  an  EMC  dipole  array 
with  radiating  elements  parallel  (collinear)  to  the  feedlines. 
Later,  Elliott  and  Stem  developed  a  Hgorous  design  theory  to 
include  the  effects  of  mutual  coupling  which  successfolly 
predicted  the  array  performance  [3],  [4].  An  efficient  way  to 
obtain  the  design  curves  was  reported  later  [S]-[7]  based  on 
solving  a  Pocklington  type  integral  equation  using  the  method 
of  moments.  The  EMC  collinear  dipole  is  ideally  suited  to  a 
corporate  feed,  and  elements  of  this  type  can  be  arranged  in 
circular  as  well  as  rectangular  grids. 

Another  dipole  antenna  of  the  Oltman  type  is  the  EMC 
transverse  dipole,  as  shown  in  Fig.  1,  where  a  dipole  is 
oriented  transverse  to  an  embedded  microstripiine.  A  string  of 
these  dipoles  above  a  common  microstripiine  becomes  a 
linear  array.  Depending  on  dipole  spacing,  one  can  obtain 
standing  wave  arrays  or  traveling  wave  arrays.  A  family  of 
these  linear  arrays  becomes  a  planar  array.  Current  excitation 
on  the  dipole  is  governed  by  the  amount  of  offset  and  the 
dipole  length.  If  the  dipole  straddles  symmetrically,  no 
excitation  of  the  dipole  occurs.  Weak  coupling  from  line  to 
dipole  can  be  achieved  through  a  slight  lateral  displacement  of 

Manuscript  received  May  5,  1988;  revised  March  IS,  1989.  This  work  was 
supported  by  U.S.  Army  Research  Office  Grant  DAAL03-86-K-0090. 

H.  Y.  Yang  was  with  the  Electrical  Engineering  Department,  University  of 
California,  Los  Angeles,  CA.  He  is  now  with  Phraxos  Research  and 
Development,  Inc.,  Santa  Monica,  CA  9040S. 

N.  G.  Alexopoulos  is  with  the  Electrical  Engineering  Department, 
University  of  California,  Los  Angeles,  CA  90U'4. 

P.  Lepeltier  is  with  Institut  National  Des  Sciences  Appliquees  De  Rennes, 
35043  Rennes,  Cedex,  France. 

G.  J.  Stem  is  with  the  Missile  Systems  Group,  Hughes  Aircraft  Company, 
Canoga  Park,  CA  91304. 

IEEE  Log  Number  8931320. 


El 

H 

^2 

; _ 

Fig.  1.  The  electromagnetically  coupled  transverse  dipole. 


the  dipole  from  its  centered  position.  Because  of  this  weak 
coupling,  large  arrays  are  feasible  [8]. 

Both  theoretical  and  experimental  studies  of  the  EMC 
transverse  dipole  have  been  reported  recently  [9],  [10],  In  this 
paper,  a  design  technique  which  includes  mutual  coupling  is 
developed  for  the  EMC  transverse  dipole  arrays.  Two  design 
equations  will  be  introduced  in  Section  II.  The  methods  for 
generating  design  curves  will  be  discussed  in  Section  III.  A 
design  example  will  be  given  in  Section  IV  together  with  the 
experimental  results.  A  numerical  verification  of  the  design 
obtained  by  solving  the  boundary  value  problem  of  the  whole 
array  system  will  also  be  provided. 

II.  Two  Design  Equations 

For  the  transverse  EMC  dipole  under  consideration,  it  can 
be  shown  from  image  theory  that  the  scattering  off  the  dipole 
is  symmetric.  In  other  words,  the  forward  and  backward 
scattering  coefficients  are  the  same.  Therefore,  in  terms  of  the 
transmission  line  equivalent  circuit,  the  dipole  can  be  approxi¬ 
mated  as  a  shunt  element  with  respect  to  the  feed  [llj.  Each 
dipole  in  the  array  environment  can  then  be  modeled  as  a  two- 
port  network  as  shown  in  Fig.  2,  and  the  whole  system  is  a 
linear  bilateral  network.  Therefore,  one  can  write 

/«=2  V„Y„„  (1) 

ma  1 

as  a.  set  of  equations  connecting  the  transmission  line  mode 
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Fig.  2.  Equivalent  circuit  of  a  dipole  seen  the  transverse  mkrostripline. 

voltages  and  cuirents  at  each  reference  port.  The  active 
admittance  of  each  dipole  seen  by  the  feed  line  can  be  defined 
from  (1)  by 


=  Y„„+  (2) 

where 

m-l  \ 

with  the  prime  on  S  indicating  the  term  n  =  mis  excluded. 
Y„„  is  the  self-admittance  and  is  referred  to  as  the  mutual 
admittance  due  to  the  mutual  coupling  between  each  dipole. 
Since  a  linear  system  is  assumed,  the  current  in  each  dipole 
can  be  written  as 


K  K 


(4) 


where  is  the  current  of  the  nth  dipole  for  a  given  mode 
voltage  in  the  absence  of  other  dipoles  and  l„„  is  the  current  of 
the  nth  dipole  due  to  the  current  in  the  mth  dipole.  The  isolated 
dipole  current  /„„  is  a  function  of  mode  voltage,  dipole  length, 
and  dipole  offset  which  can  further  be  expressed  as 


The  fundamental  design  problem  is  now  obvious.  For  a 
given  design  goal  (radiating  current  in  each  dipole  /^)  one 
wishes  to  find  a  set  of  (s„,  l„)  such  that  not  only  is  (6)  satisfied, 
but  also  the  active  admittance  seen  by  each  feed  line  is  what 
was  prescribed.  The  definition  of  the  radiating  current  in  each 
dipole  may  depend  on  the  design  goal.  For  example,  if  one 
wishes  to  design  a  specific  pattern  in  the  //-plane  (perpendicu¬ 
lar  to  the  dipole),  the  should  be  defined  as  the  current 
collapsed  in  the  feed  line.  In  other  words, 

If=^yAy)dy.  (8) 

In  the  two  design  equations  ((6)  and  (7)),  Y^,f„,  and 
Y„„  can  be  determined  by  the  method  of  moments.  This  issue 
will  be  discussed  in  the  next  section. 

Suppose  that  all  four  functions  are  known.  Further  computa¬ 
tions  are  still  required  to  find  the  dipole  lengths  and  offsets. 
Since  only  relative  currents  in  the  dipoles  are  meaningful,  one 
can  aibitrarily  choose  a  dipole,  say  Ae  mh,  with  length  l„  and 
offset  s„.  For  the  moment,  assume  that  no  mutual  coupling 
exists  and  that  the  left-hand  side  of  (6)  can  be  determined 
according  to  the  design  goal  (desired  currents  in  the  dipoles). 
One  now  can  use  the  first  design  equation  (6)  to  find  N-l 
dipole  lengths  and  offsets.  This  procedure  requires  that  a  two- . 
variable  nonlinear  equation  be  solved  N-l  times.  To  avoid  the 
stability  and  solvability  problem.of  this  nonlinear  equation,  the 
conjugate  gradient  method  [12]  can  be  used  to  provide 
optimized  solutions.  Even  if  the  above  procedures  are  com¬ 
pleted,  a  few  iterations  of  changing  the  dipole  length  or  offset 
of  the  first  selected  dipole  are  required  to  provide  the 
prescribed  input  impedance.  Now  the  design  data  is  what  one 
should  obtain  if  no  mutual  coupling  exists.  To  include  the 
effect  of  mutual  coupling,  one  can  use  the  present  design  data 
to  compute  the  mutual  admittance  and  mutual  currents  and  go 
back  to  the  two  design  equations  repeating  the  above  itera¬ 
tions.  The  whole  procedure  is  iterated  until  convergence  of  the 
design  data  is  found. 


where /„(s„,  l„)  is  a  coefficient  function  relating  the  isolated 
dipole  current  to  its  self-admittance. 

The  two  design  equations  can  be  summarized  from  the 
above  derivations  by 


/rad  y  A/  / 

T  ^  'y' 

•'n  Jn  m.i 


(6) 


and 


N  V 

y  ®  =  y  +  V '  — ^  y 


"  V 

m** !  ^ 


It  is  noted  that,  for  the  EMC  transverse  dipole,  the  current 
phase  variation  in  an  isolated  dipole  is  quite  large,  typically  5° 
to  10°;  while,  in  contrast,  the  current  due  to  mutual  coupling 
has  a  small  phase  variation.  Therefore,  in  the  design,  the 
coefficient  function  is  not  suitable  for  relating  mutual  current 
(/„„)  and  mutual  admittance  ( Y„„).  One  can  use  the  mutual 
current  term  (/„„)  directly  in  (6). 


III.  Discussion  of  the  Method  of  Moments  Solution 

In  order  to  make  an  accurate  design  possible,  information 
about  the  interaction  between  dipoles  as  well  as  dipole 
coupling  to  the  feed  line  is  required.  The  method  of  moments 
provides  a  rigorous  and  accurate  solution.  Integral  equations 
for  the  EMC  transverse  dipole  can  be  written  as 

£.(x,>-)=JJ  G^4‘Hx',r)cfx'  dy' 

+  \\G,yJ<j^Kx\y')dx'  dy'  (9) 

and 

Gy,J^^\x' ,  y')  dx'  dy' 

+  \\GyyJf)ix’,y’)dx'  dy'  (10) 

where  y')  is  the  current  in  the  microstripline  and 

JfKx',  y')  is  the  current  in  the  dipole.  The  functions  Gxx, 
Gxy,  Gyx  and  Gyy  are  the  dyadic  Green’s  function  components. 
A  nice  way  of  modeling  the  feed  line  is  to  use  a  finite  but  long 
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microstrqiline  with  a  fi-gsp  generator  [riaced  &r  from  the  line- 
dqwle  coupling  region  [S],  [6].  When  the  combination  of 
piecewise  sinusoidal  and  Maxwell  current  basis  fruictions  is 
used -in  the  method  of  moments  [S],  [6]  followed  by  the 
Galerldn  procedure,  the  matrix  equations 

IZ„,][/„1  =  [£„1  (11) 

can  be  obtained.  The  excitation  column  vector  E  has  compo¬ 
nents  =  - 1  when  the  6  gap  source  is  located  at  the  center 
of  die  Xrth  basis  function,  and  £„  =  0  anywhere  else.  The 
impedance  matrix  elements  are  in  the  form  of 

!»/2  foa 

•  Xj,)  cos  (XxAjc)  cos  (Xj,Ay)  dX  d4>  (12) 

when  the  mth  and  the  /tth  basis  functions  are  both  on  the  dipole 
or  on  the  line,  and 

!»/2  foo 

XA[/(X)-A(X)]J„(X^,  X,) 

0  *0 

•  ^(X;f,  \y)  sin  (XjfAjf)  sin  (Xj,Ay)  d\  d«f>,  (13) 

otherwise.  The  functions  J„(Xjt,  \)  and  Jm(Xjr,  X^)  are  the 
Fourier  transforms  of  the  expansion  and  testing  functions, 
respectively.  (.Ax,  Ay)  is  the  displacement  vector  of  two  basis 
function  centers,  while  X^  and  X^  are  defined  as 

Xx=Xcos<^  (14) 

and 

Xj,=Xsin0.  (15) 

The  functions /(X)  and  A(X)  are  given  explicitly  in  [13],  [14]. 
The  unknown  currents  in  the  feed  line  or  the  dipole  can  be 
obtained  by  matrix  inversion.  As  a  result,  one  can  use  the 
unimode  transmission  line  theory  to  deduce  the  circuit 
information  from  the  methou  of  moments  solution  for  the 
current  in  the  feed  line.  This  procedure  involves  finding  the 
current  maximum,  minimum,  and  their  positions.  Detailed 
procedures  for  this  are  shown  in  [S].  It  is  noted  that  in  order  to 
find  the  dipole  equivalent  admittance,  the  dipole  can  be  placed 
a  half-eledrical  wavelength  from  the  line  end  such  that  the 
stub  admittance  will  not  be  included  in  the  input  admittance 
observed  by  the  feed  line.  An  example  of  the  input  impedance 
calculation  for  an  isolated  dipole  together  with  the  experimen¬ 
tal  verificatibn  [9]  is  shown  in  Pig.  . 3.  In  the  calculation,  19 
basis  functions  are  used  in  the  dipole  and  100  basis  functions 
are  used  in  the  feed  line  of  1 .5  Xq  long.  The  comparison  shows 
rather  good  agimment. 

One  of  the  features  of  using  a  5  gap  source  is  that  the  mode 
voltage  of  the  line  changes  with  change  of  the  dipole  length  or 
offset  Ooad);  therefore,  care  rhust  be  taken  to  find  this  mode 
voltage,  since  as  shown  in  (6);  dipole  current  is  proportional  to 
the  mode  voltage  and  only  ^eir  ratio  is  useful  in  the  design. 

Another  feature  of  using  the  moment  method  in  the  array 
design  is  that  the-information  about  the  currents  in  the  dipoles 
can  be  obtained  from  the  numerical  process.  This  aspect  is 
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■1 

Fig.  3.  Comparison  between  theoretical  and  experimental  results  for  the 
input  imped^e  of  an  EMC  transverse  dipole.  «i  =  2.17,  <2  =  2.17,  H  = 
1.6  mm,  B  =  0.8  mm,  tvi  =  1  mm,  W2  =  2.2  mm.  Ax  =  0,  As  =  4.5  mm 
and  L  =  12.4  mm. 


Fig.  4.  Equivalent  circuit  of  two  parallel  dipoles  seen  by  a  transverse 
microstripline. 

particularly  helpful,  since  one  can  use  the  dipole  current 
directly  to  design  for  the  desired  excitation  instead  of  using  the 
equivalent  circuit  of  the  dipoles.  This  will  be  discussed  further 
in  the  next  section.  Other  issues  in  this  array  design  are  how 
the  mutual  coupling  information  can  be  separated  from  that  of 
the  self-term  and  how  this  can  be  achieved  without  involving 
the  whole  system  at  the  same  time.  In  order  to  .solve  these 
problems,  certain  assumptions  are  necessary.  It  is  assumed 
that  the  self-admittance  and  self-current  to  mode  voltage  ratio 
willbe  the  same  with  or  without  the  presence  of  other  dipoles, 
and  that  the  mutual  coupling  between  any  two  dipoles  is 
unaffected  by  the  rest.  These  two  assumptions  are  good  jf 
mutual  coupling  is  not  too  strong  [11],  which  is  usually  true 
for  practical  arrays. 

A  method  of  computing  mutual  coupling  of  dipoles  individ¬ 
ually  fed  by  a  microstripline  has  been  discussed  in  [7].  One  of 
the  features  of  the  array  considered  here  is  that  the  dipoles  are 
series  fed  by  transverse  microstriplines.  The  computation  of 
mutual  coupling  in  this  case  requires  a  different  approach  from 
[7].  To  find  the  mutual  coupling  information,  one  can  consider 
two  dipoles  fed  by  a  microstripline  and  follow  a  numerical 
method  similar  to  that  for  an  isolated  dipole  except  for  the 
additional  computation  of  dipole  to  dipole  reaction.  The 
equivalent  circuit  for  these  two  dipoles  is  shown  in  Fig.  4.  An 
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Fig.  S.  A  seven-clement  linear  standing  wave  array.  Elements  are  separated 
by  a  guide  wavelength. 


Real  pait  of  self  admittance 

Fig.  6.  Self-admittance  of  the  dipole  versus  offset  s  and  length  L.s  =  0.32 
+  0.12 (*  -  Dmm.A:  =  1,2,3,  9.  Z,  =  11.91  +  0.025  (n  -  1) 

mm,  n  =  1,  2,  3,  •••  7. 


asymptotic  extraction  technique  together  with  point  source 
approximation  has  been  developed  in  [15]  to  compute  effi¬ 
ciently  and  accurately  the  reaction  of  two  dipoles.  After  the 
matrix  inversion,  the  solution  for  the  current  in  the  line  and 
dipoles  provides  the  total  active  admittance  as  well  as  active 
currents  of  the  dipoles.  For  the  nth  and  mth  dipole  with 
resonant  spacing,  the  method  of  moments  allows  one  to 
compute  the  total  active  admittance 

tot  n  m 

=  Y„„+Y„„^2Y„„.  (16) 


The  minus  or  plus  signs  depqnd  on  whether  the  dipole  spacing 
is  an  odd  or  even  integer  multiple  of  a  half-guide  wavelength. 
From  (16),  if  the  self-adn<ittance  of  each  dipole  is  known,  the 
mutual  admittance  can  be  determined.  The  active  currents  in 
the  dipole  /'  and  /"  can  also  be  obtained  numerically  and  can 
be  described  as 


_  Inn  Imn 

yn~'Vn  Vn 


(17) 


and 


(18) 


From  (17)  and  (18)  together  with  knowledge  of  the  current 
excited  in  an  isolated  dipole,  the  mutual  current  (I„„)  can  be 
obtained. 


Real  pan  of  the  coefficient  function 

Fig.  7.  Coefficient  function  of  the  dipole  versus  offset  s  and  length  L.s  = 
0.32  +  0.12(A:  -  l)mm.A:  =  1,2,3,  •••  9.L--  11.91  +  0.025  (n  - 
1)  mm,  n  =  1,  2,  3,  •  •  •  7. 

IV.  A  Design  Example 

The  previous  discussion  of  the  design  theory  applies  either 
to  a  linear  or  a  planar  array.  Here,  a  seven-element  standing- 
wave  linear  array  will  be  designed  to  illustrate  the  design 
technique.  The  geometry  is  shown  in  Fig.  5.  A  sum  pattern  in 
the  ff-plane  with  a  -  20  dB  sidelobe  level  was  prescribed  for 
this  array.  The  printed  dipoles  are  series  fed  by  a  50  S 
microstripline  embedded  in  the  middle  of  a  substrate  of 
thickness  62  mil  and  permittivity  2.35.  The  design  frequency 
is  8.5  GHz  and  element  spacing  is  chosen  to  be  one  guide 
wavelength.  All  dipoles  have  the  same  width  of  1  mm  and  the 
offsets  and  lengths  are  to  be  found.  The  design  curves  for  the 
self-admittance  Y„„  and  the  coefficient  function  /„  as  a 
function  of  offset  and  length  obtained  from  the  method  of 
moments  solution  are  shown  in  Figs.  6  and  7,  respectively.  It 
is  found  that,  for  the  EMC  transverse  dipole,  many  basis, 
functions  are  required  to  obtain  adequate  convergence.  To 
obtain  each  data  point,  19  expansion  modes  are  used  in  each 
dipole,  and  piecewise  sinusoidal  modes  of  size  of  0.04  guide 
wavelength  are  used  in  the  line.  It  is  observed  from  Fig.  7  that 
for  different  dipole  offsets  and  lengths,  the  phase  of f„  is  not  a 
constant.  This  implies  that  even  for  resonant  spacing,  to  have 
in-phase  excitation,  the  dipole  cannot  be  self-resonant.  To 
obtain  a  perfect  match,  the  stub  length  Ax  in  Fig.  5  can  be 
suitably  adjusted  to  tune  out  the  total  active  susceptance. 

The  sampled  data  are  used  to  construct  the  database  such 
that  for  a  given  offset  and  length  the  function  value  can 
be  obtained  through  a  two-dimensional  interpolation  routine. 
Mutual  coupling  between  two  dipoles  is  a  function  of  dipole 
lengths  and  offsets  for  a  fixed  spacing.  It  is  found  that  mutual 
coupling  is  not  sensitive  to  a  small  change  in  dipole  length. 
Also  from  the  results  for  no  muhial  coupling,  it  is  found  that 
the  lengths  of  all  the  dipoles  are  different  by  less  than  0.2%. 
Therefore,  in  the  mutual  coupling  computation,  the  dipole 
lengths  are  held  fixed.  The  mutual  admittance  and  mutual 
current  as  a  function  of  offsets  for  one  guide  wavelength 
spacing  and  fixed  dipole  length  are  shown  in  Figs.  8  and  9, 
respectively.  A  similar  procedure  can  be  followed  for  a  two 
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Offset  of  the  nih  dipole  (nun) 

Fig.  8.  Mutual  admittance  of  two  dipoles  vetsus  their  offsets.  Element 
spacing  is  one  guide  wavelength  and  £  =  1 1.98  inm. 


Fig.  9.  Mutual  current  of  two  dipoles  versus  their  offsets.  Element  spacing 
is  one  guide  wavelength  and  £  =  1 1.98  mm. 


wavelength  spacing,  and  so  on.  As  a  result,  the  interpolation 
or  extrapolation  method  can  provide  the  mutual  coupling 
information  for  any  dipole  spacing. 

The  design  data  after  a  few  iterations,  including  the  self-  and 
mutual  admittances  are  shown  in  Tabic  I.  The  element  spacing 
is  23.6  mm  and  Ax  -  9.68  mm.  It  is  seen  that  the  mutual 
admittance  is  more  than  25%  of  the  total  active  admittance. 
Theretbre,  it  is  concluded  that  even  at  one  wavelength  spacing, 
the  effect  of  mutual  coupling  should  not  be  ignored.  To  provide 
a  confident  check  of  the  design  data,  the  method  of  moments  is 
applied  to  the  seven-element  linear  array.  The  results  of  active 
admittance  and  current  in  each  dipole  are  shown  in  Table  II 
togedwr  with  the  results  from  the  synthesis  technique.  It  is 
observed  that  the  current  amplitudes  agree  within  1  to  2%  and 
the  phases  agree  within  ±1*.  The  admittance  comparison  is 
also  good. 

The  antenna  array  was  built  on  a  12-in  square  Duroid  board 
as  shown  in  Fig.  10.  The  measured  return  loss  from  the  feed 
line  is  shown  in  Fig.  1 1 .  The  bandwidth  of  this  array  is  about 
4.2%.  At  the  designed  resonant  frequency  (8.5  GHz),  the 
measured  VSWR  is  about  1.1.  The  agreement  between  the 
theoiy  and  measurement  conrirms  the  importance  of  the 


effects  of  mutual  coupling.  Both  the  desired  and  measured 
radimion  patterns  in  the  /f-plane  are  shown  in  Fig.  12.  The 
measured  pattern  is  found  to  be  in  good  agreement  with  the 
design  criteria.  The  asymmetry  of  the  sidelobes  is  due  to  the 
fact  that  the  array  is  not  symmetrically  located  with  respect  to 
the  edges  of  the  finite  array  and  therefore  space  and  surface 
wave  diffraction  at  the  edges  contributes  nonuniformly  to  the 
radiation  pattern.  Another  contributing  factor  may  be  the 
limited  accuracy  of  the  photo-etching  process  in  producing 
identical  microstrip  dipoles  and  highly  accurate  interdipole 
spacing. 

V.  Conclusion 

A  synthesis  method  for  the  design  of  transversely  fed  EMC 
microstrip  dipole  arrays  has  been  presented.  By  using  a 
network  representation,  two  design  equations,  which  include 
the  effect  of  mutual  coupling,  are  developed.  An  iterative 
procedure  using  the  conjugate  gradient  method  has  been 
applied  to  solve  the  design  equations.  The  design  curves  are 
(Stained  numerically  by  the  method  of  moments.  The  method 
of  computing  mutual  coupling  is  also  described.  It  is  found  that 
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TABLE  1 

DESIGN  DATA  FOR  A  SEVEN-ELEMENT  LINEAR' ARRAY 


No. 

Offset  in  mm 

Dipole  Length  in  mm 

ynn 

ymulual 

1 

0.5151 

11.982 

0.0536 -i-j  0.1040 

0.014  -Fj  0.010 

2 

0.6884 

12.000 

0.0805 -Fj  0.0674 

0.038 -Fj  0.023 

3 

0.8925 

11.983 

0.1422 -Fj  0.0093 

0.056 -Fj  0.029 

4 

0.9500 

11.973 

0.1674 -j  0.00% 

0.064 -Fj  0.031 

5 

0.8925 

11.983 

0.1422 -Fj  0.0093 

0.056 -fJ  0.029 

6 

0.6884 

12.000 

0.0805 -Fj  0.0674 

0.038 -fJ  0.023 

7 

0.5151 

11.982 

0.0536 +j  0.1040 

o 

+ 

o 

o 

o 

y*,  =  0.996  +  y0.02  (result  from  iterations) 
Yi„  =  l.(X)0  -  70.04  (result  from  lES) 


TABLE  11 

DESIGN  CHECK  THROUGH  AN  INTEGRAL  EQUATION  SOLUTION 


No. 

D^ired  Current 

lES  Current 

1 

1.0000  <  0.0 

1.0000  <  0.0 

2 

1.2751  <0.0 

1.2%1  <-0.7 

3 

1.6810  <  0.0 

1.6703  <  -0.9 

4 

1.8351  <0.0 

1.8222  <  0.3 

5 

1.6810  <0.0 

1.6730  <-0.7 

6 

1.2751  <  0.0 

12971  <-0.3 

7 

1.0000  <0.0 

0.9940  <  0.7 

The  unit  of  the  phase  of  current  is  in  degrees 


Fig.  10.  A  12-in  square  microstrip  dipole  array. 


to  obtain  a  satisfactory  design,  one  needs  to  include  the  mutual 
admittance  as  well  as  the  mutual  current  in  the  design 
equations.  The  design  of  a  seven-element  standing  wave  linear 
array  is  presented.  The  design  data  obtained  from  the  synthesis 
procedure  are  implemented  in  a  numerical  experiment, 
namely,  solving  the  boundary  value  problem  of  the  whole 
array  system.  The  input  admittance  and  radiating  currents 
from  this  integral  equation  solution  are  found  to  be  in  good 
agreement  with  the  design  goal.  The  measurements  of  the 
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radiation  pattern  and  input  admittance  also  compare  very  well 
with  the  design  criteria. 
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Abstract 


A  rigorous  procedure  is  used  to  analyze  sev¬ 
eral  microstrip  line-to-line  transitions  in  a  shielded 
multi-layer  structure.  The  transitions  studied  in¬ 
clude  edge-coupled  lines,  overlay-coupled  lines  and 
coupled-to-single  lines.  A  power  conservation  check 
based  on  a  rigorous  Poynting  vector  analysis  is  also 
used  to  determine  the  accuracy  of  the  numerical 
convergence.  The  results  of  power  distributions  and 
coupling  coefficients  of  the  line-to-line  transitions 
are  studied  parametrically  to  indentify  the  proper¬ 
ties  and  applications  of  each  transition. 


Introduction 


Proximity-coupled  line-to-line  transitions  are 
important  building  blocks  for  high  frequency  inter¬ 
connects.  Applications  in  millimeter-wave  integrated 
circuits  include  high-pass  filters,  multiplexers  and 
directional  couplers.  Losch  [1]  has  designed  a  broad¬ 
band  highpass  filter  in  realization  of  an  overlay  cou¬ 
pled  line  transition  based  on  a  quasi-static  formula¬ 
tion.  A  more  rigorous  full  wave  analysis  for  coupled 
line  filters  associated  with  the  open  structure  has 
been  discussed  by  Katehi  [2]  for  an  edge-coupled 
transition  and  by  Yang  and  Alexopoulos  [3]  for  an 
overlay-coupled  transition.  In  [3]  a  spectral-domain 
approach  by  expanding  the  current  in  the  coupled 
line  section  with  a  combination  of  entire  domain  and 
subdomain  modes  is  used.  This  mode  expansion 
machanism  seems  to  be  the  most  efficient  and  fruit¬ 
ful  by  far.  For  the  advantage  of  preventing  unnec¬ 
essary  interaction  and  radiation  loss,  a  waveguide 


housing  is  sometimes  more  common  and  practical 
in  the  real  circuit  design.  In  this  work,  a  full-wave 
moment  method  is  used  to  characterize  shielded  mi¬ 
crostrip  line-to-line  transitions. 

Method  of  Moments 


Several  different  types  of  electromagnetically 
coupled  lines  as  shown  in  Figs.  1-3  are  investigated. 
The  methodology  applied  here  is  in  analogy  to  that 
reported  in  [3],  however,  the  spectral  Green’s  func¬ 
tion  and  the  numerical  procedure  arc  very  much 
different.  Since  the  line-offset  and  the  width  of  mi¬ 
crostrips  in  shielded  structures  are  comparable  to 
the  waveguide  dimensions,  the  transverse  current 
component  should  not  be  neglected  and  complete 
dyadic  Green’s  function  of  a  multi-layer  waveguide 
is  required.  The  integral  equation  after  a  Galerkin’s 
procedure  can  be  converted  into  a  set  of  linear  equa¬ 
tions,  when  expressed  in  matrix  form; 


[Z\ut  [ZViL]  (2i2..1  (2{X.] 

IZ2U.\  [ZitU  IZ22..] 


[lu] 


-r 

■  1/n.M  ■ 
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(1) 

where  each  submatrix  [  Z  ],  due  to  the  presence  of 
both  X  and  z  directed  currents,  contains  4  subma¬ 
trices,  for  example: 


[Ziu,]  = 


(2umi,]  [Ziioir] 

[2l],,x>]  [Z\\,izx] 


(2) 


and  submatrix  [  I  ]  and  [  ]  contain  two  subma¬ 

trices  ns  follows: 
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Each  element  in  these  submatrices  represents  the 
reaction  of  different  basis  functions.  For  instance, 
the  elements  of  are  the  reaction  between 

X'directed  currents  of  subdomain  mode  associated 
with  microstrip  1  and  z-directed  currents  of  subdo¬ 
main  mode  associated  with  microstrip  2.  These  sub- 
domain  modes  are  basically  either  PWS  functions  or 
pulse  functions.  The  entire  domain  modes  are  com¬ 
posed  of  the  reflected  mode,  transmitted  mode,  and 
incident  mode  which  are  distinguished  with  abbre¬ 
viation  re /,  Ira,  and  inc,  respectively.  The  computa¬ 
tion  of  each  element  requires  both  infinite  summa¬ 
tion  and  integration  in  spectral  domain  and  their 
expressions  are  in  the  general  form  of: 


nx-oe'-* 

(5) 

and 

n«»«r 

(6) 

where  is  the  spectral-domain  dyadic  Green's 

function.  Superscript  p  identifies  different  entire  do¬ 
main  modes.  Jm  is  the  Fourier  transform  of  trans¬ 
verse  dependence  and  Jm  are  Fourier  transform 
of  iongitudinai  dependence  with  respect  to  subdo¬ 
main  and  entire  domain  modes,  respectively. 


Power  Conservation  Check 

For  shielded  microstrip  transitions,  the  con¬ 
vergence  of  the  moment  method  soiution  is  very  sen¬ 
sitive  to  the  type  and  number  of  expansion  functions 
ciiosen.  Power  conservation  provides  a  nice  way  of 
checking  the  accuracy  of  the  solution.  According  to 
the  configurations  shown  in  Figs.  1-3,  the  incident 
power  should  be  equal  to  the  summation  of  reflected 
power,  transmitted  power  and  some  loss  coupled  to 
the  multi-layered  waveguide  modes.  With  proper 
waveguide  dimension,  the  loss  coupled  to  the  waveg¬ 
uide  modes  can  be  removed  and  the  expression  of 
power  conservation  can  be  simply  written  as 

irp  +  =  I 


where  Zd  and  Za  ere  the  characteristic  impedance 
of  feed  line  and  parasitic  coupled  line  respectively. 
A  frequency-dependent  method  of  computing  the 
characteristic  impedance  of  both  single  and  coupled 
microstrip  transmission  lines  described  in  [4,  5]  can 
be  used  to  determine  the  accuracy  of  the  numeri¬ 
cal  results  given  in  this  work.  In  the  present  com¬ 
putations  of  the  transition  problems,  entire  domain 
modes  of  3  guided  wavelength  and  0  to  18  subdo¬ 
main  modes  are  used  in  each  microstrip  line.  The 
convergence  has  been  checked  within  1%  accuracy. 
The  power  conservation  is  also  checked  with  good 
consistency.  An  example  of  this  check  is  shown  in 
Table  I. 

Numerical  Results  and  Discussions 

From  Figs.  5  -  11,  the  maximum  coupling 
occurs  around  ovl  =  lA,,  (A,,  is  the  wavelength  of 
the  odd-mode  guided  in  the  coupled  line  section) 
with  a  wide  frequency-insensitive  range.  This  im¬ 
plies  that  the  transitions  are  broad-band  and  are 
very  useful  in  many  MMIC  applications.  Besides, 
from  Figs.  5-8,  the  coupling  efficiency  is  better  in 
overlay  line-to-line  transition  than  in  edge  coupled 
line-to-line  transition.  This  indicates  the  former  will 
be  a  promising  element  in  realization  of  millimeter 
wave  high-pass  Alter. 

Figs.  9—10  show  the  results  for  the  case  of 
an  overlay  coupled-to-single  microstrip  transition. 
It  is  seen  that  the  even-mode  coupling  depends  less 
on  the  line-oAset  of  parasitic  coupled  line,  at  com¬ 
pared  with  single-line  coupling.  It  is  also  noted  that 
the  even-mode  of  coupled  lines  can  couple  energy  to 
a  centered  parasitic  microstrip  line  while  the  odd¬ 
mode  can  not.  This  may  find  applications  in  a  phase 
detector. 

The  frequency  response  of  an  overlay-coupled 
microstrip  transition  is  shown  in  Figs.  11  —  13.  The 
geometrical  parameters  are  specially  designed  in  the 
coupled  section  where  the  line-width  is  much  larger 
than  the  spacing  between  two  lines.  It  is  seen  that, 
in  a  wide  frequency  range,  the  coupling  coefficient  is 
almost  independent  of  frequency.  In  addition,  it  is 
possible  to  couple  more  than  95%  of  the  total  power 
through  the  discontinuity.  This  geometry  (so  calletl 


(7) 


Fig.  9  |r|  versus  overlap  for  even-mode  excitation, 

(configuration  Fig.  2.3) 
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Fig.  1 1  Msgnttude  of  T  and  T  versus  frequency. 
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The  topic  of  this  dissertation  involves  a  dynamic  modeling  of  frequency  de¬ 
pendent  passive  integrated  circuit  components.  The  first  part  of  this  disserta¬ 
tion  deals  with  the  characterization  of  microwave  and  millimeter  wave 
discontinuities.  The  analytic  approach  is  based  on  solving  integral  equations 
with  the  method  of  moments.  In  the  method  of  moments  procedure,  an  exact 
Green's  function  is  u.scd,  which  takes  into  account  all  the  physical  effects  in¬ 
cluding  radiation,  surface  waves  and  high  order  mode  effects.  In  order  to 
characterize  the  discontinuities,  a  numerical  scheme  is  developed  to  model 
coupled  semi-infinite  lines.  This  involves  using  both  entire  domain  modes  and 
subdomain  modes  as  the  expansion  functions  to  represent  the  fundamental 
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propagating  mode  and  higher  order  evanescent  rhodcs  respectively.  This  type 
of  expansion  mechanism  is  numerically  efficient  and  particularly  useful  in 
coupled  line  modeling. 

Microstrip  open-end  and  gap  discontinuities  in  layered  integrated  circuits 
arc  studied  with  emphasis  on  the  material,  radiation  and  surface  tvave  effects. 
The  analysis  is  modified  to  study  slotline  short-end  discontinuities  in  sand- 
wichco  structures.  The  advantage  of  using  a  slot  line  sandwich  is  also  dis¬ 
cussed.  With  the  understanding  of  the  simplest  discontinuities  in  microstrip 
and  slotlinc  circuits,  the  study  is  further  extended  to  the  linc-to-line  transition 
which  becomes  increasingly  important  in  monolithic  circuits.  Three  types  of 
circuit  transitions  arc  invc.stigatcd:  microstrip-slotlinc  transition,  proximity 
coupled  collinear  and  transverse  microstrip-microstrip  transition.  The  appli¬ 
cations  of  these  transitions  arc  described  and  the  validity  of  the  developed 
theory  is  verified  by  experiments. 

The  second  part  of  this  dissertation  involves  the  study  of  printed  circuit 
antennas  which  are  special  cases  of  integrated  circuit  discontinuities.  The  nu¬ 
merical  method  for  integrated  circuit  discontinuities  can  be  applied  directly  to 
the  printed  circuit  antenna  structures.  Four  potential  printed  antenna  archi¬ 
tectures  for  millimeter  wave  monolithic  phased  array  applications  arc  dis¬ 
cussed.  These  include  a  microstrip  fed  slot  antenna,  a  slotline  fed  dipole 
antenna,  an  EMC  transverse  dipole  antenna  and  a  microstrip  fed  slot  coupled 
dipole  antenna.  The  features  of  each  architecture  are  described. 

Design  techniques  and  procedures  for  microstrip  dipole  arrays  transversely 
fed  by  proximity  coupled  microstrip  lines  are  also  presented.  Two  design 
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equations  which  include  the  effects  of  mutual  coupling  are  developed.  The 
corresponding  design  curves  are  obtained  by  a  rigorous  integral  equation  sol¬ 
ution.  A  seven  element  standing  wave  linear  array  is  designed  to  illustrate  the 
developed  procedures.  The  data  is  checked  by  a  complete  integral  equation 
solution  of  the  array  with  excellent  agreement.  The  antenna  pattern  and  input 
impedance  are  also  compared  with  the  measurement  with  good  agreement. 
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Chapter  I 


Introduction 


1.1.  Integrated  Circuit  Modeling 

Integrated  circuit  structures  have  been  used  for  numerous  electronics  ap¬ 
plications  due  to  the  features  of  low  cost,  light  weight,  easy  fabrication  and 
reproducibility  [l],  [2].  The  basic  structure  is  composed  of  a  grounded 
substrate  with  various  circuit  components  integrated  on  the  air-dielectric 
interface.  In  solid  state  electronics,  bipolar  or  field  effect  transistors  usually 
require  semiconductor  layers  on  top  of  the  substrate  [3]  .  This  type  of  ge¬ 
ometry  is  usually  referred  to  as  a  layered  structure.  For  frequencies  below  1 
GHz,  passive  components,  which  include  signal  transmission  path  (wire), 
inductance,  capacitance  etc,  have  been  well  understood  using  lumped  circuit 
concepts.  However,  in  the  current  state  of  the  art,  various  applications  require 
the  operating  frequency  of  the  dc\’icc  to  go  beyond  the  microwave  wave  or 
millimeter  wave  range  [4]  -  [7].  Microwave  and  millimeter  wave  integrated 
circuits,  in  the  current  trend,  have  become  increasingly  important  for  both 
commercial  and  military  applications.  When  the  operating  frequency  is  high 
enough  such  that  the  device  dimension  is  comparable  to  a  wavelength,  the 
circuit  components  are  clectromagnelically  coupled  (EMC)  together  and  the 
design  concept  needs  to  be  completely  reformulated.  In  the  past  decade,  nu¬ 
merous  works  have  contributed  to  the  investigation  of  the  electromagnetic 
phenomena  of  microwave  integrated  circuits  (MlCs)  and  microwave 
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monolithic  integrated  circuits  (MMICs).  However,  most  of  the  published- 
works  arc  confined  to  the  characteristics  of  wa\’c  guiding  structures.  For  this 
reason,  many  discontinuity  problems  are  not  well  understood  and  remain  to 
be  solved  [7]-  [14].  From  the  design  point  of  view,  with  the  increasing 
complexity  of  integrated  circuits,  accurate  computer  aided  design  (CAD)  tools 
arc  needed  to  determine  precisely  the  circuit  performance.  Therefore,  currently 
and  in  the  near  future,  developing  accurate  and  efficient  CAD  algorithms 
should  be  a  primary  concern  for  the  microwave  society  [15].  The  early  re¬ 
searches  on  discontinuities  problems  were  mainly  based  on  the  electrostatic 
approximation  due  to  its  simplicity  [16]  —  [18]  .  This  appro.ximation  (or  so 
called  quasi-static  approach)  is  known  to  be  valid  in  the  frequency  transition 
from  low  frequencies  to  the  microwave  range. 

A  full  wave  spectral  domain  approach  (SDA)  for  analyzing  MIC  disconti¬ 
nuities  was  propo.sed  in  the  early  !9S0s  by  Jan.sen  [19]  .  Since  this  SDA  was 
originally  propo.sed  for  waveguide  applications,  the  integrated  circuits  arc  a.s- 
sumed  surrounded  by  an  enclosed  housing.  Therefore,  this  approach  has  the 
limitation  that  radiation  and  surface  wave  effects  arc  not  included 
[20]  -  [23].  Besides,  in  order  to  simulate  the  open  structure,  waveguide  di¬ 
mensions  would  usually  be  large  which  cause  numerical  convergence  problems. 
Therefore,  it  is  reasonable  to  conclude  that  the  SDA  is  good  only  for 
waveguide  discontinuity  problems  or  for  open  structures  with  frequencies  be¬ 
low  C  band,  where  the  radiation  and  surface  wave  effects  arc  less  important. 

An  integral  equation  approach  (lEA)  was  first  propo.sed  in  1985  by  Katchi 
and  Ale.xopoulos  to  characterize  MIC  and  MM  1C  discontinuities  [24]  .  This 
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approach  is  complete  and  rigorous  in  the  sense  that  it  simulates  the  physical 
structure  and  includes  all  the  physical  effects,  although  the  applications  of  this 
method  arc  still  at  the  beginning  stage  [25]  -  [27]  . 

The  first  part  of  this  dissertation  will  be  concerned  with  a  rigorous  exam¬ 
ination  of  the  method  of  moments  solution  of  integral  equations,  in  terms  of 
its  accuracy,  efficiency  and  applicability,  to  the  characterization  of  MIC  and 
MM  1C  discontinuities.  Special  numerical  treatises  will  be  developed  to  char¬ 
acterize  a  cla.ss  of  semi-infinite  line  transitions  frequently  encountered  in  inte¬ 
grated  circuits.  The  derivation  of  Green's  function  will  be  provided  in  Chap. 
II.  Microstrip  open  end  and  gap  discontinuities  in  layered  integrated  circuits 
will  be  discussed  in  Chap  Ill.  Numerical  techniques  adopted  in  this  disserta¬ 
tion  will  al.so  be  di.scu.ssed  in  detail  in  Chap.  III.  It  will  be  shown  in  Chap.  Ill 
that  the  developed  numerical  algorithm  is  general  enough  to  apply  to  many 
related  structures.  The  analy.sis  developed  in  Chap.  Ill  is  extended  in  Chap. 
IV  to  study  slotline  short-end  discontinuities  in  sandwiched  structures.  With 
enough  understanding  of  the  .simple.st  discontinuities  in  microstrip  and  .slotline 
circuits,  the  analysis  can  be  further  extended  to  the  investigation  of  the  im¬ 
portant  subject  of  line-to-line  transitions  in  integrated  circuit  structures.  In 
Chap.  V  a  dynamic  model  for  the  microstrip-slotline  transition  is  developed. 
Two  types  of  novel  transitions  arc  analyzed  in  Chap.  VI,  namely  the  EMC 
collincar  microstrip-microstrip  transition  and  the  EMC  transverse  microstrip- 
microstrip  transition.  The  potential  applications  of  thc.se  transitions  arc  also 
dc.scribed.  The  analysis  is  also  confirmed  by  experiments. 
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1.2.  Printed  Antenna  Modeling 

Integrated  or  printed  circuit  antennas  (PCA)  are  a  natural  evolution  of  in¬ 
tegrated  circuit  components  [28],  [29].  The  radiating  elements  may  consist 
of  dipoles,  slots  or  patches.  Since  the  invention  of  microstrip  antennas  baek 
to  the  early  50s  [30].  extensive  rescarehes  ha\e  been  earned  out  worldwide. 
Mathematical  modeling  was  initially  proposed  in  the  early  1970s  with  a  trans¬ 
mission  line  analogy  [31  ],  [32]  .  Later  on,  a  cavity  method  and  a  quasi-static 
approach  were  also  applied  to  certain  structures  with  some  suceess 
[33]  -  [37]  .  However,  for  frequencies  above  the  mierowave  range,  the 
above  mentioned  methods  will  not  be  accurate  due  to  their  thin  substrate  ap¬ 
proximation.  Besides,  the.se  methods  are  not  applicable  for  printed  dipoles  and 
slots,  and  do  not  include  surface  wave  effects.  A  rigorous  and  almost  exact 
dynamic  analysis  of  PCA  was  finst  performed  at  UCLA  in  1979  [38].  This 
approach  is  based  on  soh  ing  the  Pocklington  type  integral  equation  with  the 
method  of  moments.  Since  then,  an  extensive  amount  of  work  has  been  car¬ 
ried  out  toward  the  efficient  and  accurate  evaluation  of  this  Sommcrfeld-type 
integral  [39]  -  [44]  .  At  the  same  time,  fundamental  properties  of  PCAs 
such  as  gain,  efficiency,  bandwidth  and  antenna  patterns  have  also  been  in¬ 
vestigated  in  terms  of  various  device  parameters  [45]  -  [55].  Fundamental 
radiation  characteristics  of  a  printed  slot  will  be  discussed  in  Chap.  VII.  The 
properties  of  a  printed  dipole  and  a  printed  slot  arc  also  compared  in  terms  of 
the  radiation  efficiency,  surface  wave  effects  and  radiation  resistance. 

One  of  the  crucial  parts  in  the  design  of  printed  circuit  antennas  is  the  de¬ 
sign  of  feeding  structures.  Microstrip  lines  clectromagnctically  coupled  (EMC) 


4 


to  the  radiating  elements  has  proven  to  be  the  most  versatile  and  convenient 
type  of  feed  for  antenna  arrays  [563.  Rigorous  analysis  of  EMC  collinear 
dipoles  has  been  thoroughly  investigated  in  Katehi's  and  Jackson's  disserta¬ 
tions  [57],  [58].  In  the  second  part  of  this  dissertation,  several  new  feeding 
structures  will  be  in\estigated  by  a  method  of  moments  solution  of  coupled 
integral  equations.  The  current  trend  in  printed  circuit  antenna  technology  is 
to  develop  an  architecture  which  is  completely  monolithic.  In  Chap.  VIII, 
four  potential  architectures  for  millimeter  monolithic  phased  array  applications 
are  discussed.  The  features  of  each  architecture  are  decribed.  These  four 
feeding  structures  arc  a  microstrip  fed  slot  antenna,  a  slotlinc  fed  dipole  an¬ 
tenna,  an  EMC  transverse  dipole  antenna  and  a  microstrip  fed  slot  coupled 
dipole  antenna. 

The  ultimate  goal  in  the  analysis  of  printed  circuit  antennas  is  to  provide  a 
design  technique  of  microstrip  arrays.  Looking  back  to  the  past,  the  only  rig¬ 
orous  and  complete  design  of  microstrip  arrays  was  published  by  Elliott  and 
Stern  in  1981  [59],  [60]  .  In  Chap.  IX  of  this  dissertation,  a  new  type  of 
microstrip  array  will  be  investigated.  The  geometry  is  composed  of  dipoles 
parallel  to  each  other  fed  by  an  EMC  transverse  microstriplinc  embedded  in 
the  substrate.  The  design  of  an  antenna  array  requires  the  knowledge  of  each 
antenna  element's  size  and  position  with  respect  to  other  elements  and  the  feed 
such  that  the  desired  excitation  in  each  antenna  and  therefore  the  desired  an¬ 
tenna  pattern  can  be  achieved.  Two  design  equations  which  arc  closely  related 
to  those  developed  by  R.S.  Elliott  in  waveguide  slot  array  design  [61  ]  will  be 
discussed.  The  corresponding  design  curves  will  be  obtained  numerically 
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through  a  rigorous  integral  equation  solution.  In  order  to  show  how  the  design 
technique  works,  a  se\en  element  standing  wave  linear  array  will  be  designed 
to  have  the  performance  of  20  dB  Dolph-Chebychev  broadside  sum  pattern  in 
the  H-  plane.  The  importance  of  the  effects  of  mutual  coupling  will  also  be 
addressed. 
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Chapter  11 


Analysis  of  Hertzian  Dipoles  in  a  Layered  Planar  Geometry 

In  this  chapter,  the  electromagnetic  fields  due  to  an  electric  Hertzian 
(infinitesimally  small)  or  a  magnetic  Hertzian  dipole  arc  derived.  These  fields, 
or  so  called  Green's  function  must  be  known  before  the  method  of  moments 
can  be  applied.  This  Green's  function  for  a  layered  planar  geometry  is  well 
known  as  a  Sommcrfeld-typc  solution  [62].  In  section  2.1,  the  Green's  func¬ 
tion  for  an  electric  dipole  is  discussed.  In  section  2.2,  the  Green's  function  for 
a  magnetic  dipole  will  be  developed. 

2.1,  Green's  Function  for  an  Electric  Hertzian  Dipole. 

For  the  substratc-superstrate  geometry  under  consideration,  the  formu¬ 
lation  of  the  Green's  function  for  a  dipole  embedded  in  the  superstrate  (Fig. 
2-1)  or  in  the  substrate  (Fig.  2-2)  is  different  and  will  be  discussed  separately. 
The  ca.se  for  a  dipole  in  the  supenstrate  will  be  considered  first. 

(a).  An  Electric  Hertzian  Dipole  in  the  Superstrate 

The  geometry  shown  in  Fig.  2-1  contains  a  substrate  with  permittivity  e,  , 
permeability  /t,  and  thickne.ss  b,  and  a  superstrate  with  permittivity  e,  and 
permeability  The  total  thicknc.ss  of  the  material  is  h.  The  dipole  is  located 
at  z  =  Zj  with  b  <  z,  <  b  and  oriented  in  the  x  direction. 


7 


A 

-►  X 


ground  plane 


Figure  2-1.  An  electric 
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Figure  2-2.  An  electric  Hertzian  dipole  in  the  substrate. 
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The  electromagnetic  fields  in  each  region  can  be  characterized  in  terms  of  the 
Hertz  potential  in  the  following  forms: 


(2.1) 


and 


=  yw£,eQ(Vx^'^) 


(2.2) 


with  i  representing  the  region  (i).  From  Sommerfeld's  work  [62],  it  is  known 
that  the  Hertz  potential  may  contain  only  x  and  z  components.  In  other  words, 


(2.3) 


The  nice  feature  of  using  the  Hertz  potential  is  that  its  components  satisfy  the 
scalar  wave  equations,  where  the  solution  can  be  obtained  in  a  straight  for¬ 
ward  manner.  Through  the  use  of  the  Hertz  potential.  Maxwell's  equations 
can  be  reformulated  in  the  form  of 


— -  v,)(5(r  - 
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The  lateral  electromagnetic  fields  in  region  i  arc  related  to  the  Hertz  potential 
through 
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(2.7) 
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With  the  boundary  conditions  that  Ey,  Hj.  and  Hy  are  continuous  across 


each  interface,  and  with  the  radiation  condition  at  z  -»■  oo,  one  is  able  to  obtain 
an  expression  for  the  Hertz  potential  in  the  Fourier  domain  from  eight 
equations  with  eight  unknowns.  The  resulting  Hertz  potential  in  each  region 
can  be  summarized  as: 

In  region  0:  h  <  z 
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In  region  I:  0  <  2  <  /» 


(2.11) 
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In  region  2:  B  <  z  <  H 


(2.13) 


and 
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The  functions  fp.)  and  /;//)  arc  shown  in  Appendix  A.  Other  parameters  arc 
defined  as  follows: 


op.)  =  qpi  + 
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and 
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+  (r/2/ii  +  - - - )  smh  qpi 

/^2V2 


h) 

h) 


(2.16) 


where 


’  /:2  ,  ,2 

A  -  P  /.^  +  Ay  , 

(2.17) 

q  =  \/a-  -  /C(y, 

(2.1S) 

..  /32  1.2 
<7i  =  v^-  -  ^1  . 

(2.19) 

r/2  =  -  A'l , 

(2.20) 
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^’O  ~  (2.21) 

^’l  ~  ?  (2.22) 

^2  ~  ^0\  ^2^2  '  (2.23) 

fi\  =  ^/^l£l  ,  (2.24) 

/?2  =  (2.25) 

and 

V  =  An~Jo}t:2R(\-  (2.26) 


(b).  An  Electric  Hertzian  Dipole  in  the  Substrate. 

When  the  dipole  is  embedded  in  the  substrate,  as  shown  in  Fig.  2*2,  the 
eleetromagnetic  fields  will  be  different  from  the  ease  when  the  dipole  is  in  the 
superstrate,  although  the  analytic  procedure  is  very  much  the  .same.  The 
dipole  in  this  ease  is  located  at  z  =  z,  with  0  <  Zj  <  b  and  oriented  in  the  x 
direction.  The  electromagnetic  fields  in  each  region  can  be  characterized  in 
terms  of  the  Hertz  potential  given  in  Eqs.  2.1  and  2.2  which  satisfy  the 
equation 

+  kfl^  =  ■  (2-27) 

The  lateral  electromagnetic  fields  in  region  i  are  related  to  the  Hertz  potential 
through  Eqs.  2.10-2.13.  With  the  same  boundary  conditions  as  the  case  for 
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a  dipole  in  the  superstraie,  the  Hertz  potential  in  each  region  can  be  summa¬ 
rized  as: 

In  region  0:  h  <  z 

4‘i’=-|rr  r (2.28) 

y  •'-oo  *'-00 


'  •'-00  ’'-00 


In  region  \:0  <  z  <  h 


(1)  J.p  -j>.p-ys)^^  y; 


(2,30) 


4i'  =  r  Mil n 

I  •’-00  *^-00 


In  region  2:  h  <  z  <  h 


IVA- 


_  p  p"  rp-)  -j^x  -  .V,)  -  >g 


(2.32) 


~  dX  d/  (2  33) 


where 
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(2.34) 


}•'  =  4;r7a)5i6o- 

The  functions /',(/)  and  //'//)  arc  shown  in  Appendix  B.  Other  parameters  are 
found  in  Eqs.  2.15-2.26. 

2.2.  Greenes  Function  for  a  Magnetic  Hertzian  Dipole. 

Slots  in  planar  layered  structures  arc  ahso  considered  in  this  dissertation. 
The  problem  is  formulated  for  the  case  where  slots  in  a  ground  plane  are 
sandwiched  by  two  layers  on  both  sides  of  the  ground  plane.  From  the 
Stratton-Chu  .solution  of  Maxwell's  equations  [611,  one  is  able  to  see  that  the 
tangential  electric  fields  in  the  boundary  can  be  viewed  as  sources  which  gen¬ 
erate  electromagnetic  fields  inside  the  region  of  interest.  These  fictitious 
.sources,  usually  called  magnetic  currents,  when  placed  in  Maxwell's  equations 
arc  dual  to  electric  currents.  Maxwell's  curl  equations,  with  only  magnetic 
current  sources  arc 

V  X  £  =  ~  J,„  -  J(OfiH  (2.35) 

and 

VxH  =  jioeE.  (2.36) 

These  two  equations  arc  constructed  for  the  reason  of  mathematical  conven¬ 
ience.  The  procedure  of  interchanging  sources  and  boundary  conditions  is,  in 
fact,  well  known  in  the  realm  of  partial  differential  equations  [64].  For  the 
geometry  under  consideration,  one  can  argue  that  the  electromagnetic  fields 
above  the  ground  plane  arc  due  to  magnetic  currents  (physically  they  are 
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electric  fields)  at  a  slot  position  infinitesimally  close  to  the  ground  plane.  The 
same  arguments  are  valid  for  the  space  below  the  ground  plane.  In  this  situ¬ 
ation,  the  ground  plane  is  in  effect  closed  without  any  apertures.  Therefore, 
in  the  formulation  of  electromagnetic  fields  in  a  half  space  (cither  above  or 
below  the  ground  plane),  the  electromagnetic  fields  on  the  other  side  are  irrel¬ 
evant.  The  nice  feature  of  this  is  that  in  the  Green's  function  formulation,  one 
can  consider  the  geometry  above  or  below  the  ground  plane  sequentially. 
However,  if  the  magnetic  currents  in  the  slots  arc  to  be  determined  instead  of 
being  known  beforehand,  the  fields  in  the  whole  space  will  be  related  through 
an  additional  boundary  condition  across  the  slots.  This  additional  boundary 
condition  will  lead  to  an  integral  equation  which  can  be  solved  by  the  method 
of  moments. 

The  Green's  function  for  a  magnetic  Hertzian  dipole  is  obtained  by  solving 
Ma.xwell's  equations  with  magnetic  currents  replaced  by  a  delta  source  ori¬ 
ented  in  the  .v  direction.  Since  the  source  is  known,  only  the  half  space  as 
shown  in  Fig.  2-3  needs  to  be  considered.  Due  to  geometric  symmetry,  the 
formulation  for  the  other  half  space  is  the  same.  The  electromagnetic  fields  in 
each  region  can  be  characterized  in  terms  of  the  Hertz  potential  in  the  follow¬ 
ing  forms: 

//'"  =  + (2.37) 

and 

£"*=  (2.38) 
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ground  plane 


Figure  2-3.  A  magnetic  Hertzian  dipole  in  a 

substrate-superstrate  configuration. 
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where  i  represents  the  region  (i)  and 


Eq.  2.40  can  be  solved  subject  to  the  boundary  conditions  at  each  interface  and 
z-^oo  .  The  resulting  Hertz  potentials  in  each  region  are: 

In  region  0:  fi  <  z 
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C  =  (2.45) 


— oo  — oo 


Tr^JJl  =  -L  f ^  jVhrd^-)  f,-4.2-h)  -P4x-  xs)  -j>.y{y-  jg  ^  p 

y^nLj_^D,{A)D,4^k)'  (-^6) 


In  region  1:  0  <  z  <  b 


'mi_^i_^D,p.)  ^  > 


(2.47) 


,ou_Lr  r  !hMli^,-U-->=,)ri>-P-y,)ax  dX  (2^ 
"'‘  dx,dx^.{2. 


In  region  2:  b  <  z  <  li 


« ■  icr.-ra''"-'''-'''''"'".-".  <“> 


where 


y,„  =  47t7w/^,/iQ. 


(2.51) 


The  lunction.sX„(-^>)  and  //„„(2)ean  be  found  in  Appendix  C.  Other  parameters 
arc  defined  in  Eqs.  2.15-2.26. 
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Chapter  111 

MIC  Discontinuities  I  - 
Microstrip  Open>end  and  Gap  Discontinuities 
in  a  Two-layer  Configuration 


3.1.  Introduction 

Microstrip  open-end  and  gap  discontinuities  are  useful  in  the  design  of 
matching  stubs  and  coupled  filters.  In  recent  years,  layered  integrated  struc¬ 
tures  have  found  various  applications  in  MIC  and  printed  circuit  antennas, 
especially  for  monolithic  applications.  Therefore,  design  data  for  open-end 
and  gap  discontinuities  in  layered  structures  would  be  useful.  In  this  chapter, 
microstrip  open-end  and  gap  discontinuities  in  a  substrate-superstratc  config¬ 
uration  (Fig.  3-1  and  3-2)  will  be  considered.  The  characterization  of  these 
discontinuities  for  a  single  layer  has  been  performed  quite  extensively  in  the 
past.  Quasi-static  analysis  based  on  solving  Poisson's  equation  has  been  ap¬ 
plied  for  low  frequency  applications  [16]  -  [18],  [65]  .  For  higher  fre¬ 
quencies,  models  based  on  rigorous  dynamic  analysis  are  required.  A 
spectral-domain  approach  [19],[22]  has  been  used  to  characterize  the  dis¬ 
continuity  problems  with  an  enclosed  waveguide  housing.  A  dynamic  method 
based  on  solving  integral  equations  by  the  method  of  moments  has  recently 
been  applied  to  the  modeling  of  microstrip  open-end  and  gap  discontinuities 
for  a  single  layer  ca.se  [24],  [25].  This  analysis  takes  into  account  all  the 
physical  effects  including  radiation,  surface  waves  and  dominant  as  well  as 
higher  order  mode  coupling.  In  [24],  a  finite  but  long  microstrip  line  with  a 
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Figure  3-1.  An  open-end  microstrip  line  in  a  two-layer 
structure. 


Figure  3-2.  Microslrip  gap  in  a  two-layer 


structure. 


d  gap  source  is  used,  which  typicall}'  requires  many  basis  functions  and  is  nu¬ 
merically  inefficient.  In  [25],  a  more  efficient  method,  using  a  combination 
of  the  entire  domain  modes  and  subdomain  modes,  is  used.  However,  from  the 
discussion  presented  in  [25],  it  seems  that  the  method  docs  not  provide  reli¬ 
able  results  for  the  capacitance  calculation.  In  this  chapter,  a  revised  analysis 
of  [25]  together  with  a  detailed  discussion  of  the  convergence  of  the  solutions 
is  presented.  Since,  for  the  rest  of  this  dissertation,  the  methodology  adopted 
here  will  be  repeated,  the  analysis  will  be  discussed  in  extensive  detail  in  this 
chapter. 

A  crucial  step  (and  difficult)  is  choosing  the  proper  basis  functions  to  pro- 
\'idc  efficient  and  accurate  numerical  computation.  In  the  present  problem, 
modeling  of  semi-infinite  lines  is  required.  A  combination  of  semi-infinite 
traveling  wave  modes  and  local  subdomain  modes  is  fruitful  and  can  be  mod¬ 
ified  easily  to  adapt  to  different  geometries.  The  traveling  wave  mode  corre¬ 
sponds  to  the  fundamental  guided  wave  mode  of  the  microstrip  line.  The  local 
subdomain  modes  are  used  in  the  vicinity  of  the  discontinuity  region  to  take 
into  account  the  higher  order  mode  effects.  For  the  transvcr.se  dependence  of 
the  expansion  functions,  it  is  possible  that  with  sufficiently  high  frequencies, 
the  simple  Maxwellian  or  pulse  function  used  in  [24],  [25]  may  not  be  a  good 
approximation  when  the  dominant  mode  is  not  TEM-like.  Therefore  in  this 
analysis,  the  tran.svcrsc  dependence  of  the  longitudinal  current  is  obtained  by 
a  two  dimensional  infinite  line  analysis  where  three  modified  cosinc- 
Ma.xwcllian  functions  arc  u.scd.  The  characterization  of  an  open-end  disconti¬ 
nuity  is  through  the  open-end  capacitance  which  is  mainly  due  to  the  fringing 
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electric  fields  [17].  When  radiation  and  surface  wave  losses  are  considered, 
a  conductance  should  also  be  included  in  the  equivalent  circuit  model.  The 
gap  discontinuity  can  be  modeled  as  a  ;r  network  with  two  capacitances  (see 
Figs.  3-3  and  3-4).  The  loss  mechanism  can  also  be  included  by  adding  two 
conductances  in  Fig.  3-3.  The  material  effects  of  layered  structures  on  the  ra¬ 
diation  and  surface  wave  lo.ss,  and  the  fringing  fields  at  the  discontinuities  will 
also  be  discussed. 

3.2.  The  Method  of  Moments  and  Matrix  Formulation 

The  transverse  current  has  been  found  in  the  past  to  be  a  few  orders  of 
magnitude  smaller  than  the  longitudinal  current  for  a  strip  width  <  O.Ug  and 
is  neglected  for  simplicity  [8].  Under  such  circumstances,  the  integral 
equation  for  the  open  end  case  can  be  simplified  in  terms  of  the  longitudinal 
electric  field  on  the  microstrip: 

rO  rwj2 

£;^.(.Y,y',z)  =  y,  z\x,,y^,  Zj)J^Xs,ys)ily^({Xs  (3.1) 

where  is  the  electric  field  due  to  the  eurrent  at  ^  =  Zy  The  Green's  func¬ 
tion  G^jf  is  the  value  of  on  the  (.\,  y)  of  microstrip  due  to  an  x  directed  delta 
source  at  (x^y^).  This  Green's  function  has  been  de.scribcd  in  Chap.  II.  In 
the  method  of  moment  procedure,  the  unknown  current  distribution  is  ex¬ 
panded  in  terms  of  a  set  of  known  functions.  An  efficient  way  is  to  use  a 
semi-infinite  traveling  wave  mode  to  represent  the  fundamental  guide  wave 
mode  of  the  microstrip  line  and  to  u.se  subdomain  modes 


Figure  3-3.  Equivalent  circuit  of  a  microstrip  gap. 


Figure  3-4.  Top  view  of  a  microstrip  gap  discontinuity. 
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(piecewise  sinusoidal  modes)  to  describe  the  current  behavior  in  the  vicinity  of 
the  discontinuities.  Since  the  microstrip  open  end  is  a  special  case  of  the  gap 
discontinuity,  the  formulation  in  the  following  will  be  for  the  gap  case.  The 
open-end  case  will  be  discussed  later.  The  current  in  the  microstrip  gap  shown 
in  Fig.  3-4  can  be  e.xpandcd  as 

y,(.v.y)  =  y(.v)y/)-)  (3.2) 


with 

A' 

y(.v)  =  -  r  forxiO,  (3.3) 

rt=l 

A' 

Xv)  =  +  ^/„,ay.v)  /orxas,  (3.4) 

lt=  1 

and 


Jh')  = 


<7(,  +  o,  co.s(-^jO  +  ih  co.s(-^y) 


;nrv  I  -  (2y/w)“ 


(3.5) 


where  F  is  the  retlection  coefficient  from  the  discontinuity  and  T  is  the  wave 
amplitude  of  the  transmitted  wave.  The  layout  of  the  e.xpansion  modes  is 
shown  in  Fig.  3-5.  The  parameters  k„,  a^,  a,  and  ^2  are  obtained  through  an 
infinite  line  analysis,  which  involves  solving  a  characteristic  equation  in  a  ma¬ 
trix  form  [66]  .  The  piecewise  sinusoidal  modes  (PV'S)  are  defined  as 


(3.6) 
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and 


-  for  u-nd,-s\<d, 


where  </,  is  half  length  of  the  PWS  mode.  The  choice  of  can  be  quite  ar¬ 
bitrary. 

When  the  expansion  functions  arc  used  in  Eq.  3.1,  followed  by  Galerkin's 
procedure  [63],  integral  equations  arc  converted  into  a  .set  of  linear  equations. 
These  2N  +  2  equations  when  cxprcs.sed  in  matrix  form  arc 


r  p 


[/me,] 

[/me,] 


The  matrix  elements  in  each  submatrix  can  be  expressed  in  the  following 
form  after  some  tedious  algebraic  manipulations: 

[Z^f/y]  is  an  (A’  +  I)  x  N  matrix  with  matrix  elements 


=  IT  cos[;.>?  -  (3.9) 


—•>0  —CO 


[Z„^y]  is  an  (A'+  1)  x  1  column  vector  with  elements 


rr  (3.10) 


— CO'^— CO 


[2«aa]  +  Ox  A'  matrix  with  matrix  elements 
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7 — 
^ccact  ~ 


—  ->  T 

— — OO 

cos[A  ^.(«rfj  +  null  ~  » 


(3.1I) 


^^leact^  is  an  (A'+  I)  x  I  column  vector  with  elements 

poo  poo  _  ^ 

2,cac,=  \  I 

—00—00 

cos[/i^.(//i/j  -  .?)]  dXy.  dXy , 


(3.12) 


C4u-|l  (/V+  1)  X  1  column  vector  with  elements 


— >C»  — OO 


[/„„,]  is  an  (A'+  1)  x  1  column  vector  with  elements 


G.yJ'-.y.  ■>•,.)  d(’y)Qp.:,)  A, u^) 
CQnlX^Oidi  -  .9)]  dX^  dXy , 


“00—00 


(3.13) 


where 


A  ,(2 , 


1  ^  (cosA^K/l  -  cos2^^/,) 

‘A-l  -  ^ - , 


{Xl  -  kh) 


(3.14) 


j 

^  +  kn)  +  J4r-^Xy  -  ^'Tt)], 

k=() 


(3.15) 
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fO 


(3.16) 


-p-x^-  fO 

2l(4)  =  2Am  - /]  sinA'„,,-V(?^"-^  J.V 

V  ^  _ 


(3.17) 


f  .sin(A„,v)6-M/.v  =  -  ^(4  +  AJ].  (3.18) 


The  function  (7^^^  is  the  Fourier  tran.sform  of  the  Green's  function.  From  Eqs. 
3.9-3.13.  one  can  sec  that  this  Green's  function  will  play  an  important  role  in 
the  numerical  computation.s. 


3.3.  Numerical  Techniques 

Eqs.  3.9-3.13  indicate  two  types  of  integrations.  One  is  related  to  the  PWS 
and  PWS  modes  reaction  while  the  other  is  related  to  the  traveling  wave  and 
the  PWS  modes  reaction.  The  case  of  the  PWS  and  PWS  modes  interaction 
will  be  considered  first.  Eq.  3.9,  after  being  transformed  into  polar  coordi¬ 
nates.  can  be  written  as 


^■j,.) ?(4)  co.s[4(w  -  nyi^']  X  dM({> 


(3.19) 


where 


-  yZo  p  A?  -  xl 
4;r*-A„ 


Da{X)DM) 


(3.20) 
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=  X  coscj), 


(3.21) 


X^  =  X  sin  <j>.  (3.22) 

The  functions  f,  and  can  be  obtained  from  the  Green's  function  discussed  in 
Chapter  11. 

Eq.  3.19  involves  a  double  integration.  For  the  finite  integral,  a  32  points 
Gauss  quadrature  formula  is  used  [57].  For  the  infinite  integral,  speeial  nu¬ 
merical  methods  are  required.  One  can  break  the  infinite  integration  range 
into  two  parts  (0,  A)  and  (A,  oo)  such  that  in  the  first  section  the  integrand 
contains  singularities  or  derivative  singularities,  while  for  the  second  section 
the  integrand  is  well  behaved  but  slowly  convergent.  The  ehoice  of  A  is  quite 
fie.xiblc,  but  it  should  satisfy  A  >  ma.\  (A-,,  A^).  The  first  integral  contains  sur¬ 
face  wa\'e  poles  whenever  Dp.)  or  D,p.)  become  zero.  If  a  pole  e.xtraction 
technique  is  applied  [39]  in  addition  to  the  residue  and  Cauchy  Principal 
value,  four  sections  of  integrations  are  required  due  to  the  derivative 
singularities  at  X  =  Aq,  A,  and  Aj.  Another  way  to  perform  the  integration 
from  0  to  A  is  to  deform  the  contour  off  the  real  a.xis  and  apply  the  Cauchy 
Riemann  theorem  such  that  the  integrand  is  well  behaved  [67].  This  method 
is  particularly  u.seful  in  a  multi-layered  structure,  since  it  is  not  required  to 
know  the  pole  position  and  the  integration  has  no  singularities.  Both  of  the 
above  mentioned  methods  have  been  used  and  a  negligible  difference  has  been 
observed. 

The  second  integration  from  A  to  oo  is  the  so-called  tail  integration.  For 
an  asymptotic  approximation,  it  is  convenient  to  choose  A  such  that 


29 


tanh  /.b  ^  1  and  tanh  ^  1  .  for  X>  A.  If  one  studies  the  asymptotic  behav¬ 
ior  of  the  Green's  function  ?.y)  (omitting  the  term  outside  the  bracket 

of  Eq.  3.20),  one  can  find  that  as  /  -»  cxd,  G^^^  converges  slowly  when  the  test¬ 
ing  and  observation  points  are  on  the  same  plane  (same  z)  and  decays  e.xpo- 
ncntially  otherwise.  This  is  because  the  Green's  function  contains  a  singularity 
in  the  space  domain  due  to  the  delta  source  and  converges  slowly  in  the 
Fourier  domain.  The  dominant  and  first  order  asymptotic  behavior  of  can 
be  summarized  by  the  following: 


Case  (a):  z  =  =  h  =  {b  +  t) 


2c/,.  V/  X^ 


where  ^  and  q,  =  JX^  ~ 


Ca.se  (b):  h  <  z  =  Zy<  h 


where  =  Sj  and  c/,.  =  JX^  -  kjyy . 


(3.23) 


(3.24) 
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(^1  +  Si)  j  — TT~ 

Where  ^ - and  q,  =  -  k]ff. 


Case  (d):  b  <  z  =  z^<  h 


(%- 


2£ 


^ffic 


+ 


(3.26) 


where  z^jj  =  e,  and  -  k^ff . 

It  can  be  seen  that  the  dominant  term  of  the  Green's  function  is  in  the  same 
form  for  all  four  cases,  which,  when  transformed  back  to  real  space  is  due 
to  a  delta  current  in  the  x  direction  in  a  homogeneous  space  with  an  effective 
dielectric  constant  The  dominant  term  of  the  integral  will  be  computed 
separately.  This  formulation  will  be  discussed  in  Appendix  D.  It  is  seen  that 
the  first  order  term  left  in  Eqs.  3.23-3.26  will  enable  the  integrand  in  Eq.  3.19 
to  converge  either  exponentially  or  in  the  order  I/A^  and  can  be  integrated  nu¬ 
merically  in  a  straightforward  way.  It  has  been  shown  in  Appendix  D  that  if 
A\.,  in  the  PWS  modes  is  chosen  to  be  the  same  as  Ayy,  the  computations  arc 
greatly  simplified.  Therefore,  throughout  this  dissertation,  PWS  expansion 
modes  arc  chosen  in  this  way. 

The  other  type  of  integration  is  related  to  the  PWS  modes  and  traveling 
wave  mode  reaction.  When  Eq.  3.18  is  inserted  in  Eq.  3.10,  with  suitable  re¬ 
arrangements,  one  has 


r 

CO 


2k 


m 


~P-x^ 


(3.27) 
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where 


(3.28) 


At  a  first  look,  it  seems  that  Eq.  3.27  contains  poles  at  =  +  k„  which  may 
cause  trouble  in  the  numerical  integration.  However,  as  mentioned  before, 
k,„  is  obtained  by  an  infinite  line  analy.sis  where  the  characteristic  equation  is 
S(A^)  =  0.  Therefore,  one  would  have  S{k„)=  0  and  the  singularities  in  Eq. 
3.27  turn  out  to  be  removable.  The  integration  in  Eq.  3.27  after  being  trans¬ 
formed  into  polar  coordinates  may  be  cxprc.sscd  as 


=  25, ;r/,  -  ^)  + 
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It  is  obvious  from  Eqs.  3.19  and  3.30  that  the  integral  in  Eq.  3.30  can  be 
computed  numerically  in  the  same  manner  as  the  one  occurring  in  Eq.  3.19 
except  for  the  tail  integration.  This  difference  is  due  to  the  fact  that  the  trav¬ 
eling  wave  mode  is  semi-infinitely  long  so  the  technique  used  for  Eq.  3.19  can 
not  be  applied.  In  other  words,  if  the  integration  is  transformed  back  to  real 
space,  the  integration  over  the  half-infinite  microstrip  line  will  cause  other  nu¬ 
merical  convergence  problems.  However,  due  to  the  use  of  entire  domain  base 
functions,  the  integration  in  Eq.  3.30  usually  converges  better  than  the  one  in 
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Eq.  3.19  because  the  integrand  contains  the  Fourier  transform  of  the  entire 
domain  mode  which  decreases  quickly  away  from  the  point  where  is  equal 
to  the  phase  constant  k^.  Further  details  in  the  asymptotic  analysis  of  Eq. 
3,27  by  using  branch-cut  integration  is  shown  in  Appendix  E. 

Although  the  above  discussion  is  for  the  microstrip  gap  case,  the  open-end 
ca.se  is  also  included.  From  the  submatrices  and  C/,>,c,!]  in  Eq. 

3.8,  the  information  of  an  open-end  microstrip  line  can  be  obtained. 

3.4.  Results 

Usually  MIC  discontinuities  are  characterized  by  their  equivalent  circuits. 
Therefore,  in  order  that  the  characterizion  be  meaningful,  all  the  disturbances 
of  the  current  should  die  out  quickly  as  one  moves  away  from  the  discontinu¬ 
ity.  In  microstrip  structures,  any  discontinuity  will  generate  radiating  and 
surface  wa\c.s.  These  waves  will  also  propagate  along  with  the  microstrip  line 
fundamental  mode.  Therefore,  when  radiating  or  surface  waves  arc  strong 
enough  such  that  their  interactions  with  the  microstrip  guided  mode  become 
noticeable,  the  computed  equivalent  circuits  will  not  be  accurate.  This  implies 
that,  for  this  ca.se,  if  one  tries  to  measure  the  equivalent  circuits,  the  results 
will  be  different  at  different  reference  planes. 

The  equivalent  circuit  of  a  microstrip  open  end  is  computed  from  the  re¬ 
flection  coefficient  which  is  obtained  directly  from  the  matrix  inversion  of  Eq. 
3-8.  The  convergence  of  the  results  depends  on  the  size  of  the  c.xpansion 
functions  and  the  region  where  subdomain  modes  arc  u.scd  (.sec  Fig.  3-5).  An 
example  of  a  convergence  check  for  the  open-end  conductance  and  capacitance 
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for  the  number  of  modes  used  in  a  fixed  subdomain  mode  region  is  shown  in 
Figs.  3-6  and  3-7.  The  subdomain  region  is  chosen  as  10  times  the  substrate 
thickness.  It  is  found  that  with  only  two  modes  a  convergent  result  for  the 
conductance  calculation  has  already  occurred.  However,  the  convergence  of 
the  capacitance  value  is  very  slow  with  respect  to  the  number  of  expansion 
modes.  The  region  where  the  subdomain  modes  are  used  is  less  important. 
An  example  of  this  convergence  check  is  shown  in  Fig.  3-8  where  the  basis 
function  size  is  fi.xcd.  It  is  found  that  the  results  are  almost  unchanged  in  a 
wide  range.  Physically,  this  means  that  the  higher  order  modes  generated  by 
the  discontinuity  have  already  died  out  in  the  testing  region.  The  above  con¬ 
vergence  tests  are  for  the  case  that  the  radiating  and  surface  waves  are  weakly 
excited.  It  is  found  that  the  results  do  not  converge  at  all  when  the  surface 
waves  and  radiation  lo.ss  are  strong,  due  to  their  interaction  with  the  microstrip 
fundamental  mode.  For  a  circuit  to  be  useful,  radiation  and  surface  wave  lo.ss 
should  be  as  small  as  possible.  With  a  careful  convergence  study,  it  is  found 
that,  within  the  useful  frequency  range,  typically,  19  piecewise  .sinu.soidal 
modes  of  size  0.06  guided  wavelength  (0.6  guided  wavelength  in  total)  can 
provide  results  within  a  few  percent  accuracy.  The  validity  of  the  current 
analysis  is  further  checked  with  the  quasi-static  method  at  low  frequencies. 
Fig.  3-9  shows  the  comparison  between  this  analysis  and  the  quasi-static 
method  for  a  single  layer  ca.se.  The  substrate  thickne.ss  is  1%  of  a  free  space 
wavelength  which  is  thin  enough  to  insure  the  accuracy  of  the  quasi-static 
method.  The  comparison  yields  very  good  agreement. 
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Energy  loss  due  to  radiation  and  surfaee  waves  at  a  microstrip  open  end  is 
shown  in  Fig.  3-10  with  and  without  a  cover  layer.  The  microstrip  line  in  this 
case  is  embedded  between  the  substrate  and  the  superstrate.  It  is  found  that 
with  the  presence  of  the  cover  layer,  the  loss  increases  with  the  increase  of  the 
superstrate  dielectric  constant  due  to  stronger  radiation  and  surface  waves. 
The  length  e.xtension  (or  capacitance)  at  the  open  end  due  to  the  fringing  field 
for  the  same  geometry  as  that  in  Fig.  3-10  is  shown  in  Fig.  3-11.  It  is  found, 
by  adding  a  cover  layer,  the  e.xcess  length  to  substrate  thickness  ratio  (or  end 
capacitance)  is  larger  due  to  a  stronger  fringing  field.  In  general,  the  e.xcess 
length  increa.ses  with  the  increa.se  of  effective  dielectric  constant  and  is  insen¬ 
sitive  to  frequency  c.xccpt  when  the  surface  waves  arc  strong.  The  excc.ss 
length  values  for  microstrip  in  a  composite  substrate  arc  shown  in  Fig.  3-12. 
Two  different  material  arrangements  were  investigated  which  include  the  case 
of  a  large  pcrmiiliN’ity  on  the  top  with  a  lower  one  on  the  bottom  and  the  case 
of  the  other  way  around.  It  is  found  that  the  exce.ss  length  value  for  the  first 
ca.se  is  significantly  larger  than  the  second  one.  This  implies  that,  when  the 
microstrip  line  is  on  the  larger  dielectric  constant  material,  the  fringing  electric 
field  is  stronger. 

For  the  gap  case,  after  a  matrix  inversion  in  Eq.  3-8  is  performed,  the  re¬ 
flection  coefficients  F  and  transmission  coefficient  T  arc  5,,  and  5,2  respec¬ 
tively.  Therefore,  the  admittance  matrix  of  the  gap  discontinuities  can  be 
obtained  by  the  following  transformation 

[y]  =  ([to  -  [S])([t0  +  [S])"'  (3.31) 
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where  [L’]  is  the  unitary  matrix.  By  comparing  the  two-port  n  network  and 
Fig.  3-3.  one  has 


^0 


=  -  >'l2 


and 


(3.32) 


+  jo)C^ 


(3.33) 


The  results  for  the  gap  discontinuity  arc  first  compared  with  those  obtained 
by  the  quasi-static  method  [18]  and  those  obtained  by  measurements  [68] 
and  are  shown  in  Fig.  3-13.  Since  the  gap  capacitance  is  small  and  is  sensitive 
to  the  device  tolerances,  the  measurement  is  inherently  difficult  to  perform 
accurately.  In  this  analysis,  the  frequency  is  chosen  to  be  5  GHz.  It  is  found 
that  this  dynamic  model  agrees  well  with  the  quasi-static  approach.  Some 
discrepancies  for  large  gap  spacing  may  be  due  to  the  fact  that,  in  such  cases, 
the  amount  of  energy  coupled  through  the  gap  is  comparable  to  the  energy 
losses  due  to  surface  waves  and  radiation,  and  this  aspect  is  not  included  in  the 
quaskstatic  approach.  The  gap  conductances  Gp  and  G^  arc  shown  in  Fig. 
3-U  for  the  .same  material  arrangements  as  those  shown  in  Fig.  3-10.  It  is 
found  that  the  cover  layer  (superstrate)  will  increase  conductance  due  to 
stronger  fringing  field  and  more  energy  losses.  Also  since  Gp  +  jojCp  repres¬ 
ents  the  input  admittance  when  a  perfect  magnetic  wall  is  in  the  middle  of  the 
gap,  it  is  c.\pei.*ed  that,  due  to  image  cancellation,  G^  will  be  much  larger  than 
Gp  .  As  shown  in  Fig.  3-14,  G^  is  about  two  orders  larger  than  Gp.  For  a 
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narrower  gap.  C  will  be  larger  than  due  to  stronger  end  coupling.  How¬ 
ever  for  wide  gap  spacing,  the  input  admittance  seen  from  cither  side  of  the 
gap  is  mainly  the  open-end  admittance,  so  will  be  larger  than  Cg.  The  re¬ 
sults  for  the  microstrip  gap  on  top  of  the  superstrate  arc  shown  in  Fig.  3-15. 
It  is  found  that  Cg/w  is  more  frequency  dependent  than  Cp/w,  when  the 
dispersion  is  strong. 
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r,  =  9/i  and  h  -  0.01  A„. 
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Chapter  IV 

MIC  Discontinuities  11- 
Slot  Line  Short-end  Discontinuities 


4.1  Introduction 

Slot  line  was  first  proposed  by  Cohn  in  1968  [69]  as  a  waveguiding  struc¬ 
ture  for  MIC  applications.  The  basic  structure  consists  of  a  narrow  slot  in  a 
conducti\  c  coating  on  one  side  of  a  substrate,  the  other  side  of  .the  substrate 
being  interfaced  with  air.  The  geometry  of  a  slot  line  is  shown  in  Fig.  4-1. 
In  order  that  electromagnetic  waves  be  confined  in  the  vicinity  of  the  slot,  the 
permittivity  of  the  substrate  is  usually  high  ( >  10).  This  slot  line  structure  has 
been  used  in  filters,  couplers  and  circuits  containing  semiconductor  devices 
[70],  [I3].  Since  the  electric  field  of  the  guide  wave  is  approximately  trans¬ 
verse  to  the  wave  propagation  direction,  the  fundamental  mode  is  TE-like. 
There  is  no  cut-off  frequency  for  this  structure,  since  the  slot  separates  two 
semi-infinite  ground  planes.  Slot  lines  can  also  be  included  in  microstrip  cir¬ 
cuits  by  etching  the  slot  circuits  in  the  ground  plane.  With  the  slot  and 
microstrip  combinations,  many  circuits  have  been  realized;  for  example,  hybrid 
branchline  directional  couplers  [71]  and  microstrip  bandstop  filters  [16]. 
Other  discontinuities,  such  as  the  microstrip-slot  transition  for  a  two-level  cir¬ 
cuit  design,  and  re.sonant  slots  for  antenna  applications,  will  be  discussed  in 
Chaps.  V  and  VII  respectively.  The  characteristics  of  a  slot  line  including  its 
propagation  constant  and  characteristic  impedance  have  been  studied  exten¬ 
sively  [71]  -  [74];  however,  analytic  methods  for  slot  line  discontinuities  are 
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scarce,  especially  in  terms  of  including  the  radiation  and  surface  waves  effects. 
This  slot  line  short-end  discontinuity,  as  shown  in  Fig.  4-2,  is  the  most  com¬ 
monly  seen  discontinuity  in  slot  line  circuits  [74].  It  can  be  used  in  the  design 
of  matching  networks,  filters  and  couplers.  Since  the  current  will  flow  near  the 
end  of  the  slotline,  there  is  an  apprecialble  amount  of  energy  storage  beyond 
the  termination.  The  end  discontinuity  will  cause  radiation  and  generates 
surface  waxes.  Without  the  image  cancellation,  energy  storage  and  losses  are 
more  .severe  than  in  the  microstrip  line  case,  and  this  end  effect  has  to  be  ac¬ 
counted  for  in  accurate  circuit  designs. 

In  slot  line  structures,  the  larger  the  effective  dielectric  constant,  the  more 
the  energy  is  confined  to  the  vicinity  of  the  slot.  If  the  other  side  of  the  ground 
plane  is  coxcred  by  a  dielectric  material  (instead  of  being  free  space),  the  cf- 
fectix'e  dielectric  constant  of  the  slot  line  will  increase.  This  offers  the  advan¬ 
tages  of  shorter  wax'elcngth.  greater  confinement  of  electromagnetic  fields,  and 
stronger  coupling  between  slot  line  circuits  [16].  The  geometry  of  a  short-end 
slot  line  sandwich  is  shown  in  Fig.  4-3.  The  comparison  of  the  end  effects  in¬ 
cluding  fringing,  radiation  and  surface  waves  between  a  slot  line  and  a  slot  line 
sandwich  will  be  provided  in  this  chapter. 

The  analysis  of  slotline  dhscontinuiiics  is  almost  a  one-to-one  correspond¬ 
ence  with  microstrip  discontinuities.  The  numerical  methods  discussed  in  the 
prexlous  chapter  can  be  modified  and  applied  to  the  slot  line  structure. 
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4.2.  Integral  Equation 


Since  the  one  layer  slot  line  is  the  case  of  the  slot  line  sandwich  with 
22  =  1 ,  only  the  analysis  for  the  latter  case  will  be  discussed.  For  the  short- 
end  slot  line  under  consideration,  an  integral  equation  can  be  formulated  in 
terms  of  a  transverse  electric  field  in  the  slender  slot,  viz.. 


pO  rwjz 
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y,  z  =  0“ )  =  G , ^(.v,  y |.v„  y^)  £  (.v^,  y^)  ci y^d . 
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Hp.y,z=‘Qr)  -  //,.(.v,y,r  =  0~)  =  0 


on  the  slot  line,  where  is  the  magnetic  field  due  to  the  electric  field  on  the 
slot  line.  The  Green's  function  can  be  found  from  Eqs.  2.47  and  2.48  of 
Chap.  11,  while  due  to  .symmetry  is  equal  to  G,^  e.xcept  the  sign  and  mate¬ 
rial  parameters.  Integral  equation  in  Eq.  4.3  can  be  solved  by  the  method  of 
moments.  This  procedure  is  the  same  as  that  for  the  open-end  microstrip  line. 
4.3.  Results 

From  Schelkunoff's  equivalence  principle,  the  electric  field  in  the  slot  is 
equivalent  to  the  magnetic  source  in  the  ground  plane.  From  image  theory 
[61],  the  magnetic  .source  and  its  image  have  the  same  magnitude  and  phase. 
On  the  other  hand,  in  the  microstrip  structure,  the  electric  current  and  its  im- 


age  arc  180^’  difference  in  phase.  Therefore,  the  slot  line  circuits  tend  to  radiate 
more  energy  than  microstrip  circuits.  An  example  of  radiation  and  surface 
wa\  e  losses  at  a  slot  line  short-end  in  terms  of  total  incident  power  versus  fre¬ 
quency  for  three  different  slot  widths  is  shown  in  Fig.  4-4.  It  is  seen  that  the 
energy  losses  increase  with  the  increase  of  slot  line  width  or  frequency.  This 
is  because,  by  increasing  the  frequency  or  slot  width,  less  energy  is  confined 
near  the  slot.  The  above  discussion  also  explains  the  result  in  Fig.  4-5  where 
it  is  shown  that  the  normalized  reactance  due  to  fringing  increases  with  the 
increase  of  slot  width  or  frequency.  It  is  seen  from  Fig.  4-4  that,  for  the 
substrate  with  £,  =  12,  when  the  substrate  thickness  is  about  O.O6A0,  more 
than  25"o  of  the  incident  power  is  lost  at  the  short  end  for  all  the  three  slot 
widths.  Fig.  4-4  provides  an  idea  of  the  upper  frequency  where  a  slot  line 
circuit  is  still  useful.  It  is  also  observed  that  even  for  a  thin  substrate  thickness 
(0.01/.„),  the  radiation  and  surface  wave  losses  are  5  to  10  %  of  the  total  inci¬ 
dent  power.  This  is  one  of  the  main  disadvantages  of  this  one-sided  slotlinc 
structure. 

By  adding  a  cover  layer  on  the  frec-space  side  of  the  slot  line,  the  guide 
wavelength  will  decrease  and  more  energy  will  be  confined  near  the  slot  region. 
This  implies  that  the  characteristic  impedance  will  decrease.  The  result  is,  that 
with  a  unit  voltage  wave  incident  to  the  slot  line,  the  total  incident  power 
will  increase  quite  noticeably  (typically,  5  to  10  %).  On  the  other  hand, 
for  a  unit  voltage  wave,  the  radiation  and  surface  wave  lo.s.ses  at  the  short  end 
increase  only  slightly  when  a  layer  is  added.  Therefore,  it  is  expected  that  the 
amount  of  the  energy  loss  at  the  discontinuities  can  be  reduced  by  using  a  slot 
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line  sandwich.  An  c.xample  of  this  result  is  shown  in  Fig.  4-6.  It  is  seen  that 
the  radiation  and  surface  wave  los.ses  are  greatly  reduced.  It  is  also  seen  from 
Fig.  4-7  that,  by  using  the  slot  line  sandwich,  the  fringing  fields  and  normal¬ 
ized  reactance  increase  due  to  the  higher  effective  dielectric  constant  and  to 
more  energy  being  confined  near  the  slot  line.  This  implies  that  the  coupling 
between  slot  line  circuits  is  enhanced. 
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Chapter  V 


A  Dynamic  Model  for  a  Microstrip-Slotline  Transition 

5.1  Introduction 

Based  on  the  approach  described  in  the  Chap.  Ill,  a  dynamic  model  for  a 
fnicrostrip-slotlinc  transition  and  its  related  structures  (such  as  a  microstrip  fed 
slot  and  a  slotline  fed  printed  dipole)  is  proposed  in  this  chapter.  The  devel¬ 
oped  model,  with  some  modifications,  can  be  applied  to  other  types  of  transi¬ 
tions  in  MIC  and  MM  1C  design.  Some  of  these  examples  will  be  discussed  in 
the  next  chapter. 

In  a  microstrip-slotlinc  transition,  a  short-circuit  slotline  which  is  etched  on 
one  side  of  the  substrate  is  cro.sscd  at  a  right  angle  by  an  open-circuit  micro¬ 
strip  on  the  opposite  side.  This  type  of  trdnsition  makes  a  two-level  circuit 
design  possible  [16].  Some  experimental  work  has  been  reported  [13],  [75] 
and  a  transmission  line  circuit  model  has  been  described  in  [76].  In  the 
present  approach,  the  radiation  and  surface  waves  due  to  the  cro.ss-ju notion, 
the  line  di.scontinuities,  and  all  the  mutual  coupling  due  to  the  dominant  mode 
as  well  as  higher  order  modes  of  each  line  arc  included  in  the  method  of  mo¬ 
ments  solution.  The  VSWR  and  input  impedance  of  the  transition,  can  be 
determined  by  the  current  distribution  on  the  microstripline  in  conjunction 
with  transmission  line  theory.  In  the  formulation  procedure,  certain  impor¬ 
tant  problems  in  MIC,  MM  1C  or  printed  antennas  design  can  also  be  solved. 
This  aspect  will  be  discussed  in  Chapter  Vlll. 


5,2.  Theory 

(5.2.1)  Green's  Function  Formulation. 

The  microstrip  to  slotline  transition  is  shown  in  Fig.  5-1  and  the  cross  sec¬ 
tion  is  shown  in  Fig.  5-2,  where  the  lines  are  e.xtcndcd  a  certain  distance  be¬ 
yond  the  cross-junction,  so  that  their  e.xtension  may  act  as  a  tuning  stub.  Due 
to  the  assumption  that  the  strip  or  slot  is  slender,  the  transverse  vector  com: 
ponents  {Jy  and  on  the  lines  arc  a  second  order  effect,  and  are  neglected 
for  simplicity.  Therefore,  only  the  Jc-dircctcd  electric  surface  current  on  the 
strip  and  the  v-directed  magnetic  surface  current  M„y  (.v-directed  electric  field 
arc  considered.  Under  the  above  assumptions,  the  coupled  integral 
equations  can  be  formulated  in  terms  of  and  Hy.  As  a  result,  the  electric 
Ey  at  (x,  y,  d)  due  to  the  prc.scncc  of  both  strip  and  slot  is 

E^  =  J  +  J  j*  6'jy,  ,  (5.1) 

and  the  difference  of  the  magnetic  field  Hy  at  (x,  y,  0+)  and  at  (x,  y,  0")  is 

A//^  =  j*  J  (7^,^  y,.  +  J  J  </.Vj ,  (5.2) 

where  and  Gy^  arc  the  dyadic  Green's  function  components  due  to  an  Jc  - 
directed  infinitesimal  electric  dipole  at  z  =  d  and  G^y  and  Gyy  arc  the  dyadic 
Green's  function  components  due  to  a  >•  -  directed  infinitesimal  magnetic 
dipole  at  z  =  0.  is  the  current  on  the  microstrip  s„  while  M^y  is  the  magnetic 
current  (electric  field)  on  the  slotlinc  . 
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Figure  5-1.  Top  view  of  microstrip-slot  line  transition. 
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Figure  5-2.  Cross  section  of  microstrip-slot  line  transition. 
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The  dyadic  Green's  function  components  G^y,  Gy^  and  Gyy  can  be  derived 
from  Eqs.  2.13-2.14  and  2.47-2.48.  The  results  can  be  c.xpresscd  as 

3«)  =  I  I  G^^a^,  (5.3) 

— oo  oo 

C.y<-v.  vl  v  JV)  =  r  r  (5.4) 

OO  *^—00 

=  -  Gy,^{x,  vljc^^yg  (5.5) 

and 


3K,  .\g  =  j  j  GyJ^Xy,,  Xy)eJ'-^^  .I'o)  (5,6) 


— OO  OO 


where 


4n“A'„£^ 


,  A^<7i(1 -e^)sinh<7,r/ 
+  — 


DJ^X) 


DmDJX) 


sinh  <7if/,(5.7) 


C.v,«.v  V  = 

4rt“ 


<7l  .  Kx^^r  -  l)sinh(<7i^/) 

T 


D,{X)  DJiX)DJX) 


(5.8) 


GyyiX,,Xy) 


-J 


An'-Z.k,  L 


2  _  ^2.  <71  cosh  <71^7  4-  £^.sinh<7|<7 

'  ^  D,IX) 

*  DM)D„IX)  ? 


(5,9) 


and  with 

Z)^A)  =  q  sinh(<7j</)  -1-  r/,  cosh(<7jrO,  (5.10) 
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sinh((7,^/)  +  cosh(<7ji/). 


(5.11) 


Other  pertinent  parameters  have  been  defined  in  Chap.  II. 

(5.2.2)  The  Choice  of  Expansion  Modes. 

In  the  method  of  moments  procedure,  and  are  expanded  in  terms 
of  a  set  of  known  functions.  For  the  transition  under  consideration,  the  mod¬ 
eling  of  two  half  infinite  lines  is  necessary,  in  which  .several  mechanisms  are 
possible.  If  sub.sectional  expansion  modes  are  used  in  the  two  finite  lines  with 
a  d  gap  source,  in  order  to  characterize  the  cross-  junction  of  this  resonator, 
two  sets  of  wa\’c  amplitudes  on  each  line  corresponding  to  two  different  length 
of  the  parasitic  line,  are  required.  This  scattering  matrix  formulation  has  been 
found  to  be  very  sensitive  to  error.  A  more  reliable  method  is  to  simulate  the 
physical  situation  where  both  the  micro.striplinc  and  the  slotlinc  arc  terminated 
by  a  matched  load.  In  this  .scheme,  subsection  expansion  modes  are  used  in 
both  lines  near  the  cross-junction  region,  while  entire  domain  travelling  waves 
arc  used  to  rcpre.sent  the  transmitted  wave  in  the  parasitic  line  (slotlinc)  and 
the  incident  wave  and  the  reflected  wave  on  the  feed  line  (microstripllnc). 
Other  choices  of  expansion  modes  arc  also  possible,  such  as  using  subscctional 
basis  functions  in  the  feed  line  and  subscctional  modes  and  traveling  wave 
modes  in  the  parasitic  line,  or  the  traveling  wave  modes  starting  away  from  the 
cross-junction  such  that  the  mutual  coupling  of  traveling  wave  and  PWS 
modes  on  different  lines  is  negligible.  These  types  of  e.xpansion  modes  have 
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some  advantages  in  numerical  analysis  and  will  be  discussed  in  the  next  sec¬ 
tion. 

The  traveling  wave  modes  used  involve  microstriplinc  and  slotlinc  propa¬ 
gation  constants  A'„,  and  From  the  knowledge  of  and  k^,  the  unknown 
current  distribution  on  the  microstripline  can  be  expanded  as 

;V 

;(.V)  =  ^  i„Mx)  (5.12) 

while  the  unknown  electric  field  distribution  in  the  slotlinc  can  be  expanded 
as 


M 

ki}')  =  I ''  +  (5.13) 

where 

^ ^  ,5,14, 

-re-'*"-’'  (5,15) 

and 

j-t  =  7’^.yV  (5.16) 


Piecewise  sinusoidal  (PWS)  modes  are  used  as  subscctional  e.xpansion  modes 
and  are  defined  starting  from  the  end  of  each  line.  These  modes  were  shown 
in  Eqs.  3.6  and  3.7. 
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(5.2.3)  The  Method  of  Moments  and  Matrix  Formulation 

The  coupled  integral  equations  can  be  obtained  from  Eqs.  5.1  and  5.2  by 
forcing  the  boundary  conditions  that  the  E-  field  must  be  zero  on  the 
microstriplinc  and  the  H-  field  must  be  continuous  across  the  slotlinc.  When 
when  the  expansion  modes  arc  substituted  into  Eqs.  5.1  and  5.2  followed  by 
Galcrkin's  procedure,  one  has 

^cmact^  ’^icmact  J  [Ije//] 

where  the  submatriccs  arc  the  reaction  between  different  expansion  modes, 

[/]  is  an  /V  X  I  column  vector  with  elements  /,,  /j,  .  ,  Is  ,  and  [£]  is  an 

.-\/x  1  column  vector  with  elements  . .  £u  . 

(5.2.4)  Some  Aspects  of  the  Numerical  Analysis. 

The  formulation  in  the  la.st  section  is  quite  flexible  and  can  be  easily  modi¬ 
fied  to  other  types  of  mode  expansion  mechanisms.  For  example,  the  traveling 
wave  modes  may  start  more  than  a  wavelength  away  from  the  cross-junction, 
which  modifies  the  Fourier  transform  of  the  traveling  wave  mode  in  the  above 
formulation.  This  type  of  expansion  has  the  advantage  that  the  mutual  cou¬ 
pling  between  traveling  wave  modes  and  PWS  modes  of  the  other  line  ( 
{.Tfcmact^  and  is  negligible.  Therefore,  the  computation  effort  in  Eqs. 

5.17  can  be  reduced.  Besides,  the  solutions  arc  automatically  convergent  in 
the  sense  of  the  number  of  c.xpansion  modes.  Another  type  of  expansion  is  also 
possible  where  only  PWS  modes  arc  used  in  the  feed  line.  This  has  the  ad- 
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(5.17) 
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vantage  of  providing  some  insight  into  the  current  distribution  on  the  feed  line 
and  of  avoiding  the  computations  of  the  submatriccs  [  and  [  in 

Eq.  5.17.  The  above  different  mode  expansion  mechanisms  may  also  be  used 
to  check  the  convergence  and  stability  of  the  solution. 

5.3.  Numerical  results 

The  results  for  a  microstrip-slotlinc  transition  have  been  obtained  based  on 
the  developed  algorithm.  The  numerical  analysis  was  performed  on  the  IBM 
3090  system.  Typically,  for  each  data  set  it  takes  about  one  minute  and  thirty 
.seconds  of  computer  time  in  contrast  to  a  half  second  to  obtain  the  propa¬ 
gation  constant  although  a  lot  of  effort  has  been  made  to  reduce  the  com¬ 
puter  cost.  An  example  of  a  50  H  microstripline  to  a  80  Q  slotline  transition 
is  given.  The  results  of  the  VSWR  and  input  impedance  arc  shown  in  Figs. 
5-3  and  5-4  respectively.  The  results  for  the  VSWR  arc  first  checked  by 
interchanging  the  feed  line  and  parasitic  line.  The  differences  in  in  are  within 
2'’/o,  which  is  consistent  with  the  property  of  low  loss  two  port  networks.  The 
complex  rcfiection  coefficient  that  is  obtained  is  checked  further  by  changing 
the  number  of  modes  and  different  mode  expansion  mechanisms  as  described 
in  the  last  .section.  Two  sets  of  input  impedances  with  different  numbers  of 
expansion  modes  and  base  function  size  arc  shown  in  Fig.  5-4  to  illustrate  a 
convergence  test  example.  With  the  particular  device  parameters  chosen  it  is 
found  that  both  the  magnitude  and  phase  of  the  reflection  coefficient  converge 
very  well  (2  %  in  IFl  and  5°  in  phase  )  for  (7<  0.036/o.  However,  for  higher 
frequencies,  the  results  arc  more  unstable,  and  typically  the  results  are  5-10 
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%  accurate  in  |r|  and  10-15  degrees  in  phase  before  higher  order  modes  turn 
on.  This  behavior  may  be  due  to  the  reasons  that  1),  when  the  radiated  and 
surface  waves  are  not  weakly  excited,  the  transmission  line  theory  applied  to 
the  microstripline  or  slotline  is  only  an  approximation  and  the  mode  expansion 
approach  is  somewhat  of  a  brute  force.  2),  the  transverse  vector  components 
{Jy  and  which  are  neglected  in  the  present  investigation  become  more  im¬ 
portant  as  frequency  increases.  The  VSWR  obtained  by  the  transmission  line 
circuit  model  [76]  and  the  measurement  [75]  arc  also  shown  in  Fig.  5-3  to 
provide  a  comparison.  In  the  transmission  line  circuit  model  the  stub  length 
is  assumed  to  be  measured  from  the  center  of  each  line  and  the  propagation 
constants  and  and  the  excess  length  are  obtained  from  the  current  anal¬ 
ysis,  It  is  seen  from  Figs.  5-3  and  5-4  that  the  present  method  agrees  very  well 
with  the  circuit  model  in  the  low  frequency  range.  The  discrepancy  for  higher 
frequencies  is  probably  due  to  the  higher  order  modes,  surface  waves,  and  ra¬ 
diation  effects  which  arc  neglected  in  the  circuit  model.  The  measurements 
reported  in  [75]  show  a  wider  bandwidth  than  that  of  cither  the  circuit  model 
or  the  present  analysis.  It  is  believed  that  the  accuracy  of  the  device  param¬ 
eters,  the  non-ideal  matched  load  and  (c.spccially),  the  coaxial  to  microstripline 
transition  affected  the  frequency-dependent  results  in  the  measurement.  Be¬ 
sides,  the  material  u.scd  in  [75]  is  Custom  HiK  707-20  {c,  =  20)  which  is 
usually  very  lossy  especially  for  higher  frcq.uencics.  These  effects  may  explain 
the  discrepancy  between  theory  and  experiment. 
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Figure  5-3.  VSWR  versus  frequency  for  microstrip-slot  line  transition. 
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AAA  transmission  line  circuit  model 
o  o  o  present  theory,  M=N=4  and  h=0.03 

■*  A  A  present  theory,  M=N=«7  and  h=0.025  Xq 


Chapter  VI 


Proximity  Coupled  Microstrip  Transition 
in  Double  Layer  Integrated  Circuits. 


6.1.  Introduction 

Analytic  and  numerical  methods  described  in  the  previous  chapter  are  e.x- 
tended  here  to  study  two  types  of  proximity  coupled  open  end  microstriplines 
in  a  double  layer  planar  structure.  These  pro.ximity  coupled  transitions  con¬ 
stitute  potentially  important  components  for  MIC  and  MMIC  design.  Fig.  6-1 
shows  two  semi-infinite  coHinear  microstriplines  at  different  levels.  This  type 
Of  transition  has  the  advantage  over  the  end  coupled  lines  in  that  the  overlap 
distance  4/  may  be  used  to  control  the  coupling.  Also  this  transition  provides 
a  wider  range  of  coupling  coefficient  with  a  reasonably  large  bandwidth  and 
it  therefore  is  useful  in  coupler  or  filter  design.  Fig.  6-2  shows  two  EMC 
transxerse  microstriplines.  In  this  type  of  transition,  an  open-circuit 
microstripline  printed  on  top  of  the  superstrate  is  crossed  at  a  right  angle  by 
another  open-circuit  microstriplinc  embedded  on  the  substrate.  These  two 
lines  are  extended  a  certain  distance  beyond  the  cross-junction  to  provide 
tuning  stubs.  This  type  of  transition  has  the  properties  of  broadband  and  good 
match  due  to  the  presence  of  the  double  stub.  The  materials  in  the  substrate- 
superstrate  configuration  may  greatly  affect  the  coupling  in  the  transition  and 
this  issue  has  also  been  investigated.  In  Section  6.2,  the  method  of  moments 
solution  of  integral  equations  is  formulated.  In  Section  Ill,  the  results 
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microstrip  2 


Figure  6-1.  Proximity  coupled  collinear  microstrip-microslrip  transition. 


Figure  6-2.  Proximity  coupled  transverse  microstrip-microstrip  transition. 
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from  the  numerical  analysis  arc  presented  and  some  interesting  properties  of 
the  above-mentioned  transitions  arc  discussed. 

6.2.  Analysis 


(6.2.1).  EMC  collincar  microstriplines. 

Integral  equations  for  EMC  collincar  microstriplincs  are 

2 

y=i 

where  £<'>  is  the  electric  field  at  microstrip  i,  i  =  1  or  2.  Here  microstrip  1  is 
at  z  =  b  while  microstrip  2  is  at  z  =  h.  The  Green's  function  (7,y  is  at  (x, 
y)  of  microstrip  i  due  to  an  .v  directed  delta  source  at  of  microstrip  j. 
This  Green's  function  has  been  derived  in  Chap.  11  and  can  be  cxprcs.scd  as 

Gij  =  I"’  r  (^-2) 

— OO  — OO 


where 


Dp.^..Ay) 


~J^()  r^2 


(6.3) 


The  function  /y  and  gy  can  be  identified  from  the  Green's  function  described 
in  the  Chap.  11.  Other  parameters  arc  defined  in  Eqs.  2.20-2.31.  When  the 
expansion  modes  (similar  to  those  in  Chaps.  Ill  and  V)  arc  used  in  Eq.  6.1, 
followed  by  Galcrkin's  procedure  in  the  same  way  as  described  in  the  last 
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chapter,  integral  equations  are  converted  into  a  set  of  linear  equations.  These 
M  +  N  +  2  equations  when  expressed  in  matrix  form  are 

The  computation  of  each  matrix  clement  in  Eq.  6.4  requires  a  double  infi- 

\ 

nitc  integration  where  the  integrand  contains  the  corresponding  Green's  func¬ 
tion  /.y)  and  the  Fourier  transform  of  the  current  expansion  functions. 
For  example  is  an  (/V  +  1)  x  M  matrix  with  matrix  elements 

i-ooX-oo  ^  (65) 

cosC-?.^n</j  +  mdi  - 

where  is  the  microstrip  overlap  length,  and  all  other  parameters  in  Eq.  6.5 
has  been  described  in  Chap.  111. 

(6.2.2).  EMC  transverse  microstriplincs. 

The  analysis  of  EMC  transverse  microstriplines  is  almost  one  to  one  in 
correspondence  to  the  collincar  case.  In  matrix  formulation,  if  local  coordi¬ 
nates  are  u.scd,  it  can  be  easily  shown  that  the  self  reaction  in  each  microstrip 
is  identical  for  the  collincar  and  transverse  ca.scs.  To  be  more  specific,  the  in¬ 
tegral  equations  for  the  transverse  microstriplincs  arc 

4'*  =11  +  J  j  Vf  (6-6) 


L/,J 

-r 

i:/2] 

T  . 


(6.4) 
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and 


4'*  =  J  J  +  j  j  (6.7) 

where  £y>  and  B/*  are  the  longitudinal  electric  fields  at  microstrip  1  (at  z  = 
b)  and  microstrip  2  (at  z  =  h),  respectively.  The  function  equals  to  G,,, 
while  Gyy  is  the  same  as  except  and  Xy  arc  interchanged  in  Z)22-  Other 
Green's  functions  arc  ^ 

G^=r  r  (6.8) 

oo  oo 

and 

^yx  =  (6-^) 

where 


r 

2  ^ 
47r  A'q£2 


W” 


^x^y^2 

D^mX) 


(6.10) 


The  mode  expansion  mechanism  and  the  method  of  moments  procedure  follow 
in  the  same  manner  as  for  the  collinear  case.  The  final  matrix  is  in  exactly  the 
same  form  as  Eq.  6.4.  The  submatriccs  ,  [Z„e/y;.]  with  i  =*  1  or  2  and 

arc  identical  for  the  collinear  and  transverse  cases.  All  other  subma¬ 
triccs  can  be  obtained  in  a  similar  way  as  for  the  longitudinal  coupling  case. 
For  example, 
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Too  roo 

ZS,,  =  -l 

•'-oo’'-oo  (w-ll) 

sin[/j^rt«/,  -  4,)]  sinC/jXmrf,  “  » 

where  is  an  (A/  +  1)  x  N  matrix,  The  parameters  /„,  and  4,2  are  the 

stub  length  for  microstrips  1  and  2  respectively.  Sine  function  in  £<].  6.11  in¬ 
stead  of  a  cosine  in  Eq.  6.5  is  because  the  Green's  function  in  Eq.  6.11  is  an 
odd  function  of  cither  or  Ay . 

6.3.  Results 

Although  the  impedance  matrix  in  Eq.  6.4  looks  formidable,  the  computa¬ 
tion  can  be  simplified  further  based  on  some  physical  insight.  For  c.xample, 
the  Green's  function  and  basis  functions  are  the  same  in  each  submatrix  except 
for  a  translation  in  reaction  center.  Therefore  in  the  numerical  process,  these 
common  factors  need  to  be  computed  only  once.  Also  due  to  reciprocity,  only 
a  fraction  of  the  impedance  elements  need  to  be  computed.  In  the  computa¬ 
tions  for  the  transition  problem,  entire  domain  modes  of  three  and  a  half 
guided  wavelength  long  and  eight  to  thirteen  PWS  modes  (depending  on  the 
overlap  or  stub  length)  arc  used  in  each  microsiriplinc.  The  convergence  has 
been  checked  for  5,,  (F),  to  within  3%  in  magnitude  and  3®  in  phase.  The 
magnitude  of  the  reflection  coefficients  is  shown  in  Fig.  6-3  as  a  function  of 
overlap  for  the  collincar  transition  with  three  types  of  material  arrangements. 
The  corresponding  microstrip  width  is  chosen  such  that  the  microstrip  lines 
have  a  50  Q  characteristic  impedance.  For  the  case  of  a  large  dielectric  con¬ 
stant  material  in  the  substrate  and  a  smaller  one  in  the  superstrate,  since  en- 
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ergy  is  moslly  confined  In  the  substrate,  less  power  is  transmitted  than  that  of 
"Other  types  of  material  arrangements.  This  behavior  is  observed  in  Fig.  6-3. 
From  these  two  figures,  one  can  sec  that  the  coupling  coefficient  depends  on 
the  ^amount  of  energy  of  an  embedded  microstrip  stored  in  the  superstrate. 
The  relationship  between  overlap  length  and  power  distribution  obser\ed  in 
Fig.  6-3  is  not  obvious.  It  can  be  explained  empirically  as  follows.  The 
amount  of  current  induced  in  the  parasitic  line  is  mainly  due  to  the  longitudi¬ 
nal  electric  field  generated  by  the  open  end  feed  line.  This  current  varies 
sinusoidally  and  decreases  as  the  observation  point  moves  further  away  from 
the  open  end.  Therefore,  as  the  two  microstrips  arc  brought  closer,  increased 
coupling  occurs.  As  the  coupling  gets  stronger,  the  electric  field  due  to  the 
parasitic  line  will  interact  with  the  feed  line  field.  Since  these  two 
microstriplincs  have  different  fundamental  modes,  as  the  overlap  increases 
further,  the  coupling  starts  to  decrease  due  to  wave  dcstructi\c  interference. 
It  is  interesting  to  sec  that  with  a  particular  overlap  length,  very  little  coupling 
occurs.  This  behavior  is  found  to  be  related  to  the  superstrate  thickness  and 
dielectric  constant.  It  is  also  ob.scrvcd  that  as  the  overlap  gets  larger,  the  re¬ 
flection  becomes  smaller.  This  implies  that  in  such  a  ca.se,  the  guided  funda¬ 
mental  mode  is  more  like  the  coupled  line  mode.  It  is  further  found  from  Fig. 
6-3,  that  in  a  certain  region  where  coupling  reaches  a  local  ma.ximum,  the 
scattering  matrix  is  insensitive  to  overlap  length.  Since  the  line  impedance  and 
effective  dielectric  con.stant  arc  also  frequency  insensitive,  in  this  microstrip 
transition,  only  the  effective  overlap  length  will  be  frequency  sensitive.  This 
implies  that  this  transition  is  broadband.  An  e.xamplc  is  shown  in  Fig.  6-4. 
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One  can  see  that  the  scattering  parameters  change  no  more  than  3%  in  mag¬ 
nitude  fur  this  particular  X  band  computation. 

For  the  transverse  microstrip  transition,  the  material  also  has  a  strong  effect 
on  the  coupling  mechanism,  as  shown  in  Figs.  6-5  and  6-6.  The  coupling  be¬ 
tween  two  transN  crse  microstrips  is  further  complicated  by  the  presence  of  the 
two  tuning  stubs.  The  effect  of  these  two  stubs  is  very  much  different  from 
those  in  the  microstrip-slotlinc  transition  where  optimum  coupling  occurs  when 
both  stubs  arc  about  a  quarter  wavelength  long.  For  the  transverse  microstrip 
transition,  the  coupling  is  minimum  when  either  stub  is  about  a  quarter 
wavelength  and  is  maximum  when  both  stubs  arc  a  half  wavelength  long. 
This  phenomenon  is  due  to  the  fact  that  the  parasitic  line,  from  a  circuit  point 
of  view,  is  a  shunt  element  to  the  feed  line,  and  vice  versa.  Therefore,  when 
both  stubs  arc  a  half  wavelength  long,  the  circuit  (looking  from  the  cross 
junction)  is  in  resonance,  while  when  either  stub  is  about  a  quarter  wavelength 
long,  the  circuit  is  in  effect  shorted.  For  the  parameters  in  Fig.  6-5,  the  guide 
wavelength  for  each  microstrip  is  approximated  as  =  0.547^o  ^nd 
/.,„2  =  0,507/.o  ,  while  in  Fig.  6-6,  the  guide  wavelength  is  =  0.32Uo  for 
microstrip  I  and  X„2  ~  0.360Ao  for  microslrip  2.  In  Fig.  6-5,  maximum  cou¬ 
pling  occurs  when  both  stubs  and  arc  about  a  half  guide  wavelength 
long.  In  this  particular  ease,  one  can  see  that  the  VSWR  of  this  transition  can 
be  as  small  as  1.1.  Therefore,  this  transition  can  be  potentially  useful  in  a  two 
level  circuit  design. 

To  verify  the  analysis,  a  3  inch  by  3  inch  circuit  was  built  and  tested  for  the 
ease  of  a  transverse  microstrip  transition.  Duroid  materials  with  permittivity 
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2.2  and  10.2  were  the  substrate  and  the  superstrate  respectively.  The  dimen¬ 
sions  of  the  dc\  ice  were  chosen  to  be  the  same  as  those  in  Fig.  6-5  except  that 
stub  lengths  of  about  a  half  guide  wavelength  were  used  (0.75  cm  and  0.81  cm 
for  the  top  and  bottom  microstrip  rcspccthcly).  The  circuit  was  made  using 
a  standard  photo-etching  technique  and  was  measured  on  an  HP-8510  net¬ 
work  analyzer.  Both  computed  and  measured  results  for  VSWR  arc  shown  in 
Fig.  6-7.  The  comparison  shows  good  agreement.  The  ripple  observed  in  the 
measurement  may  be  due  to  an  imperfect  match  at  the  coaxial-microstrip 
transitions.  One  can  observe  that  the  VSWR  is  less  than  1 .8  from  7  to  1 1 
GHz.  Such  a  broadband  transition  would  be  very  useful  in  circuit  design. 
This  broadband  property  is  mainly  attributed  to  the  double  resonance  due  to 
the  presence  of  double  stubs.  In  this  inve.stigation,  50  Q  microstrip  lines  were 
used.  The  impedance  levcl  will  affect  coupling  in  the  transition.  Therefore,  the 
results  presented  here  may  not  be  optimal.  The  choice  of  the  impedance  level 
may  depend  upon  the  purpose  of  the  circuit  design. 
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Figure  6-4,  IH  and  \7\  versus  frequency.  =  0.15  cm,  Ei  »  2.2  and  £2  =  10.2, 
h  «  50  mil,  b  ■  25  mil,  h-,  s  42  mil  and  w,  *  76  mil. 


Figure  6-5.  IFj  as  a  function  of  stub  length. 

30n-30n  transverse  microstrip  transition, 
f  «  10  GHz,  h  ■  30  mil,  b  «  23  mil,  =  42  mil  and  =  76  mil. 
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Figure  6-6.  |71  as  a  function  of  stub  length  for  a  500-300  transverse  transition, 
f «  10  GHz,  h  «  50  mil,  b  >  25  mil,  m'|  s  16  mil  and  ^2  -  50  mil. 
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Chapter  VII 

Printed  Slot  Characteristics 


7.1.  Introduction 

The  printed  slot  antenna  is  a  radiating  element  within  the  category  of  inte¬ 
grated  circuit  antennas.  The  basic  structure  consists  of  a  grounded  substrate 
with  a  slot  etched  in  the  ground  plane.  The  geometry  is  shown  in  Fig.  7-1. 
This  printed  slot  structure  has  the  features  of  producing  bidirectional  radiation 
patterns  and  offering  an  additional  degree  of  freedom  when  combined  with 
microstrip  dipoles  or  patches  [79]  .  The  combination  of  printed  slots  and 
dipoles  also  enables  the  realization  of  circularly  polarized  radiation  patterns 
[55].  In  this  chapter,  some  characteristics  of  printed  slot  antennas,  such  as 
efficiency,  power  distribution  in  each  region,  and  the  effects  of  materials  are 
studied.  In  Section  II,  Green's  function  of  an  infinitesimally  small  rnagnctic 
dipole  in  the  ground  plane  is  used  to  study  the  material  efTects  on  the  radiated 
power  in  cither  side  of  the  ground  plane  and  the  surface  wave  power  in  the 
substrate.  In  Section  III,  substrate  effects  on  the  resonant  length  and  the  input 
admittance  of  a  center-fed  printed  slot  are  investigated  through  a  method  of 
moments  .solution  of  a  magnetic-type  integral  equation, 

7.2.  Substrate  ElTects  on  Printed  Slot  Properties 

The  slot  antenna  geometry  under  consideration  is  shown  in  Fig.  7-2.  The 
problem  consists  in  its  simplest  forms  of  an  Jr-  directed  infinitesimal  magnetic 
dipole  skintight  against  the  ground  plane.  The  electromagnetic  fields  of  this 


77 


A 

z 


Figure  7-2.  A  Hertzian  magnetic  dipole  in  the  ground  plane. 
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geometry  follow  from  the  derivation  of  the  Green  function  presented  in  Chap. 
II  and  can  be  described  in  term  of  a  Sommerfcld-typc  [62]  integral. 

The  radiating  and  surface  waves  of  interest  here  arc  the  electromagnetic 
waves  in  the  far  zone,  where  radiating  fields  arc  spherical  waves  with 
R  =  z~  -» oo,  and  where  surface  waves  are  cylindrical  waves  with 

p  -*oo.  The  radiating  fields  can  be  obtained  by  the  method  of  steepest  descent 
[so];  while  the  surface  waves  can  be  found  by  applying  the  Residue  and 
Cauchy  Ricmann  theorems  [SO].  The  final  results  for  the  far  zone  radiating 
fields  can  be  summarized  as 


Eff  =  —j<opQ  sin  (f>  cos  dkQ 
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-JIcqR 


DJ^k^sind)  R 


(7.1) 


g-Jki)R 

cos  cos  OAtq  A/(A:q  sin  0) — - — 


(7.2) 


M{X  =  A'o  sin  (?)  =  cos  ^k^f 


2n}  2/1  -  n?)  sinh  <716 

=  CHS  dk,i - - ! - 


D^X)  D„j,X)DJ,X) 


(7.3) 


J±. 

mu 


and 


I20n  ' 


The  radiating  power  above  the  substrate  is 


(7.4) 
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(7.6) 
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The  0  integration  is  carried  out  anahlically  while  6  is  computed  numerically. 

It  can  be  shown  by  reciprocity  that  each  surface  wave  mode  propagates  in¬ 
dependently  [49].  The  surface  wave  modes  propagate  along  the  dielectric 
surface  and  decay  exponentially  toward  free  space.  For  the  TM  surface  wave 
modes  (transverse  to  z),  the  total  surface  wave  power  Pj-^  is  the  sum  of  the 
surface  wave  power  in  region  0,  {Ptmo)  and  in  region  I,  (P7■^|,)  where 
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and 

li,p.)  =  <7|  coshr/|6  4-  sinh<7|/>.  (7.9) 

For  the  TE  surface  waNc  modes,  the  total  surface  wave  power  Pj^  is  the  sum 
of  surface  wave  power  in  region  0,  (Fy^o).  and  in  region  1,  {Pte\)^  where 
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The  total  surface  \va\e  power,  is  the  sum  of  Pj^i  and  Pj^.  The  power  in 
the  free-space  side  of  the  ground  plane,  Pf,  can  be  obtained  from  the  expression 
in  Eq.  7.3  as  a  special  case 


An 

Pi,  =  ^(120;r)-. 


(7.12) 


For  printed  slot  antennas,  it  is  important  to  know  how  the  radiated  powers 
on  either  side  of  the  ground  plane,  and  the  power  losses  due  to  surface  waves 
arc  affected  by  the  substrate  thickness,  permittivity  and/or  permeability.  The 
radiation  efficiency  in  the  upper  half  of  the  ground  plane  (r  >  0)  is  defined  as 
the  ratio  of  radiated  power  to  total  power,  i.c 


= 


P  +  P 

*  <1  ~  '  .r.H'. 


(7.13) 


The  powers  Pg,  Fj  and  Fr.»,.as  a  function  of  the  substrate  thickness  with 
c,  =  4  and  /<,  =  I  are  shown  in  Fig.  7-3,  where  each  power  component  has 
been  normalized  to  Fj,  It  is  seen  that,  as  substrate  thickness  increases  from 
zero,  both  Pg  and  increase.  The  radiated  power  reaches  a  maximum  when 
the  first  TE  surface  wave  mode  turns  on.  When  the  substrate  thickness  in¬ 
creases  further,  both  Pg  and  decrease  until  the  next  surface  wave  mode 
nearly  turns  on.  The  local  maximum  of  radiated  power  is  due  to  the  fact  that 
when  a  surface  wave  mode  turns  on,  there  is  a  considerable  amount  of  power 
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radiating  along  the  horizon  [53].  The  oscillatory  behavior  of  the  surface  wave 
variation  with  substrate  thickness  can  be  explained  in  terms  of  ray  optics 
[57]  and  is  mainly  due  to  the  constructive  and  destructive  interference  of 
wavefronts  when  rays  bounce  back  and  forth  inside  the  substrate.  In  Fig.  7-4, 
the  power  distribution  versus  substrate  thickness  for  e,  =  12.5  and  /t,  =  1  is 
shown.  By  comparing  Fig.  7-3  and  Fig.  7-4,  one  can  see  that,  for  a  larger 
substrate  permittivity,  higher  order  surface  wave  modes  turn  on  at  smaller 
substrate  thickness.  It  is  also  observed  that  materials  with  larger  permittivity 
can  support  more  surface  wave  power  and  allow  more  radiation.  The  results 
show  that  e\en  though  both  radiated  and  surface  wave  power  in  the  region 
i  >  0  increase  with  the  increase  of  substrate  permittivity,  the  increase  of  the 
surface  wave  power  is  more  marked.  The  power  distribution  versus  substrate 
thickness  for  a  magnetic  material  e  =  I  and  /t,  =  4  is  shown  in  Fig.  7-5.  It 
is  obserxed  that  for  a  magnetic  substrate  P/,  is  greater  than  Pg.  Generally 
speaking,  in  printed  slot  structures,  more  power  radiates  in  the  substrate  side 
of  the  ground  plane  (P^)  than  that  in  the  free  space  side  (F*),  when  the 
substrate  is  a  dielectric;  while  the  reverse  is  true  for  a  magnetic  material.  It 
is  also  observed  that  the  surface  wave  power  for  magnetic  materials,  in  con¬ 
trast  to  dielectric  materials,  increases  even  after  the  first  TE  surface  wave 
mode  turns  on.  The  radiation  efficiency  versus  substrate  thickness  is  shown  in 
Fig.  7-6.  For  a  magnetic  material,  with  a  small  substrate  thickness,  surface 
waxe  poxver  increases  and  radiated  poxver  decreases  xvith  the  increase  of 
substrate  thickness.  Therefore,  it  is  seen  from  Fig.  7-6  that  efficiency  decreases 
drastically  versus  substrate  thickness  for  a  thin  substrate  of  magnetic  material. 
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For  /i,  =  4.£,  =  1,  less  than  20%  efficiency  is  observed  when  the  substrate 
thickness  is  only  0.1  Aq  . 


7.3.  A  Center-fed  Printed  Slot 

A  printed  slot,  center-fed  by  an  ideal  d-gap  source,  is  considered  in  this 
section.  The  geometry  is  shown  in  Fig.  7-7.  A  method  of  moments  solution 
of  integral  equation  is  used  here  to  study  the  characteristics  of  a  slender  rec¬ 
tangular  slot.  The  procedure  described  in  Chap.  V  for  the  slot  can  be  used  to 
obtain  the  following  linear  simultaneous  equations:  . 


cner:  =  [/]. 


(7.14) 


where  1',.  . .  I'y  are  the  amplitudes  of  the  e.xpansion  functions.  Sup¬ 

pose  that  an  odd  N  is  u.scd,  then  for  a  center-fed  slot, 
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and  /*  =  0  for  k  ^  (/V  +  1  )/2.  The  matri.x  elements  can  be  formulated  as 
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and  is  the  Fourier  transform  of  the  PVVS  expansion  mode.  The  method 
of  computing  Eq.  7.16  has  been  discussed  in  the  previous  chapters.  The  input 
admittance  is  defined  as 


^’in  =  j^s  = 


V, 


(s+m 


(7.19) 


Based  on  the  above  analysis,  the  input  admittance  of  a  center-fed  slot  is 
examined  as  a  function  of  various  structure  parameters.  In  Fig.  7-8,  the  input 
admittance  versus  slot  length  is  shown  with  two  different  material 
pcrmitti\itics.  It  is  observed  from  Fig.  7-8  that  when  the  permittivity  in¬ 
creases.  the  slot  resonant  length  decreases  and  resonant  resistance  increases. 
The  center-fed  slot  bandwidth  (BW)  can  be  defined  [46]  as 
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It  is  found  from  Fig.  7-8  that  with  h  =  0.02Ao  and  w  =  O.OUq  »  "’hen  e,  in¬ 
creases  from  2.2  to  5.0,  the  bandwidth  (BW)  decreases  from  21.5%  to  13.2%. 
The  input  admittance  versus  slot  length  with  two  different  substrate  thickness 
is  shown  in  Fig.  7-9.  It  is  observed  that  increase  of  the  substrate  thickness 
decreases  the  rc.sonant  length,  the  rc.sonant  resistance  and  the  bandwidth.  The 
effect  of  slot  width  on  the  input  admittance  is  shown  in  Fig.  7-10.  It  is  found 
that  the  increase  of  slot  width  causes  the  slot  resonant  length  and  resonant  re¬ 
sistance  to  decrease  and  the  bandwidth  to  increase. 
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Giapter  Vlll 

Printed  Antenna  Feeding  Structures 


8.1.  Introduction 

The  design  of  feeding  structures  is  one  of  the  most  important  parts  in  the 
design  of  printed  circuit  antennas.  Since  the  geometry  of  printed  circuit  an¬ 
tennas  involves  grounded  substrates,  the  design  of  feeding  networks  usually 
requires  integrated  circuit  technology.  For  monolithic  phased  array  applica¬ 
tions,  the  transition  from  feeding  networks  to  radiating  elements  is  the  key  to 
the  success  of  the  design.  Classical  feeds  of  printed  circuit  antennas  contain 
cither  a  microstrip  transmission  line  in  physical  contact  with  the  radiating  ele¬ 
ments  or  a  coaxial  line  penetrating  through  the  ground  plane  [81],  [28].  In 
the  former  case,  feed  and  antenna  together  form  a  resonator.  Due  to  the  high 
Q  of  this  type  of  resonator,  the  bandwidth  is  very  narrow.  Undesired  spurious 
radiation  at  the  microstrip-antenna  junction  may  also  be  a  serious  problem. 
The  probe  fed  microstrip  antenna  is  not  suitable  for  millimeter  wave  or 
monolithic  applications  mainly  due  to  the  presence  of  the  probe. 

To  overcome  the  disadvantages  of  classical  feeds,  a  new  fccd-antcnna  ar¬ 
rangement  using  electromagnetic  coupling  (EMC)  has  been  proposed 
[8l],  [56],  In  this  type  of  structure,  there  is  no  physical  connection  between 
feed  lines  and  radiating  elements.  The  microstrip  line  can  be  very  close  to  the 
ground  plane  to  reduce  spurious  radiation  and  the  antenna  is  away  for  the 
ground  plane  to  increase  the  bandwidth.  GaAs  substrate  (c,  =  12.8)  is  very 
suitable  for  active  devices,  but  its  dielectric  constant  is  too  high  for  radiating 
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elements  when  efficiency  and  bandwidth  are  considered.  EMC  antenna 
structures  are  excellent  for  thc.se  monolithic  phased  array  applications,  where 
active  devices  can  be  printed  on  a  GaAs  substrate  and  antennas  on  a  low 
permittivity  supcrstratc.  The  architectures  using  EMC  antennas  in  layered 
configurations  seem  to  be  the  ultimate  choice  for  monolithic  phased  arrays. 

The  EMC  collincar  dipole  has  been  c.xtensively  studied  in  the  past  few  years 
[59].  [60],  [42],  [43]  .  In  this  chapter,  four  EMC  antenna  feeding  struc¬ 
tures  arc  discussed.  These  include  a  microstrip  fed  slot,  a  slotline  fed  dipole, 
an  EMC  transverse  dipole  and  a  microstrip  fed  slot  coupled  dipole.  These 
structures  have  the  common  feature  that  a  string  of  the  antenna  clement  above 
a  common  feed  line  becomes  a  linear  array.  The  analyses  discussed  in  Chaps. 
V  and  VI  may  also  be  u.scd  to  .study  the  printed  circuit  antenna  structures. 
The  features  of  the  four  feeding  structures  together  with  numerical  results  will 
be  discussed  in  Sections  8.2-8.5. 

8.2.  A  Microstrip  Fed  Slot 

The  radiation  characteristics  of  a  printed  slot  antenna  were  discu.sscd  in  the 
last  chapter.  The  features  of  the  slot  antenna  as  compared  to  the  strip  dipoles 
are  that  it  is  broadband  and  gives  bidirectional  radiation.  A  simple  way  to 
design  its  feeding  structure  is  to  use  a  microstrip  line  electromagnetically  cou¬ 
pled  to  the  transverse  slot,  as  shown  in  Fig.  8-1.  The  analysis  of  this  structure 
has  been  reported  using  an  equivalent  circuit  model  [82]  ,  a  waveguide  model 
[83],  and  via  reciprocity  [84].  However  up  to  now  there  are  no  rigorous  re¬ 
sults  available.  In  this  section,  a  rigorous  result  based  on  solving  an  exact 
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integral  equation  by  the  method  of  moments  is  presented.  The  equivalent 
circuit  of  a  slender  rectangular  slot  perpendicular  to  the  microstrip  is  a  scries 
element.  In  order  to  understand  the  property  of  the  slot,  the  equivalent  circuit 
as  a  function  of  various  device  parameters  should  be  in\'cstigatcd.  The  nor¬ 
malized  equivalent  resistance  and  reactance  as  a  function  of  slot  length  are 
shown  in  Figs.  8-2  and  8-3  respectively,  for  two  different  substrate  dielectric 
constants.  It  is  observed  that  for  a  fixed  substrate  thickness,  the  resonant 
length  dccrca.ses  and  maximum  coupling  (resonant  resistance)  increases  with 
the  increase  of  the  substrate  permittivity.  These  phenomena  are  due  to  the 
fact  that  the  increase  of  substrate  permittivity  will  decrease  the  effective 
substrate  thickness  (effective  distance  between  microstrip  and  slot)  and  de¬ 
crease  the  effective  length  of  the  structure;  as  a  result,  the  resonant  length  dc- 
crca.ses  and  maximum  coupling  increases.  It  is  also  seen  that  for  smaller 
substrate  permittivity,  the  equivalent  circuit  is  less  sensitive  to  the  change  of 
the  structure  parameters. 

Fig.  8-4  shows  the  Smith  chart  plot  of  the  equivalent  impedance  of  the  slot 
as  a  function  of  offset  (<5)  for  three  dificrent  slot  lengths.  The  offset  6  is  de¬ 
fined  as  the  distance  between  the  center  of  the  slot  and  the  center  of  the 
microstrip  line.  It  is  ob.scrvcd  that  the  coupling  between  the  slot  and  the 
microstrip  line  decreases  monotonically  as  offset  increases.  It  is  interesting  to 
sec  that,  for  small  ofTset,  the  impedance  is  insensitive  to  the  offset  change. 
This  implies  that  the  alignment  between  the  slot  and  the  microstrip  line  is  less 
critical  than  for  EMC  dipoles  [60],  where  the  offset  control  is  very  important 
to  the  accuracy  of  a  design. 


Nomulizcd  equivalent  reactance  Normalized  equivalent  resistance 


rifiutc  8-2.  Norninli7.cd  tcMSinncc  of  n  inicto^dip  fed  slot. 
/i  =  0.()2X„  nml  *«•,  =  O.CMXn- 


rifiire  8  NoiinaliKd  icaclanre  of  a  50U  miciosliip  fed  slot. 
h  -  0.02).o  and  =  0.0 1  X„. 
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Figure  8-4.  Smith  chiit  plot  of  the  impedance  of 
a  microstrip  fed  slot  versus  offset, 
e,  =  2.2,  h  s  0.02X^,  w,  *  O.OIX^  and  *  3.125h. 
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The  input  impedance  of  a  stub-tuned  slot  obtained  from  information  above 
the  reflection  coefficient  in  the  microstrip  line  is  shown  in  Fig,  8-5  as  a  func¬ 
tion  of  frequency.  It  is  seen  that  the  bandwidth  is  mainly  determined  by  the 
tuning  stub  since  the  resistance  is  quite  insensitive  to  the  frequency.  There¬ 
fore,  to  increase  the  bandwidth,  the  stub  length  should  be  chosen  such  that  at 
the  resonant  frequency  the  change  of  stub  impedance  with  frequency  is  as 
small  as  possible.  Another  way  to  increase  bandwidth  is  to  control  the  device 
parameters  such  that  resonance  occurs  even  without  the  tuning  stub.  Fig.  8-5 
shows  a  typical  e.xample  for  this  design  where  the  bandwidth  is  6%  with  a  stub 
length  0.02A(,.  The  bandwidth  will  also  increase  if  the  slot  width  is  enlarged. 

8.3.  A  Slot  Line  Fed  Dipole 

A  new  feeding  structure  for  printed  dipole  antennas  is  shown  in  Fig.  8-6, 
where  a  slot  line  is  etched  in  the  ground  plane,  and  c.xcites  a  printed  dipole 
through  proximity  coupling.  In  order  that  the  slot  line  be  an  effective 
waveguiding  structure,  the  substrate  permittivity  .should  be  large  (>I0)  to 
confine  energy  along  the  slot.  The  equivalent  circuit  of  the  dipole  seen  by  the 
slot  line  is  a  shunt  element.  An  e.xample  of  this  result  as  a  function  of  dipole 
length  with  GaAs  substrate  (e,  =  12.8)  is  shown  in  Fig.  8-7.  It  is  observ'ed 
that  the  admittance  changes  more  quickly  with  antenna  length  than  for  printed 
slots.  The  equivalent  admittance  of  a  dipole  as  a  function  of  frequecy  is  shown 
in  Fig.  8-8  for  a  substrate  with  £t  =  12.  It  is  observed  that  the  admittance 
changes  drastically  for  a  small  change  of  frequency.  This  antenna  structure  is 
inherently  narrow  band.  Another  problem  is  that  the  efficiency  of  the  antenna 
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is  small  for  a  large  substrate  permittivity.  To  overcome  these  disadvantages, 
dipoles  can  be  printed  in  a  low  permittivity  substrate  and  at  the  same  time  a 
high  permittivity  substrate  can  be  used  on  the  other  side  of  the  ground  plane 
to  confine  energy  near  the  slot  line.  This  type  of  slot  line  sandwich  structure 
has  been  discussed  in  Chap.  IV,  Another  advantage  of  this  structure  is  that, 
for  monolithic  applications,  active  devices  can  be  made  on  a  GaAs  substrate 
and  isolated  from  the  antenna  elements  by  a  ground  plane.  This  type  of  ar¬ 
chitecture  can  reduce  spurious  radiation  due  to  active  devices.  An  example  of 
the  results  of  the  equivalent  admittance  of  a  dipole  in  a  two-side  slot  line 
structure  is  shown  in  Fig.  8-9.  The  results  are  for  the  case  when  the  dipole  is 
printed  on  the  substrate  with  e,  =  2.54,  while  the  other  substrate  is  Cj  ==  12.. 
From  the  results  in  Figs.  8-8  and  8-9,  it  is  seen  that  bandwidth  is  indeed  im¬ 
proved  by  using  a  slotlinc  sandwich.  It  should  also  be  mentioned  that  the  ra¬ 
diation  and  surface  wave  lo.sscs  due  to  the  fccdlinc  may  also  be  reduced  by 
using  a  slot  line  sandwich. 

8.4.  An  EMC  Transverse  Dipole 

A  dipole  that  is  clectromagnctically  coupled  to  microstrip,  but  oriented 
transverse  to  the  microstrip,  as  shown  in  Fig.  8-10,  is  one  of  the  printed  an¬ 
tenna  feeding  structures  proposed  by  Oltman  [56].  One  of  the  features  of  this 
EMC  transverse  dipole  is  the  light  coupling  between  the  dipole  and  the 
microstrip  line.  Therefore,  it  is  fairly  well  suited  for  radiating  element  in  large 
antenna  array.  A  theoretical  and  experimental  investigation  of  this  EMC 
transverse  dipole  has  been  reported  in  [77]. 
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transverse  dipole 


Figure  8-10.  The  electromagnetically  coupled  transverse  dipole. 
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Figure  8-11,  Comparison  between  theoretical  and  experimental  results 
for  the  input  impedance  of  an  EMC  transverse  dipole. 

C|  s  2.17,  C2  s  2.17,  At|  s  Aj  >  0.8  mm,  Wi  s  I  mm,  m-j  s  2.2  mm, 
Ar  =  0,  Av  =  4.S  mm  and  L  -  12.4  mm. 
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In  the  next  chapter,  array  designs  of  this  antenna  architecture  will  be  dis¬ 
cussed.  In  order  to  check  the  analysis  of  this  work,  the  result  of  the  present 
analysis  is  compared  against  the  experimental  results  presented  in  [77].  The 
comparison  is  shown  in  Fig.  8-11.  It  is  found  that  the  agreement  is  quite 
satisfactory. 

8.5.  A  Microstrip  Fed  Slot  Coupled  Dipole 

A  common  problem  for  the  antenna  architectures  discu.ssed  in  the  last  few 
sections  is  that  the  feed  line  and  radiating  elements  are  on  the  same  side  of  the 
ground  plane,  and  spurious  radiation  due  to  the  transmission  line  can  not  be 
eliminated  completely.  In  millimeter  wave  applications,  feed  line  radiation 
may  cause  severe  cross  polarization.  A  nice  way  to  alleviate  this  difficulty  is 
to  use  a  two-sided  structure,  as  shown  in  Fig.  8-12,  where  any  radiation  due 
to  the  feeding  network  and  active  devices  is  isolated  from  the  antenna  by  a 
ground  plane  [85],  [86],  The  dipole  is  actually  fed  by  a  secondary  source, 
namely  a  slot  in  the  ground  plane.  This  type  of  aperture-coupled  antenna 
architecture  seems  to  be  a  promising  candidate  in  millimeter  wave  monolithic 
phased  array  applications.  Since  the  combination  of  a  slot  and  a  dipole  pro- 
N’idcs  sufficient  degrees  of  freedom,  the  dipole  is  centered  to  the  slot  and  the 
slot  is  centered  to  the  microstrip  line.  Also,  since  the  main  radiation  comes 
from  the  dipole,  the  slot  length  can  be  chosen  far  below  its  resonant  length. 
The  equivalent  circuit  of  a  slot  coupled  dipole  seen  by  the  microstrip  line  is  a 
scries  clement.  An  example  of  this  equivalent  circuit  as  a  function  of  dipole 
length  for  two  different  slot  lengths  is  shown  in  Fig.  8-13.  It  is  observed  that 
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the  increase  of  slot  length  will  decrease  the  dipole  length  to  obtain  resonance. 
This  is  because  when  the  dipole  length  decreases,  resonance  will  be  due  to  the 
resonance  of  the  slot,  whereas  when  the  slot  length  decreases,  the  resonance 
will  be  due  to  the  resonance  of  a  center-fed  dipole.  It  is  also  observed  that  the 
resonant  resistance  is  larger  for  larger  slot  length. 
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ground  plane 


Figure  8*12.  A  microsUip  fed  slot  coupled  dipole. 

Microsirip  line  width:  w  ^  ,  slot  length  :  L^, 
slot  width:  w  ,  dipole  length:  L  and 

9  0 

dipole  width:  w . 
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Figure  8-13.  Normalized  impedaiKC  of  a  micros: 
versus  dipole  length.  e|  «  C}  =  2.54, 
=  1  mm,  s  4.45  mm,  -  h,. 
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Chapter  IX 


Design  of  Transversely  Fed  EMC  Microstrip  Dipole  Arrays 
Including  Mutual  Coupling 


9.1  Introduction 

Oilman  introduced  a  class  of  electromagneticaly  coupled  (EMC)  dipole  an¬ 
tennas  [81],  [56].  The  advantages  of  EMC  dipoles  are  greater  bandwidth, 
higher  efficiency  and  more  easily  matched  to  the  feed  lines,  when  compared  to 
classically  fed  printed  antennas.  Based  on  the  transmission  line  circuit  model, 
Oilman  and  Huebner  [56]  built  a  longitudinal  EMC  dipole  array.  Later, 
Elliott  and  Stern  developed  a  rigorous  design  theory  to  include  the  effects  of 
mutual  coupling  which  successfully  predicted  the  array  performance 
[59],  [60].  An  efficient  way  to  obtain  the  design  curves  was  reported  later 
[42],  [87]  based  on  solving  a  Pocklington  type  integral  equation  using  the 
method  of  moments.  The  EMC  collinear  dipole  is  ideally  suited  to  a  corporate 
feed,  and  elements  of  this  type  can  be  arranged  in  circular  as  well  as  rectan¬ 
gular  grids. 

Another  dipole  antenna  of  the  Oilman  type  is  the  EMC  transverse  dipole, 
as  shown  in  Fig.  9-1,  where  a  dipole  is  oriented  transverse  to  an  embedded 
microsirip  line.  A  string  of  these  dipoles  above  a  common  microstrip  line  be¬ 
comes  a  linear  array.  Depending  on  dipole  spacing,  one  can  obtain  standing 
wave  arrays  or  traveling  wave  arrays.  A  family  of  these  linear  arrays  becomes 
a  planar  array.  Current  e.xcitation  on  the  dipole  is  governed  by  the  amount 
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Figure  9-1.  The  electromagnetically  coupled  transverse  dipole. 


of  offset  and  the  dipole  length.  If  the  dipole  straddles  symmetrically,  no 
excitation  of  the  dipole  occurs.  Light  coupling  from  line  to  dipole  can  be 
achieved  through  a  slight  lateral  displacement  of  the  dipole  from  its  balance 
position.  Because  of  this  light  coupling,  large  arrays  arc  feasible  [SSj. 

Both  theoretical  and  experimental  studies  of  the  EMC  transverse  dipole  are 
reported  recently  [77].  In  this  chapter,  a  design  technique  which  includes 
mutual  coupling  is  developed  for  the  EMC  transverse  dipole  arrays.  Two  de¬ 
sign  equations  will  be  introduced  in  Section  9.2.  The  methods  for  generating 
design  curves  will  be  discussed  in  Section  9.3.  A  design  example  will  be  given 
in  Section  9.4  together  with  the  experimental  results.  A  numerical  verification 
of  the  design  by  solving  the  boundary  value  problem  of  the  whole  array  system 
will  also  be  provided. 

9.2.  Two  Design  Equations. 

For  the  transverse  EMC  dipole  under  consideration,  it  can  be  shown  from 
image  theory  that  the  scattering  off  the  dipole  is  symmetric.  In  other  words, 
the  forward  and  backward  .scattering  coefficients  arc  the  same.  Therefore,  in 
terms  of  the  transmission  line  equivalent  circuit,  the  dipole  can  be  approxi¬ 
mated  as  a  shunt  clement  with  respect  to  the  feed  [61].  Each  dipole  in  the 
array  environment  can  then  be  modelled  as  a  two  port  network  and  the  whole 
.system  is  a  linear  bilateral  network  [59]  .  Therefore,  one  can  write 

N 

I,  -  Z 
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as  a  set  of  equations  connecting  the  transmission  line  mode  voltages  and  cur¬ 
rents  at  each  reference  port  [61  ].  The  active  admittance  of  each  dipole  seen 
by  the  feed  line  can  be  defined  from  Eq.  9.1  by 


=  }’  +  Y 

* nn  ^  ‘  n 


where 


with  the  prime  on  X  indicating  the  term  n  =  m  is  excluded.  is  the  self- 
admittance  and  is  referred  to  as  the  mutual  admittance  due  to  the  mutual 
coupling  between  each  dipole.  Since  a  linear  system  is  assumed,  the  current  in 
each  dipole  can  be  written  as 
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where  is  the  current  of  the  nth  dipole  for  a  given  mode  voltage  without  the 
prt'sence  of  other  dipoles  and  is  the  current  of  the  nth  dipole  due  to  the 
current  in  the  mth  dipole.  The  isolated  dipole  current  is  a  function  of  mode 
voltage,  dipole  length  and  offset  which  can  further  be  written  as 
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where  l„)  is  a  coefficient  function  relating  the  isolated  dipole  current  to  its 
self-admittance. 

The  two  design  equations  can  be  summarized  from  the  above  derivations 
by 
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It  is  noted  that,  for  the  EMC  tranverse  dipole,  the  current  phase  variation  in 
an  isolated  dipole  is  quite  large,  typically  5  to  10  degrees;  while,  in  contrast,  the 
current  due  to  mutual  coupling  has  a  small  phase  variation.  Therefore,  in  the 
design,  the  cocfilcient  function  is  not  suitable  for  relating  mutual  current  (/„,„) 
and  mutual  admittance  (}'„,„).  One  can  use  the  mutual  current  term  (/„,„)  di¬ 
rectly  in  Eq.  9.6. 

The  fundamental  design  problem  is  now  obvious.  For  a  given  design  goal 
(radiating  current  in  each  dipole  /"‘O  one  wishes  to  find  N  set  of  (5„,  Q  such 
that  not  only  Eq.  9.6  is  satisfied  but  also  the  active  admittance  seen  by  each 
feed  line  is  what  was  prescribed.  The  definition  of  the  radiating  current  in 
each  dipole  may  depend  on  the  design  goal.  For  e.xample,  if  one  wishes  to  de¬ 
sign  a  specific  pattern  in  the  H  plane  (perpendicular  to  the  dipole),  the  M 
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should  be  defined  as  the  current  collapsed  in  the  feed  line  [61  ].  In  other 
words, 

(9-8) 

In  the  two  design  equations  in  Eqs.  9.6  and  9.7,  Y ImJK  and  'Y„„  can  be 
determined  by  the  method  of  moments.  This  issue  will  be  discussed  in  the  next 
section. 

Suppose  that  all  four  functions  arc  known.  Further  computations  are  still 
required  to  find  the  dipole  lengths  and  offsets.  Since  only  relative  currents  in 
the  dipoles  are  meaningful,  one  can  arbitrarily  choose  a  dipole,  say  the  nth, 
with  length  l„  and  offset  s„.  For  the  moment,  assume  that  no  mutual  coupling 
exists  such  that  the  left  hand  side  of  Eq.  9.6  can  be  determined  according  to 
the  design  goal  (desired  currents  in  the  dipoles).  One  now  can  use  the  first 
design  equation  (Eq.  9.6)  to  find  N-1  dipole  lengths  and  offsets.  This  proce¬ 
dure  requires  that  a  two-variable  nonlinear  equation  be  solved  N-1  times.  To 
avoid  the  stability  and  .solvability  problem  of  this  nonlinear  equation,  the 
conjugate  gradient  method  [89]  cun  be  u.scd  to  provide  optimized  solutions. 
Even  if  the  above  procedures  are  completed,  a  few  iterations  by  changing  the 
dipole  length  or  offset  of  the  first  selected  dipole  are  required  to  provide  the 
prescribed  input  impedance.  Now  the  design  data  is  what  one  should  obtain 
if  no  mutual  coupling  e.xists.  To  include  the  effect  of  mutual  coupling,  one  can 
use  the  present  design  data  to  compute  the  mutual  admittance  and  mutual 
currents  and  go  back  to  the  two  design  equations  repeating  the  above  iter- 
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aliens.  The  whole  procedure  is  iterated  until  the  convergence  of  the  design 
data  is  found. 

9.3.  Discussion  of  the  Method  of  Moments  Solution 

In  order  to  make  an  accurate  design  possible,  information  about  the  inter¬ 
action  between  dipoles  as  well  as  dipole  coupling  to  the  feed  line  is  required. 
The  method  of  -moments  provides  a  rigorous  and  accurate  solution.  Integral 
equations  for  the  EMC  transverse  dipole  can  be  written  as 

£^.v,.v)  =  I  J  C^,4'V,y)<iV  <*''  +  II y)<ix’  4'  (9.9) 

and 

j  Cy^ji'\x',y)d.x'  4'  +  ||c,^.^V./)*'4''  (9.I0) 

where  is  the  current  in  the  microstrip  line  and  is  the  cur- ' 

rent  in  the  dinole.  The  functions  Gxp  Gyx  and  Gyy  arc  the  dyadic  Green's 
function  components.  A  nice  way  of  modeling  the  feed  line  is  to  use  a  finite 
but  long  microstrip  line  with  a  5-gap  generator  placed  far  from  the  linc-dipolc 
coupling  region  [42],  [43].  When  the  combination  of  piecewise  sinusoidal 
and  Maxwell  current  basis  functions  is  used  in  the  method  of  moments 
[42],  [43]  followed  by  the  Galerkin  procedure,  the  matrix  equations 

[z«]  [/»]  =  [£■»]  (’ll) 
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can  be  obtained.  The  e.xcitation  column  vector  E  has  components  £*  =  -1 
when  the  <5  gap  source  is  located  at  the  center  of  the  kth  basis  function,  and 
£„  =  0  any^vhcrc  else.  The  impedance  matri.x  elements  are  in  the  form  of 


Z„,„  =  /1(A)]  JMx’  K) 

••o  •'0 

Xy)  cos(y^.v)  cos{XyAy)dXd<i> 


(9.12) 


when  the  mth  and  the  nth  basis  functions  arc  both  on  the  dipole  or  on  the  line, 
and 

rnfl  roo  _ 

Jq  Jq  _  (7.13) 

Xy)  sin(/.^.Y)  sin{XyAy)dXd(p  , 

otherwise.  The  functions  7„{Xx,Xy)  and  J„{Xx,Xy)  arc  the  Fourier  transforms 
of  the  expansion  and  testing  functions  respectively.  (Ax,  Ay)  is  the  displace¬ 
ment  vector  of  two  basis  function  centers, 

Xy.  =  ycos0  (9.14) 

and 

Xy  =  X  .sin  (/>.  (9.15) 


The  function  J(X)  and  h{X)  arc  related  to  the  Hertz  potential  described  in  Chap. 
II.  The  unknown  currents  in  the  feed  line  or  the  dipole  can  be  obtained  by 
matrix  inversion.  As  a  result,  one  can  u.sc  the  unimode  transmission  line  the¬ 
ory  to  deduce  the  circuit  information  from  the  method  of  moments  solution  of 
the  current  in  the  feed  line.  This  procedure  involves  finding  the  current  maxi- 
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mum.  minimum  and  their  positions.  Detailed  procedures  of  this  are  shown  in 
C42].  It  is  noted  that  in  order  to  find  the  dipole  equivalent  admittance,  the 
dipole  can  be  placed  a  half  electrical  wave  length  from  the  line  end  such  that 
the  stub  admittance  will  not  be  included  in  the  input  admittance  observed  by 
the  feed  line. 

One  of  the  features  of  using  b  gap  source  is  that  the  mode  voltage  of  the  line 
changes  with  the  change  of  the  dipole  length  or  offset;  therefore  care  must  be 
taken  to  find  this  mode  voltage,  since  as  shown  in  Eq.  9.6,  dipole  current  is 
proportional  to  the  mode  voltage  and  only  their  ratio  is  useful  in  the  design. 

Another  feature  of  using  the  moment  method  in  the  array  design  is  that  the 
information  about  the  currents  in  the  dipoles  can  be  obtained  from  the  nu¬ 
merical  process.  This  aspect  is  particularly  helpful,  since  one  can  use  the 
dipole  current  directly  to  design  for  the  desired  excitation  instead  of  using  the 
cqui\alcnt  circuit  of  the  dipoles.  This  will  be  discussed  further  in  the  next 
section.  Other  issues  in  this  array  design  arc  how  the  mutual  coupling  infor¬ 
mation  can  be  .separated  from  that  of  the  self  term  and  how  this  can  be 
achieved  without  involving  the  whole  system  at  the  same  time.  In  order  to 
solve  these  problems,  certain  assumptions  arc  necessary.  It  is  assumed  that 
the  self  admittance  and  self  current  to  mode  voltage  ratio  will  not  change  with 
or  without  the  presence  of  other  dipoles,  and  that  the  mutual  coupling  between 
any  two  dipoles  is  unaffected  by  the  rest.  These  two  assumptions  are  good  if 
mutual  coupling  is  not  too  strong  [61  ],  which  is  usually  true  for  practical  ar¬ 
rays. 
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The  method  of  computing  mutual  coupling  of  dipoles  individually  fed  by  a 
microstrip  line  has  been  discussed  in  [87].  One  of  the  features  of  the  array 
considered  here  is  that  the  dipoles  arc  scries  fed  by  transverse  microstrip  lines. 
The  computation  of  mutual  coupling  in  this  case  requires  a  different  approach 
from  [87].  To  find  the  mutual  coupling  information,  one  can  consider  two 
dipoles  fed  by  a  microstrip  line  and  follow  a  similar  numerical  method  as  for 
an  isolated  dipole  case  except  for  the  additional  computation  of  dipole  to 
dipole  reaction.  An  asymptotic  c.xtraction  technique  together  with  point 
source  approximation  has  been  developed  in  [44]  to  compute  efficiently  and 
accurately  the  reaction  of  two  dipoles.  After  the  matrix  inversion,  the  solution 
of  the  current  in  the  line  or  dipoles  provide  the  total  active  admittance  as  well 
as  active  currents  of  the  dipoles.  For  the  nth  and  mth  dipole  with  resonant 
spacing,  the  method  of  moments  allows  one  to  compute  the  total  active 
admittance 
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The  -  or  +  signs  depend  on  whether  the  dipole  spacing  is  an  odd  or  even  in¬ 
teger  of  a  half  guide  wavelength.  From  Eq.  9.16,  if  the  self  admittance  of  each 
dipole  is  known,  the  mutual  admittance  can  be  determined.  The  active  cur¬ 
rents  in  the  dipole  and  can  also  be  obtained  numerically  and  can  be  de¬ 
scribed  as 
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From  Eqs.  9.17  and  9.18  together  with  the  information  of  the  current  excited 
in  an  isolated  dipole,  the  mutual  current  can  be  obtained. 


9.4.  A  Design  Example 

The  previous  discussion  of  the  design  theory  applies  either  to  a  linear  or  a 
planar  array.  Here,  a  se\en  elements  standing  wave  linear  array  will  be  de¬ 
signed  to  illustrate  the  design  technique.  The  geometry  is  shown  in  Fig.  9-2. 
A  sum  pattern  in  the  H  -plane  with  a  -20  dB  side  lobe  level  was  prescribed  for 
this  array.  The  printed  dipoles  arc  scries  fed  by  a  50  D  microstrip  line  em¬ 
bedded  in  the  middle  of  the  substrate  of  thickness  0.16  mm  and  permittivity 
2.17.  The  design  frequency  is  8.3  GHz  and  element  spacing  is  chosen  to  be  one 
guide  wavelength.  All  dipoles  have  the  same  width  of  1  mm  and  the  offsets 
and  lengths  are  to  be  found.  The  design  curves  for  the  self  admittance  Y„„  and 
the  coefficient  function  f„as  a  function  of  offset  and  length  obtained  from  the 
method  of  moments  solution  are  shown  in  Figs.  9-3  and  9-4  respectively.  It  is 
found  that,  for  the  EMC  transverse  dipole,  many  basis  functions  arc  required 
to  obtain  adequate  convergence.  To  obtain  each  data  point,  19  expansion 
modes  arc  used  in  the  dipole  and  piecewise  sinusoidal  modes  of  size  of  0.04 
guide  wavelength  arc  used  in  the  line.  It  is  observed  from  Fig.  9-4  that  for 
dificrent  dipole  offsets  and  lengths,  the  phase  of  f„  is  not  a  constant. 
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Figure  9-2.  A  seven  element  linear  standing  wave  array. 
Elements  are  spaced  by  a  guide  wave  length. 
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This  implies  that  even  for  resonant  spacing,  to  have  in-phase  excitation,  the 
dipole  can  not  be  self  resonant.  To  obtain  a  perfect  match,  the  stub  length 
Ax  in  Fig.  9-2  can  be  suitably  adjusted  to  tune  out  the  total  active  susceptance. 

The  sampled  data  arc  used  to  construct  the  data  bank  such  that  for  a  given 
offset  and  length  the  function  value  can  be  obtained  through  a  two  dimen¬ 
sional  interpolation  routine.  Mutual  coupling  between  two  dipoles  is  a  func¬ 
tion  of  dipole  lengths  and  offsets  for  a  fixed  spacing.  It  is  found  that  mutual 
coupling  is  not  sensitive  to  a  small  change  of  dipole  length.  Also  from  the 
results  of  no  mutual  coupling,  it  is  found  that  the  lengths  of  all  the  dipoles  are 
different  within  0.2%.  Therefore,  in  the  mutual  coupling  computation,  the 
dipole  lengths  arc  fixed.  The  mutual  admittance  and  mutual  current  as  a 
function  of  offsets  for  one  guide  wavelength  spacing  and  fixed  dipole  length 
arc  shown  in  Figs.  9-5  and  9-6  respectively.  A  similar  procedure  can  be  fol¬ 
lowed  for  a  two  wavelength  spacing,  and  so  on.  As  a  result,  the  interpolation 
or  extrapolation  method  can  provide  the  mutual  coupling  information  for  any 
dipole  spacing. 

The  design  data  after  a  few  iterations,  including  the  self  and  mutual 
admittances  arc  shown  in  Table  I.  The  clement  spacing  is  23.6  mm  and 
Ax  =  1 1  mm.  It  is  .seen  that  the  mutual  admittance  is  more  than  25%  of  the 
total  active  admittance.  Therefore,  it  is  concluded  that  even  at  one  wave 
length  spacing,  the  effect  of  mutual  coupling  should  not  be  ignored.  To  pro¬ 
vide  a  confident  check  of  the  design  data,  the  method  of  moments  is  applied 
to  the  seven  clement  linear  array.  The  results  of  active  admittance  and  cur¬ 
rent  in  each  dipole  arc  shown  in  Table  11  together  with  the  results  from  the 
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synthesis  technique.  It  is  observed  that  the  current  amplitude  agrees  within  1 
to  2%  and  the  phase  agrees  in  ±  1  degree.  The  admittance  comparison  is  also 
good. 

The  antenna  array  was  built  on  a  10-in  square  Duroid  board  [90].  The 
measured  return  loss  from  the  feed  line  is  shown  in  Fig.  9-7.  The  bandwidth 
of  this  array  is  about  3.75%.  The  frequency  for  a  perfect  match  is  found  to 
be  8.14  GHz  which  deviates  160  MHz  (2%)  from  the  design  frequency.  This 
result  is  rather  good,  considering  that  this  array  is  very  sensitive  to  the  toler¬ 
ances  of  the  device  parameters.  Both  the  desired  and  measured  radiation 
patterns  in  H-plane  arc  shown  in  Fig.  9-8.  It  is  observed  that  the  main  beam 
as  well  as  the  first  side  lobe  agree  well  with  the  design  criteria.  Other  side 
lobes  in  the  measured  pattern  arc  a  little  too  high,  which  is  probably  due  to  the 
finite  size  of  the  ground  plane. 
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Design  data  of  n  seven  element  linear  array 

No. 

offset  in  mm 

I 

dipole  length  tn  mm 

ymumol 

1 

Q.5()91 

12.72.1 

0.0580  +  j  0.0914 

0.015  +  j  0.001 

2 

0.1425 

12.74.1 

0.0825  4  j  0.0519 

0.018  +  j  0.006 

0.9467 

12.716 

0.14.50  +  j  0.0020 

0.051  +  j  0.006 

4 

1.0200 

12.701 

0.1769  -  j  0.0215 

0.060  +  j  0.005 

5 

0.9467 

12.716 

0.1450  +  j  0.0020 

0.051  +  j  0.006 

6 

0.7425 

12.741 

0.0825  +  j  0.0519 

0.018 +  j  0.006 

7 

0.5697 

12.721 

0.0580  +  j  0.0914 

0.015  +j  0.001 

}■„,  =  1.02  +  y0.002  (ic.sull  from  iterations) 
=  1 .08  +  y0.006  (result  from  I.E.S.) 


TABLE  II 


Design  check  through  an  integral  ^uation  .solution  (I.E.S.) 

No. 

desired  current 

I.E.S.  current 

1 

1.0000  <  0.0 

1.0000  <  0.0 

2 

1.2751  <0.0 

1.264  <-1.2 

1 

1.6810  <  0.0 

1.6.50  <-0.2 

4 

1.8.151  <  0.0 

1.824  <  1.4 

.5 

1.6810  <  0.0 

1.6.54  <  0.4 

ft 

1.2751  <  O.O 

1.266  <-0.8 

7 

l.(X)(8)  <  0.0 

().9W  <  0.8 

The  unit  III  the  ph;ise  of  cuiiciit  is  in  dcgiee 


yang  2  vr/2/%B  Sll  dB-LOSS 
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Figure  9-7.  Measured  return  loss  of  the  designed  antenna  array. 


scan  angle  in  the  1 1-  plane 


Figure  9-8.  Compari.son  between  the  dc.sircd  and  measured  antenna 
pattern  in  the  I  i-plane. 
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Chapter  X 


Conclusion 

In  this  dissertation,  a  generalized  frequency  dependent  method  has  been  de¬ 
veloped  to  characterize  integrated  circuit  discontinuities  and  printed  circuit 
antennas.  The  analysis  has  taken  into  account  all  the  physical  effects  includ¬ 
ing  radiation,  surface  waves  and  higher  order  modes.  The  method  of  moments 
has  been  used  to  soKc  coupled  integral  equations.  In  the  procedure,  a  mode 
e.xpansion  mechanism  using  the  combination  of  the  entire  domain  traveling 
wave  modes  and  the  piecewise  sinusoidal  subdomain  modes  has  been  devel¬ 
oped.  This  scheme  is  particularly  useful  for  the  analysis  of  coupled  infinite 
lines. 

The  geometries  of  microstrip  open-end  and  gap  discontinuities  in  two  layer 
structures  are  analyzed.  The  fringing  effect  and  radiation  and  surface  wave 
losses  as  a  function  of  various  device  parameters  are  studied.  The  resuits  of 
this  study  have  been  obtained  through  a  careful  convergence  test  and  com¬ 
pared  against  the  quasi-static  results  in  low  frequency  range  with  excellent 
agreement.  A  short-end  slot  line  discontinuity  in  a  one-side  or  a  sandwiched 
structure  is  also  studied.  The  advantages  of  using  a  slot  line  sandwich  have 
been  demonstrated. 

Line  to  line  transitions  which  arc  increasingly  important  in  VLSI  and 
monolithic  circuits,  have  been  investigated.  The  transition  circuits  studied 
include  microstrip-slot  line  transition,  proximity  coupled  collincar  and  trans- 


verse  microstrip  transitions.  The  analysis  of  microstrip-slot  line  transition 
which  includes  the  coupling  between  microstrip  and  slot  line,  reveals  the  ver¬ 
satility  of  the  developed  numerical  methods.  The  analysis  can  be  easily  im¬ 
plemented  in  a  computer  aided  design  of  both  microstrip  circuits  and  slot  line 
circuits.  The  analysis  of  proximity  coupled  microstrip  transition  has  shown 
that  these  transitions  arc  broadband  and  provide  a  wide  range  of  coupling  co¬ 
efficients,  and  have  potential  importance  in  the  design  of  microwave  filters  and 
couplers.  The  experiment  for  the  case  of  a  trans\crsc  transition  has  been 
performed  to  verify  the  analysis. 

The  method  for  characterizing  integrated  circuit  discontinuities  is  general 
enough  to  analyze  a  class  of  printed  circuit  antenna  feeding  structures.  These 
structures  include  a  microstrip  fed  slot,  a  slot  line  fed  dipole,  an  EMC  trans¬ 
verse  dipole  and  a  microstrip  fed  slot  coupled  dipole.  The  main  features  and 
potential  applications  for  each  antenna  architecture  in  monolithic  phase  arrays 
have  been  discussed. 

Design  techniques  and  procedures  for  microstrip  dipole  arrays  transversely 
fed  by  proximity  coupled  microstrip  lines  have  been  presented.  Two  design 
equations  which  include  the  effects  of  mutual  coupling  arc  developed  and  the 
corresponding  design  cur\es  arc  obtained  by  a  rigorous  integral  equation  sol¬ 
ution.  A  seven  clement  standing  wave  linear  array  has  been  designed  to  illus¬ 
trate  the  developed  design  procedures.  The  design  data  is  checked  by  a 
complete  integral  equation  .solution  of  the  array  with  c.xccllcnt  agreement.  The 
radiation  pattern  and  input  impedance  measurements  are  also  compared  with 
theory. 
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The  integrated  circuit  components  studied  in  this  dissertation  are  either 
isolated  line  discontinuities  or  proximity  coupled  transition  from  line  to  line. 
The  circuit  components  arc  slender  strips  or  slots  in  rectangular  shape.  These 
considerations  lead  to  the  approximation  of  a  one  current  component  in  the 
strip  and  a  one  electric  field  component  in  the  slot.  This  simplified  analysis 
has  greatly  reduced  the  CPU  time  without  losing  desired  accuracy.  With  a 
straightforward  (but  tedious)  modification  of  the  present  analysis  which  con¬ 
siders  both  current  or  field  components  in  the  circuit,  one  is  able  to  character¬ 
ize  microstrip  and  slot  junction  discontinuities.  The  extension  of  the  present 
research  to  coplanar  wa\cguidc  or  coplanar  coupled  line  is  also  possible.  The 
analysis  of  the  printed  antenna  elements  can  be  used  directly  to  design 
monolithic  antenna  arrays  without  considering  mutual  coupling.  The  method 
of  incorporating  mutual  coupling  in  the  array  designs  has  also  been  demon¬ 
strated  with  an  EMC  transverse  dipole  array  as  an  example. 
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Appendix  A 


The  pertinent  parameters  in  Eqs.  2.9-2.14  are  expre.ssed  as 
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for  b  <'  z  ^ 

M^')  =  ^4(^)sinh  <72(^  -  b)  +  A^X)co'nh  qj^z  -  b)  (/i.5) 

for  z^  <  z  <  htf2{/.)  is  .similar  to  Equation  (A.5)  c.xccpi  z  and  Zj  arc  inter¬ 


changed  due  to  the  reciprocity. 
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Appendix  B 


The  pertinent  parameters  in  Eqs.  2.28-2.33  are  e.xpressed  as 
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for  0  <  r  <  Zj,f\{X)  is  similar  to  Eq.  (B.3)  except  z  and  z,  are  interchanged 
due  to  the  reciprocity. 
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(5.9) 
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Appendix  C 


A- 

Thc  pertinent  parameters  in  Equations  (2.85)-(2.90)  are  expressed  as 
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Appendix  D 


and  J^')  is  shown  in  Eq.  (3.5).  The  quadruple  integral  in  Eq.  D.l  can  be  fur¬ 
ther  reduced  to  a  double  integral  with  integrand  containing  the  Fresnel  func¬ 
tions  if  kg  =  kg\  .  Since  the  integration  range  is  usually  small,  the  integration 
converges  quite  easily. 


Appendix  E 


The  dominant  term  of  the  asymptotic  form  of  Eq.  3.30  is 

MQk.ykf,  r'/2 

£  =  — : -  — — ~ —  ( cos  cos  cos  dX  d(f>  (£.  1 ) 

J()  It  a;  -  kl  ' 

E  when  transformed  back  to  Fourier  domain  can  be  expressed  as 

£  =  0.5  £i(</i  +  nd^)  +  0.5  £|(</i  —  —  cos  kgd^  E^(nd{) 

+  0.5  Eid^  +  nd^)  +  0.5  E^di  -  nd^)  -  cos  k^^  E^nd^) 

where 
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"  '■>  -  *»)  +  4  ■ 


and 


4  =  (£.5) 

The  integrations  in  Eqs.  E.3  and  E.4  both  converge,  however,  in  order  to  ob» 
tain  more  accurate  and  efficient  computation,  further  analysis  may  be  re¬ 
quired.  For  the  Xx  integration  in  Eq.  E.3,  the  integration  contour  in  a  complex 
plane  can  be  deformed  (referring  to  Fig.  E-1)  along  the  branch  cut.  This  will 
help  the  integration  decay  exponentially  when  the  cos  argument  is  not  zero. 


Using  the  branch  cut  integration  and  a  secant  transformation.  £,(»)  can  be 
written  as 
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The  con\  crgcncc  of  the  integration  in  Eq.  E.4  can  be  improved  for  a  7^  0  by 
deforming  the  integration  contour  in  the  complex  plane.  The  new  contour  is 
shown  in  Fig.  E-2.  The  final  form  of  the  integration  in  Eq.  E.4  can  be  written 
as 


(£'.7) 


where 

/.(«)  =  yjA~  -  xj  +  ju  (£.8) 


It  should  be  noted  that  the  .singularities  in  either  Eq.  E.I  or  E.2  arc  ignored 
since,  as  pointed  out  in  Chap.  Ill,  these  singularities  arc  removable. 
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branch  cut 


Figure  E-1,  Integration  contour  of  Eq.  E.3. 
branch  point  at 


Figure  E-2.  Integration  contour  of  Eq.  E.4. 
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Abstract-  A  full-wave  spectral-domain  analysis  is  applied 
to  the  characterization  of  multi-port  microstrip  discontinu¬ 
ities.  This  approach  employs  the  moment  method  to  find 
the  currents  in  the  microstrip  circuits  and  subsequently,  the 
scattering  parameters  of  the  junctions.  In  this  approach, 
all  the  physical  cITects  are  considered,  including  radiation 
and  surface  waves.  The  numerical  results  for  a  tee  and  a 
cross  junction  are  presented  and  agree  well  with  the  quasi¬ 
static  values  at  low  frequencies.  The  S  parameters  of  a 
tee  junction  are  further  compared  against  the  measured  re¬ 
sults  with  e.vccllent  agreement.  The  utilization  of  a  shaped 
T-junction  as  a  hroad-band  equal-power  divider  is  also  dis¬ 
cussed. 

I.  Introduction; 

The  full-wave  analysis  which  deals  with  microstrip 
discontinuitic's  in  an  open  geometry  has  been  applied  to  a 
variety  of  problems.  This  approach  based  on  the  moment 
method  .solution  of  an  exact  integral  equation  involves  the 
compulation  of  a  continuous  plane  wavcspectrum  such  that 
the  effects  of  radiation,  .surface  waves  and  the  higher-order 
modes  are  included.  This  full-wave  analysis  has  been  ap¬ 
plied  to  microslrip  open-ends  and  gaps  (1-3),  steps  (4),  and 
bends  (-i-Tj. 

From  the  leview  of  the  past  work,  one  finds  the  full 
wave  analysis  u))  to  now  is  limited  to  two-port  structures. 
In  this  paper,  a  full-wave  analysis  up  to  four  ports  is  pre¬ 
sented.  The  spectral  domain  dyadic  Green's  functions  are 
adopted  in  electric-field  integral  equations  (EFIE)  to  han¬ 
dle  double-layer  substrate  (troblems.  Uolh  longitudinal  and 
transverse  current  components  on  the  microslrip  are  taken 
into  account  and  are  the  solution  of  the  method  of  mo¬ 
ments.  In  Section  ll,  the  method  of  moments  formulation 
of  the  EFIE  is  brielly  outlined.  Mode  expansion  utilizing 
the  combination  of  semi-infinite  line  modes  on  transmis¬ 
sion  lines  and  piecewise  sinusoidal  basis  functions  on  the 
vicinity  of  the  discontinuity  are  also  discussed.  In  Section 
HI,  numerical  results  of  the  scattering  parameters  of  var¬ 
ious  nricrostrip  discontinuities  such  as  basic  tee  and  cross 
junctions  as  well  as  a  shaped  T-junction  are  discussed  and 


compared  with  available  measurements  and  quasi-static  re¬ 
sults. 

n.  Analysis: 

Microstrip  discontinuities  can  be  looked  upon,  from 
a  circuit  point  of  view,  as  N  transmission  lines  (N  ports) 
jointed  in  a  common  region.  The  modeling  involves  finding 
the  current  distribution  not  only  in  the  junction  region  but 
also  in  the  N  semi-infinite  microstrip  transmission  lines. 
The  approach  to  solve  the  problem  is  based  on  the  moment 
method  solution  of  the  exact  integral  equations. 

A  generic  four-port  microslrip  discontinuity  is  pre¬ 
sented  in  Fig.  1.  Pour  semi-infinite  microslrip  transmission 
lines  which  extend  to  d:  oo  in  x  or  y  direction  are  jointed 
in  a  common  block  (dash  line  box  in  Fig.  1 ).  The  widths 
of  these  four  transmission  lines  arc  not  necessary  the  same. 
The  planar  configuration  of  the  microstrip  discontinuity  in¬ 
side  the  common  block  is  a  state  of  the  art  to  design  the 
desired  performance  of  this  junction  in  a  specified  frequency 
range. 

For  microstrip  junction  problems,  the  concept  of  a 
module  can  be  used.  A  module  encloses  the  region  at  or 
near  the  junction  where  higher  order  modes  are  generated. 
The  region  otherwise  consists  of  purely  uniform  transmis¬ 
sion  lines.  The  currents  inside  the  module  are  expanded 
in  terms  of  piecewise  sinusoidal  basis  functions,  while  the 
currents  outside  the  module  are  uniform  transmission  line 
currents  (semi-infinite  mode  SIM). 


Fig.  1  A  generic  structure  of  a  four  port  microetrip  dismniinuiiy 


(This  research  was  supported  by  U.S.  Army  Research  Grant  DAAL  03-86-K-0090 


ir  a  local  coordinate  system  is  used  at  cacli  line  and 
the  excitation  is  in  the  it'''  port,  then  away  from  the  junc¬ 
tions,  the  current  in  the  n"'  port,  with  longitudinal  com¬ 
ponent  oriented  in  x  direction,  is 

Jj-  =  -  r„  (1) 

while  the  current  in  the  port  in  x  direction,  with  p  /  n, 
is 

=  (2) 

In  addition,  the  current  in  the  q"'  port  in  y  direction  is 


Sn 


(ji  is  a  sign  index  function  of  k"'  port. 


a  = 


( 


+l 

-1 


;  k'''  line  extends  to  -I-oo  in  x  or  y  direction 
;  k"'  line  extends  to  — oo  in  x  or  y  direction 

(4) 


where  k  could  be  p,  q  or  n.  /*  and  /?*  are  the  prc-calculated 
current  transverse  dependence  and  the  propagation  con¬ 
stant  on  the  k'''  microstrip  transmission  line,  respectively. 
It  is  noted  that,  in  terms  of  scattering  parameters,  I’*  is  5t„. 
This  aspect  describes  a  uniipie  feature  of  the  approach. 


The  basis  function  of  the  module  is  chosen  as  a 
piecewise  sinusoidal  (PWS)  function  in  the  longitudinal  di¬ 
rection  (the  direction  of  current  flow)  and  a  pulse  function 
in  the  transverse  direction.  Mathematically,  the  current 
inside  the  module  can  be  expressed  as 


nsl 


M 


E  l?l'7W(y)\  y  (5) 

.ma| 


The  closed  form  of  the  Fourier  transform  of  the  basis  func¬ 
tion  can  be  found. 


nearly  Galcrkin  method  is  applied  to  transform 
the  integral  eiinalions  into  a  matrix  equation.  Inside  the 
module,  the  tc'sting  functions  arc  in  the  same  form  as  the 
current  basis  functions.  The  testing  function  for  the  semi- 
infinite  line,  which  is  chosen  as  a  pitxrewisc  sinusoidal  func¬ 
tion  in  direction  of  current  (low  and  with  the  same  trans¬ 
verse  dependence  as  the  transmission  line  current,  is  ap¬ 
plied  adjacent  to  the  module.  Finally,  the  matrix  cipiation 
is  in  the  form  of 


lCn-1 

K, 

1 

.K 

y)/r  1  \yyy 

1^:1 


(6) 


Each  element  in  the  submatrices  of  [Z\  and  [K„)  is  the  re¬ 
action  between  the  basis  function  and  the  testing  function 
and  is  in  the  form  of  a  double  integration  in  a  spectral 
domain.  The  superscripts  in  [Z]  indicate  the  orientations 
of  the  corresponding  testing  and  basis  functions;  and  the 
subscripts  represent  their  locations.  For  instance, 

(7) 

is  the  reaction  between  the  basis  function  (current  in 
the  X  direction),  with  the  testing  function  /"'  (current  in 
the  y  direction).  [K-„]  is  in  the  same  mathematical  form  as 
[Z]  except  that  the  basis  function  is  the  semi-infinite  mode 
of  the  incident  wave. 

III.  Results  and  Discussions: 

In  this  research,  several  microstrip  discontinuities 
are  analyzed  by  the  method  described  above.  The  data 
generated  by  TOUCHSTONE  version  1.7  ,  which  are  es¬ 
sentially  quasi-static  results,  are  also  presented  for  com¬ 
parison. 
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(d).  the  substrate  structure 

Fig.  2  Layout  of  a  variety  of  juiictioiis 


A  basic  T-junction  with  three  identical  semi-inrinite 
transmission  lines  is  sliown  Fig.  2(a).  The  numerical  results 
shown  in  Fig.  3  and  4  arc-converged  within  0.2dB  in  mag^ 
nitude  and  ‘2.b°  in-  phase.  The  results  of  the  magnitude  of 
the  scattering  coefficients  are  compared  against  with  mea¬ 
surements  (Sj  and  the  quasi-static  values  (TOUCHSTONE 
data),  and  are  shown  in  Fig.  3.  It  it  seen  that  the  present 
full-wave  results  are  iii  excellent  agreement  with  the  mea¬ 
sured  data,  but.agree  well  only  in  the  low  frequency  range 
with  the  TOUCHSTONE  results.  In  the  high  frctpiency 
region,  the  unequal  power  transmitted  on  S21  and  S31  ob¬ 
served  in  both  theory  and  measurement  is  more  significant 
than  what  TOUCHSTONE  predicts.  It  is  noted  that  the 
TOUCHSTO.N'E  results  are  from  a  quasi-static  analysis, 
which  are  not  as  accurate  at  high  frequencies.  The  full- 
wave  results  and  the  TOUCHSTONE  results  for  the  phase 
of  533  arc  also  compared  and  arc  shown  in  Fig.  4.  Good 
agreement  is  found  below  10  GHz,  but  more  than  45*  dis¬ 
crepancy  is  found  at  24  GHz  in  this  particular  case. 


Fix.  3  Mtgnitudc  of  5  pwiictcn  of »  biuic  T-junction 
(i»  «  9  9 ,  A  •  3$  mil.  «  34  mtl) 


Fic.  4  Phane  o(  533  of  a  bwk  T-junction 
•  9  9 ,  A  «  3S  nul ,  u*  ■  34  mil) 


Fig.  5  Magnitude  of  5  parameten  of  a  right 
angle  croaaing  junction 
{u  «  10.3  ,  A  H  25  mi) ,  w  ■  34  mil) 

N«U  IkM  ia  TOUCHSTONE  |5n| « {.s'lli 

A  symmetrical  cross  junction  is  shown  in  Fig.  2(b), 
where  two  identical  transmission  lines  are  crossed  at  a  right 
angle.  With  the  same  numerical  convergence  criteria  as  in  a 
basic  T-junction,  the  magnitudes  of  scattering  coefficients 
are  shown  in  Fig.  5.  TOUCHSTONE  predicts  equal  power 
distribution  for  a  cross  junction  in  a  wide  frequency  range; 
however,  the  results  of  this  study  indicate  there  is  an  un¬ 
equal  power  distribution  for  the  cross  junction.  The  phe¬ 
nomenon  of  unequal  power  distribution  is  more  significant 
for  higher  frequencies.  From  a  distributed  circuit  point  of 
view,  this  phenomenon  is  obvious  since  port  11  and  port 
111  are  not  symmetric  and  the  current  tends  to  go  straight 
through  the  cross  junction.  For  a  quasi-static  calculation, 
the  cross  junction  is  like  two  wires  jointing  together  and  in 
terms  of  the  lumped  circuit  concept,  the  power  distributed 
in  each  port  is  certainly  identical. 

The  problem  for  a  basic  junction  is  that  the  power 
transmission  distributed  in  each  port  is  usually  restricted. 
For  example,  in  a  basic  T-junction,  there  is  always  more 
power  transmitted  in  port  11  than  in  port  111,  and  the  re¬ 
flection  coefficient  increases  as  the  frequency  increases.  In 
order  to  make  the  circuit  design  more  flexible,  modified  or 
compensated  discontinuities  are  usually  used  [9],  A  shaped 
T-junction  shown  in  Fig.  2(c)  is  an  example  of  this  mod¬ 
ification.  In  the  present  full-wave  analysis,  the  advantage 
of  using  piecewise  sinusoidal  basis  functions  can  be  seen  in 
this  particular  application.  By  using  piecewise  sinusoidal 
basis  functions  inside  the  module,  the  shape  of  the  junc¬ 
tion  can  be  quite  flexible  in  the  modeling.  The  design  of  a 
shaped  T-junction  shown  in  Fig.  2(c)  is  intended  for  equal- 
power  transmission  in  ports  11  and  HI  and  small  reflection 
coefficient  at  port  HI.  If  the  difference  of  the  transmitted 
power  for  ports  II  and  III  is  required  to  be  less  than  0.5  dD, 
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Fid.  6  Madnilude  of  5  pArurncten  of  o  ihoped  T<junction 
(<.  ■  10.2  .  A  ■  25  nil) ,  w  ■  24  nil) 

the  results  in  Fig.  6  show  that  equal  power  transmission  is 
very  broad  band  (from  D.C.  up  to  16  GHz).  In  compari¬ 
son,  for  a  basic  T-junction,  equal  power  transmission  only 
valid  up  to  6  GHz.  Besides,  533  of  a  shaped  T-Junction  is 
usually  a  few  dU  lower  than  that  of  a  basic  T-junction. 

IV.  Conclusion; 

A  full  wave  analysis  of  a  multi-port  network  is  car¬ 
ried  out  by  the  moment  method.  It  has  the  capability  of 
analyzing  four-port  irregular  shaped  junctions.  The  results 
for  a  tee  and  a  cross  junctions  are  presented  and  found  in 
good  agreement  with  the  quasi-static  results  at  low  frequen¬ 
cies,  An  example  of  using  a  shaped  T-junction  to  improve 
the  performance  of  a  basic  T-junction  was  given.  This  ex¬ 
ample  also  illustrated  the  flexibility  and  the  CAD  potential 
of  the  full-wave  analysis  presented  in  this  paper.  The  com¬ 
puted  results  arc  further  compared  against  the  measured 
data  for  a  tee  junction.  The  comparison  shows  excellent 
agreement. 
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Abstract 

A  generalized  method  for  calculating  both  ra¬ 
diation  and  surface-wave  losses  is  developed  for  mi¬ 
crostrip  discontinuities.  The  losses  are  determined 
by  a  rigorous  Poynting  vector  analysis  where  the 
current  distribution  over  the  entire  microstrip  dis¬ 
continuities  is  a  result  of  a  full-wave  moment  method 
solution.  It  is  found  that  above  a  certain  frequency, 
the  surface-wave  loss  becomes  more  important  than 
the  radiation  loss.  A  self-consistency  check  of  the  re¬ 
sults  based  on  power  conservation  is  also  presented. 

Introduction 

Radiation  and  surface  waves  are  unavoidable 
physical  effects  of  microstrip  discontinuities  associ¬ 
ated  with  an  open  structure.  In  recent  years,  a 
full-wave  analysis  that  includes  these  physical  effects 
has  been  developed  for  various  microstrip  disconti¬ 
nuities  [1-5],  Although  both  radiation  and  surface 
waves  are  included,  the  analysis  only  provides  the 
total  losses.  One  still  can  not  distinguish  the  per¬ 
centage  of  power  losses  due  to  radiation  and  due 
to  surface  waves.  The  full-wave  analysis  empioys 
the  moment  method  to  find  the  currents  in  the  ,i. 
crostrip  circuits  and  subsequentiy,  the  circuit  pa¬ 
rameters  of  the  discontinuities.  From  an  antenna 

i  This  research  was  supported  under  U.S. 


point  of  view,  once  the  currents  in  a  conductor  are 
known,  the  time-harmonic  fieids  can  be  computed. 
EVom  the  theory  of  printed  circuit  antennas  [6],  the 
radiated  space  waves  are  sphericai  waves  in  the  hemi¬ 
sphere  above  the  substrate;  whiie  surface  waves  are 
cylindrical  waves  guided  along  the  planar  direction 
of  the  substrate  and  decay  exponentially  toward  the 
free  space.  The  power  due  to  radiation  and  sur¬ 
face  waves  can  therefore  be  computed  separately 
through  a  rigorous  Poynting  vector  analysis.  In  this 
work,  radiation  and  surface-wave  losses  for  several 
types  of  micrpstrip  discontinuities  are  investigated. 
These  discontinuities  include  open-end,  right-angle 
bend,  rectangular  patch  and  overlay  electromagnet- 
ically  coupled  (EMC)  lines. 

Analysis 

In  the  integral  equation  formulations,  ail  field 
components  can  be  expressed  in  terms  of  the  dyadic 
Green’s  functions  and  the  current  components.  There¬ 
fore,  the  current  distribution  over  an  entire  discon¬ 
tinuity  can  be  treated  as  a  basic  block  (7)  for  this 
analysis.  The  radiation  loss  can  be  calculated  by  in¬ 
tegrating  the  Poynting  vector  over  an  infinite  plane 
(shown  in  Fig.  1)  in  the  free  space  above  the  entire 
microstrip  circuits.  The  expression  is  written  as 
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(1) 


Fi(.  1.  Integration  plane  for  calculating  radiation  loss. 


Kig.  2.  Inlcgralion  cylinder  fur  calculating  surfacc'wave  loss. 
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Surface  wavea  can  be  obtained  from  the  residues 
of  the  Fourier  integrals  in  the  spectral  domain  ap¬ 
proach.  With  the  characteristics  that  surface  waves 
propagate  along  the  surface,  the  surface-wave  loss 
can  be  found  by  integrating  the  Poynting  vector  over 
a  cylinder  (shown  in  Fig.  2)  of  large  radius  p.  The 
expression  of  the  surface-wave  power  is 

Pm  ~  ^  ii't,)  •  ppd^z  (2) 

Numerical  Results  and  Discussions 

Due  to  power  conservation,  the  summation 
of  each  power  should  be  equivalent  to  the  incident 
power.  An  example  of  this  check  is  shown  in  Table  I. 
The  incident  power  is  normalized  to  1.  Due  to  power 
conservation,  total  power  which  is  the  summation 
of  reflected,  transmitted,  radiated  and  surface-wave 
power  should  be  equal  to  1.  The  result  in  Table  I 
shows  excellent  agreement. 

Figs.  5-5  show  the  percentage  of  radiation 
loss,  surface-wave  loss,  and  total  power  loss  as  a 
function  of  frequency  for  open-end,  right-angle  bend, 
and  overlay  electromagnetically  coupled  lines  respec¬ 
tively.  It  is  seen  that  the  loues  due  to  both  radia¬ 
tion  and  surface  waves  increase  with  frequency.  At 
low  frequencies,  the  losses  are  mainly  due  to  radi¬ 
ation.  When  the  frequency  increases,  surface-wave 
loss  increases  faster  than  the  radiation  loss.  Above 
a  certain  frequency,  the  surface-wave  loss  is  more 
signiflcant  than  the  radiation  loss.  Fig.  0  shows  the 
power  distributions  of  a  rectangular  patch  antenna. 
It  can  be  seen  that  the  maximum  radiation  efficiency 
is  about  05%  in  a  very  narrow  band  around  7.2GHz 
with  about  10%  surface  wave  loss  and  25%  return 
loss. 


Table  I.  Ko«er  conservaiion  check  for  riglu-angle  bend  disconiinuiiy. 
I  Parameieri  uc  (he  same  as  (hose  in  Fig.  4. ) 
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Fig.-O.  Power  diitributioiis  versus  frcrjueiicy  for  a 
r«ct«u|ular  patch  aiiteliiia.  (c^«3.33,  wl  = 
SOiiiil.  Iis63uiil.  w2sS00iiii|.  IsSOOiiiil) 


A  generalized  method  to  distinguish  power 
losses  due  to  radiation  and  surface  waveris  presented. 
This  method  should  aid  in  CAD  for  minimizing  power 
losses  launched  into  radiated  space  waves  and  sur> 
face  waves  from  arbitrary  aiicrostrip  discontinuities 
and  maximizing  the  radiation  efficiency  for  arbitrary 
patch  antennas.  Power  loss  mechanum  in  other 
types  of  microstrip  discontinuities,  such  as  gap,  step 
and  stubline  will  be  easily  implemented  in  this  anal¬ 
ysis. 
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This  thesis  is  mainly  composed  of  two  parts.  The  first  part  is  to  deter¬ 
mine  the  frequency-dependent  characteristics  of  microstrip  lines  in  a  three  layer 
shielded  structiue.  Numerical  results  given  in  this  part  include  the  propagation 
constant,  characteristic  impedance  and  transverse  dependence  of  longitudinal  and 
transverse  currents  of  both  microstrip  single  and  coupled  lines  for  a  variety  of  di¬ 
electric  configurations.  The  second  part  is  to  analyze  the  behavior  of  several 
types  of  microstrip  line-to-line  transitions.  The  methodology  employs  a  mo¬ 
ment  method  procedure  where  the  combination  of  subdomain  expansion  modes 
and  entire  domain  microstrip  modes  is  used.  The  transitions  studied  include 
edge-coupled  lines,  overlay-coupled  lines  and  coupled-to-single  lines.  A  power 
conservation  check  based  on  a  rigorous  Poynting  vector  analysis  is  also  used  to 
determine  the  accuracy  of  the  numerical  convergence. 


Chapter  I 


Introduction 

Shielded  microstrips  are  commonly  encountered  in  microwave  and  mil¬ 
limeter  integrated  circuits.  The  advantages  include  relatively  low  losses  and 
fairly  high  Q-factor  as  compared  with  open  microstrips[l].  Many  applications 
such  as  high  directivity  directional  couplers  and  broadband  filters  have  been  re¬ 
alized  in  shielded  structures[2][3].  In  antenna  applications,  shielded  microstrips 
can  provide  nice  features  when  used  to  feed  a  linear  array  of  slots  in  the  outer 
conductor  wali[4].  For  sufficiently  low  frequencies,  the  quasi-static  theory  can 
be  employed  to  characterize  shielded  microstrip  lines  and  their  discontinuities 
reasonably  well[5]-[8].  However,  as  the  frequency  increases,  the  deviation  from 
quasi-static  behavior  becomes  significant  and  a  more  rigorous  full-wave  analysis 
is  necessary. 

To  calculate  accurately  the  dispersion  characteristics  of  shield^  microstrip 
lines,  the  spectral-domain  analysis  proposed  by  Itoh  and  Mittra[9]  is  numerically 
simpler  and  niore  efficient  than  the  conventional  space-domain  techniques.  This 
method  allows  one  to  convert  convolutions  into  algebraic  equations  in  the  Green’s 
function  formulation;  thus  avoiding  the  necessity  of  the  evaluation  of  complicated 
integrals.  Since  the  strip-width  and  strip-offset  are  comparable  to  the  waveguide 
dimensions,  the  transverse  current  component  should  not  be  neglected  and  a 
complete  dyadic  Green’s  function  is  required.  After  deriving  the  spectral-domain 
Green’s  function  for  a  three  layer  shielded  structure,  boundary  conditions  on  the 
strips  are  employed  to  obtain  integral  equations.  As  a  result,  the  propagation 
constant  and  current  distributions  on  the  strips  can  be  determined  by  the  method 
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of  moments.  An  impedance  definition  using  power  and  total  currents  are  chosen 
to  d^cribe  the  exchange  of  power  with  TEM  lines  [10].  Once  the  current  dis¬ 
tributions  on  the  strips  are  found,  the  six  field  components  can  be  obtained  and 
the  characteristic  impedance  can  be  calculated  through  performing  the  surface 
int^al  of  the  Poynting  vector  in  the  propagating  direction. 

Based  on  the  analysis  described  in  Chapter  II,  the  computer  program 
developed  for  calculating  the  dispersion  characteristics  of  shielded  microstrip  lines 
can  offer  the  capacity  to  analyze  the  shielded  microstrip  single  and  coupled  lines 
with  arbitrary  strip-width  and  strip-offset,  while  only  limited  results  in  some 
special  cases  were  reported  by  a  few  authors  [12]-[15]. 

Several  types  of  shielded  microstrip  line-to-line  transitions  are  also  studied 
in  this  thesis  as  the  important  building  blocks  for  high  frequency  inter-connects. 
Applications  in  nun-wave  integrated  circuits  include  high-pass  filter,  multiplex¬ 
ers  and  directional  couplers.  Losch  [17]  has  designed  a  broadband  highpass  filter 
in  realization  of  an  overlay-coupled  line  transition  based  on  a  quasi-static  for¬ 
mulation.  A  more  rigorous  full-wave  analysis  for  coupled  line  filters  associated 
with  the  open  structure  has  been  discussed  by  Katehi  [18]  for  an  edge-coupled 
transition  and  by  Yang  and  Alexopoulos  [20]  for  an  overlay-coupled  transition. 
In  [20]  a  spectral-domain  approach  by  expanding  the  current  in  the  coupled  line 
section  with  a  combination  of  entire  domain  and  subdomain  modes  is  used.  This 
mode  expansion  mechanism  seems  to  be  the  most  efficient  and  fruitful  by  far.  For 
the  advantage  of  preventing  unnecessary  interaction  and  radiation  loss,  a  waveg¬ 
uide  housing  is  sometimes  more  practical  to  the  real  circuit  design.  In  Chapter 
III,  a  full-wave  moment  method  is  used  to  characterize  shielded  line-to-line  tran¬ 
sitions.  The  transitions  including  edge-coupled  lines,  overlay-coupled  lines  and 
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coupled-to-single  lines  are  investigated. 


For  shielded  microstrip  transitions,  the  incident  power  should  be  equal  to 
the  summation  of  reflected  power,  transmitted  power  and  some  losses  coupled 
to  the  higher  order  modes  and  the  multi-layered  waveguide  modes.  With  proper 
waveguide  dimensions,  the  losses  coupled  to  the  higher  order  modes  and  the 
multi-layered  waveguide  modes  can  be  removed.  Therefore,  a  power  conservation 
among  the  incident,  reflected  and  transmitted  powers  can  be  checked  to  determine 
the  numerical  accuracy.  In  Chapter  III,  the  scattering  coefficients  are  obtained 
through  the  method  of  moments.  However,  due  to  the  different  characteristic 
impedance  between  the  feedline  and  the  parasitic  coupled  line,  these  coefficients 
can’t  contain  the  information  of  power  conservation.  With  the  characteristic 
impedance  obtained  through  the  Poynting  vector  analysis  described  in  Chapter 
II,,  the  reflected  power  associated  with  the  feedline  and  the  transmitted  power 
associated  with  the  parasitic  coupled  line  can  be  normalized  respectively  by  the 
product  of  the  characteristic  impedance  and  the  square  of  the  absolute  value  of 
scattering  coefficient  for  each  line. 

The  results  of  power  distributions  and  scattering  coefficients  of  line-to-line 
transitions  can  identify  the  properties  and  applications  of  each  transition. 
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Chapter  II 


The  Dispersion  Characteristics  of  Three  Layer 
Shielded  Microstrip  Lines 

2.1  Greenes  Function  for  Three  Layer  Shielded  Microstrips 

In  this  section,  the  dyadic  Green’s  function  for  an  electric  Hertzian  dipole 
embedded  in  a  three  layer  shielded  structure  is  derived  in  a  spectral  domain. 
The  geometry  shown  in  Fig.  2.1  contains  three  different  dielectric  layers  w'ith 
permittivity  £i,«2,£3,  and  thickness  respectively.  The  permeability  in 

each  region  is  assumed  fio,  the  free  space  permeability.  The  total  thickness  is  b 
and  the  distance  between  two  sidewalls  is  a.  The  dipole,  oriented  in  the  planar(x 
or  z)  direction,  is  located  either  at  y  =  hj  or  at  y  hj  +  hj. 


Fig.  2.1  An  electric  Hertzian  dipole  in  a  three  layer  shielded  structure. 
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and  at  y  =  Aj, 


*|n,.  =  *'n,„ 

(2.14) 

dy 

(2.15) 

callfy  =  Cilliy, 

(2.16) 

5y 

^nfy  an^M  ana, 

5y  5x  5x  ’ 

(2.17) 

where  11, •  =  n,*x  + 11, yy  +  n,-,i,  n,y  =  Ilfy  +  Hfy  and  fc,-  is  equal  to  Uy/lim. 

By  matching  the  boundary  conditions  that  the  tangential  electric  fields 
are  zero  at  x  =  0,  x  =  a,  y  =  0  and  y  =  6,  the  two-dimensional  Fourier  transform 
of  the  Hertz  potential  in  each  region  can  be  defined  as, 


in  region  3,  . 

n3*(®,y,«)  =  ^  J]  J  >l*8inh  93(6 -y)  sin 

||S«*00 

n3,(«,y,«)  =  r--  J]  /  A,8inh93(6-y)co8Q(„xe-’»'**d^, 

nya;,y,z)  =: -i-  J]  [  A',  cosh  93(6 -  y)  sin  OnXc"^'**d)S, 

2S'0  n.-oo-'-* 

n3»(®,  y,  ^  53  j  ^  ?3(^  -  y)  sin  o„xe"''^*d^, 


(2.18) 

(2.19) 

(2.20) 
(2.21) 


in  region  2, 

n3,(x,  y,  z)  =  53  /  (B,8inh92y  +  C,co8h9ay)sinonXc"-''’*di3,  (2.22) 

‘"■a  n— oo''-«> 

1  ^00 

nj*(x,  y,  z)  s=  5 —  53  /  (^»  sinh  93y  +  Cm  cosh  9ay)  cos  (2.23) 

„,_oo  •'-oo 

n3y(x,y,z) » -i-  2  /  (^i»inh92y4'C'ico8h9ay)8inofnxe“^'’*<//9,  (2.24) 
njy(x,y,z)  =  -J-  f;  r(B;sinh9ay  +  C;co8h9ay)8ino„xe--»^*d/?,  (2.25) 
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and  in  region  1, 

1  O®  ^OO 

niz(x,y,z)  =  T —  53  /  ■D*sinh9iysina„xe"'''**d)3, 

®  n=— oo  ’'~°® 

1  ®®  1*00 

ni*(x,y,x)  =  - —  V]  /  Dj.sinh^iycosanxc'-'^V^, 


(2.26) 

(2.27) 


L  yoo 

nyx,y,x)  =  5 —  53  /  cosh  9iy  sin  (2.28) 

nss-oo  •'-«> 

1  ®®  TOO 

nfy(x,y,x)  =  5 —  5]  /  i?*  cosh  9iy  sin  a„xc~''^*d/3,  (2.29) 

^_oo  •'-00 

where  g,-  =  and  an  =  ~.  The  coefficients  A,,  B*,  B* ,  •  •  • , 

are  determined  by  substituting  Eqs.  [2.18]-[2.29]  into  Eqs.  [2.2]-[2.17].  The  results 
are  listed  in  Appendix  A. 

The  lateral  electromagnetic  fields  in  region  i  are  related  to  the  Hertz  po¬ 
tential  through 


F<  —  fc^TT  4.  ^  I  I  \ 

B«-*,n„+^^(  +  Qy  +  Qg  )l 

F-  —  Pn-  4.  ^  I  5n,tv 

F  —  t^TT  4.  ^  4.  4. 

E»  -  k,  n,-,  +  _(._  +  -^  +  -—). 


.  vaait 

dy 

dz  ' 

e 

an.., 

dz 

dx  ‘ 

dUiy 

an.., 

(2.30) 

(2.31) 

(2.32) 

(2.33) 
(2.3-1) 

(2.35) 


Taking  the  two-dimensional  Fourier  transform  of  the  electric  field  compo¬ 
nents  in  Eqs.  [2.30]-[2.32],  the  spectral-domain  dyadic  Green’s  function  at  each 
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interface  is  then  obtained  as  follows 


Eix{y  =  hi)  =  Giux cos QnX  COS +  Guxz cos a„x  sinoffiXi,  (2.36) 
Eiz{y  =  hi)  =  (5ii« sin OnX sin a„X|  +  sin a„x cos Q(„Xj,  (2.37) 
^uiy  =  hi  +  h2)  =  G2i*x  cos  a„x  cos  a„Xi  +  G2irzCosa„xsinanaJi>  (2.38) 
Eiz{y  =  hi  +  h2)  =  G2i*,sinanxsinanxi  +  G2izxSina„xcosan®l-  (2.39) 

For  a  dipole  oriented  in  the  planar  direction  and  located  at  x  =  Xj,  y  = 
hi  +  h2  and  x  =  0,  the  expressions  for  the  dyadic  Green’s  function  are  the  same 
except  replacing  Ci,hi,xi  by  csiha.Xj,  respectively.  The  expressions  are 

E2x{y  =  hi)  =  Cli2xx  cos  a„X  cos  OnXj  +  Gi2*,  cos  OnX  sin ttnXj,  (2.40) 
E2»{y  =  hi)  —  Gi2i<sinanXsinanX3  +  Gmx sin anX cos anXj,  (2.41) 
&2x{y  =  hi  +  h2)  =  CI22XX  cos  OnX  cos  QnX2  +  ^23**  COS  0«X  sin  Q!n»2*  (2.42) 
^3«(y  »  hi  +  ha)  =  (^33„  sin  OnX  sin  anX2  +  Otux  sin  a^x  cos  OnX^.  (2.43) 

In  Eqs.  [2.36]-[2.43]  the  functions  are  defined  as 

£  r£4(a.,Ax,»')e-J«— 'MyJ  (2.44) 

««-oo 

and  the  spectrahdomain  Green’s  functions  found  in  Appendix  B. 

2.2  Integral  Equations  for  Calculating  Propagation  Constants 

By  superposition,  the  total  electric  fields  due  to  the  currents  on  the  mi¬ 
crostrips  are 

E(x,  y,  r)  »  jj  ^"(x,  y,  x;  x,  x')  •  J{x\  z)dxdz\  (2.45) 
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where  the  spatial  domain  dyadic  Creen’s  function  is  equivalent  to 


in  1  ^  fooq?  fsmanXi  sinttni  I  \ 

T^{x,y,z;x,z)  =  —^  I  CI{anJ,y)\  (2.46) 

°  n=-ocf  “<»  ^cos  o„x y  ycos  o„x  y 

The  spectral-domain  dyadic  Green’s  functions  G  at  each  interface  are  expressed 
in  Appendix  B.  For  microstrip  lines  located  at  y  =  hi  and  hi  -f  ^2,  the  currents 
on  strips  are  traveling  along  the  z-axis  with  propagation  constants  and  can  be 


written  as 


Jiix,  y,  z)  =  {Jitx{x)x  +  Juzix)z)e~ 


(2.47) 


Ji{x,  y\  z)  =  ( J2te(x')*  +  J2txix)z)e~^^'*  , 


(2.48) 


respectively.  With  the  boundary  conditions  that  tangential  electric  fields  are  zero 
on  the  strips,  the  integral  equations  can  be  obtained  by  substituting  Eqs.  [2.47], 
[2.48]  into  Eqs.  [2.45],[2.46].  For  each  (Xf  z)  on  the  microstrips,  these  integral 
equations  are 

•  «  P  * 

^u(»,y  *  hi,*)  /  Jiu(**)»inon*  d** 

£^ir(x,y  =  hi,z)  1  -is.,  ^  T/  13  N  fJitx(x)co8anx'dx' 

“  -e  L  Oi(Qnf^p,x) 

E2z{x,y  =  hx  +  h2,*)  ®  »»-oo  f  Jng(x  Js’inOnX  dx 

E2x(x,  y  =  hi  +  ha,  z)  f  datx(®  )  cos  Onxdx’ 


(2.49) 
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where 


Cjn«sma„x  Gn,*smanX  Gi2«sina„x  Gi2*xSinanX 
(5iix*COSa„X  Gii„C08anX  Gi2x*C0Sa„X  Gi2xxCOSanX 

G2i«sina„x  G2i,xsina„x  G22*»sina„x  G222xsina„x 
G2ix*C08a„X  G21*xC0SanX  G22xxCOSa(„X  G22x»COSa„X 


(2.50) 


The  method  of  moments  can  be  used  here  to  solve  the  integral  equations 
accurately  for  propagation  constants  and  simultaneously  obtain  a  good  ap¬ 
proximation  of  the  actual  current  distributions.  A  proper  choice[16]  of  expansion 
functions  to  satisfy  the  asymptotic  requirement  of  the  current  distribution  near 
the  edge  is 

,  at, 

•At.(ar<)  =  2  m  ss  1, 2,  •  •  • ,  (2.51) 

mat 

and 

,  N, 

•A«x(®j)  *  2)  •  •  • ,  iVx,  (2.52) 

mal 


where 


(2.53) 


(2.54) 


and  x)  =  x*-^,-.  Si  and  ib<  are  strip-offset  and  strip- width,  respectively,  associated 
with  microstrip  i.  The  functions  Tm  and  Un  are  Chebyshev  polynomials  of  the 
first  and  second  kinds,  respectively.  These  expansion  functions  [2.53],[2.54]  are 
shown  in  Fig.  2.2.  Following  the  Galerkin  procedure  applied  in  the  spectral  do¬ 
main,  electric  fieid(s  in  Eq.  [2.49]  are  tested  with  the  expansion  functions 
and  This  results  in  a  set  of  linear  equations.  These  linear  equations,  when 
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expressed  in  matrix  form,  are 


where 

[ZiUx\  [^llxs]  [Zuxx]  [2i2*x] 

[^nxx]  [■^11**1  [Zl2xx]  l^lTxx] 

[ZiXxx]  [Z2lxx\  [^22*»]  [^22*»] 

[^21**1  [Z2\xx]  [^22x*]  [^22xx] 

The  matrix  elements  in  each  submatrix  are  in  the  general  form  of 

z"*4(«=  E 

n=-oo 

where  J  is  defined  as 

Jt’;(a:;)sinan(x;  +  Si)dx\, 

•'“T 

^Tx  =  / 1  “«(*»•  + 

With  the  Bessel  function  identities 


7(«  = 


=  0, 


and 


e^'^*7;(x) 
\/l'  — 


dx 


Z  </-l 

Eqs.  [2.58], [2.59]  can  be  further  expressed  in  closed  forms,  which  are 


and 


+ 


cos  a„5,  /m 

fOnWi 

J  ‘'m-Xl 

^  2  ■ 

sin  On^i  Rx 

■m-1 J  / 

OlnWi. 

J  •/m-li 

2  ‘ 

c«a,<i  R,  [r-’mJ™(2^)^] 


(2.55) 

(2.56) 

(2.57) 

(2.58) 

(2.59) 

(2.60) 

(2.61) 

(2.62) 

(2.63) 
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In  solving  Eq.  [2.55]  to  obtain  the  coefficients  Am  and  Bm,  the  determinant 
should  be  set  to  zero  for  non-trivial  solutions.  Therefore,  the  matrix  elements 
[2.57]  can  be  evaluated  numerically  for  any  ^  and  the  values  which  force  the 
determinant  of  ^  to  zero  are  the  propagation  constants.  In  other  words,  the 
propagation  constant  must  satisfy  the  characteristic  equation 

det  l{0p)  =  0.  (2.64) 


2.3  Spectral  Domain  Calculation  of  Characteristic  Impedance 

Once  the  current  distributions  are  obtained,  the  characteristic  impedance 
of  three  layer  shielded  microstrip  lines  can  be  calculated  using  the  ratio  of  average 
power  to  half  of  the  square  of  the  total  longitudinal  current[10].  The  formula  of 
the  characteristic  impedance,  for  a  single  line,  is 


(2.65) 


while  for  symmetric  coupled  lines,  is 


,  rj  ~  p  '  (2.66) 

The  average  power  is  the  surface  integral  of  the  z  component  of  the  Poynting 


vector 


(2.67) 


It  is  noted  that  the  axial  currents  of  asymmetric  coupled  lines  /,i  and  /jj  are  not 
the  same.  Therefore,  in  order  that  the  definition  of  characteristic  impedance  be 
useful,  the  formula  described  in  [11]  can  be  used: 


„  2Pt  FU  /o*  X  H*pi)  •  zdydx 
Lpi  —  -=5—  ss  - . 


(2.68) 
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where  and  /*p<  are  respectively  the  average  partial  power  and  axial  current 
associated  with  the  microstrip  i  for  the  mode  p.  In  the  above,  Ept  is  the  total  elec¬ 
tric  field  and  Hpi  is  the  magnetic  field  generated  by  the  current  on  the  microstrip 
i.  With  similar  approach  in  deriving  Eq.  [2.49],  the  transverse  components  of 
electromagnetic  fields  in  region  i  can  be  expressed  as 

/  Jux{x)  sin  anx'dx 

I  Jttx{x')co3anxdx' 

S  sin  anx'dx' 

S  J7Ue{x)  cos  an^t'dx' 
f  Jux{x)  sin  a^xdx 
f  J\tx{x')  cos  QnXdx’ 
f  Jitxix)  sin  OnX  dx 
f  J2tx{x)  cos  Onx'dx 

f  Jux(^')9inanxdx' 
f  Jifx(x)  cos  a„xdx 
/  d3tx(x')  sin  OnXdx 
f  J2tx{x)cOS  OnXdx 

I  Juxix)sinanxdx 

f  Jltx{x)  cos  Qnxdx 
J  Jitxix)  sin  Onxdx 
J  Jiuix)  cos  anx'dx 

where 


E,  = 

[  Eiigg 

^i\x* 

Ei2xx 

^i2xx  ]  * 

(2.73) 

’Ey  = 

[  ^ilyx 

Eilyx 

Ei2yx 

iiityx  ]  > 

(2.74) 
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(2.75) 


^  S[ilxx  ^ilxx  ^i2xz  Hi2xx  j  > 


^  Ailyx  Hilyg  Hnyx  Hflyx  j  *  (2.76) 

Substituting  Eqs.  [2.69]-[2.72]  into  E)q.  [2.67],  one  can  proceed  the  integration 
with  respect  to  x  to  yield 


/  <J\tt{x)  smonxdx 

D  1  f  J  I  Jiix{^')coBa^xdx 

n=-oo  f  J2u(x)  sin  Qnx'dx 

^  f  J7tx(x')  cos  Qnx'dx'  ^ 

I  dut{x)  sin  anxdx 

t  S  Jitx{x)co8anxdx 

X  L  ^  .  (2-77) 

*•  Hy  /Jjt,(x)8inan®<^® 

^  J  J2txW)  cos  OnXdx  ^ 

where  Pi  is  the  average  power  in  region  *.  The  integration  over  y  can  be  also 
solved  in  closed  form  and  the  results  are  shown  in  Appendix  C. 


(2.77) 


The  average  partial  power  can  be  calculated  in  a  similar  way  by  con¬ 
sidering  only  Jux  and  Ju,  in  Eqs.  [2.71],[2.72j. 


2.4  Numerical  Results  and  Discussions 


The  shielded  microstrip  lines  considered  in  this  analysis  are  shown  in 
Figs.  2.3(a)-(d)  where  infinitesinudly  thin  strips  and  the  ground  plane  are  as¬ 
sumed  to  be  perfect  conductors.  It  is  also  assumed  that  the  substrate  material  is 
lossless.  The  propagation  constant  associated  with  these  structures  can  be  deter- 


Fig.  2.3(a)  Two  layer  single  line 


h 

K—  a  — H” 

Fig.  2.3(c)  Symmetric  coupled  lines 


Fig.  2.3 


Fig.  2.3(b)  Three  layer  single  line 
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Fig.  2.3(d)  Asymmetric  coupled  lines 
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mined  by  solving  Eq.  [2.64]  and  the  approximated  current  distributions  can  be 
obtained  simultaneously.  The  longitudinal  current  distribution  is  normalized  such 
that  it’s  integration  over  the  entire  strip-width  is  equal  to  unity.  The  transverse 
current  distribution  is  also  normalized  in  relative  magnitude  compared  to  the 
longitudinal  current  distribution.  Once  both  current  distributions  are  known,  one 
can  solve  the  six  field  components,  and  therefore  the  characteristic  impedance  can 
be  computed  through  a  power-current  definition.  The  dispersion  characteristics 
for  both  single  and  coupled  lines  are  evaluated  numerically  and  their  dependence 
on  dielectric  constant,  strip-width,  strip-offset  and  frequency  will  be  discussed  in 
this  section. 

Dispersion  characteristics  for  the  single-line  configuration  are  shown  in 
Figs.  2.4-2.7.  In  Fig.  2.4,  the  dispersion  characteristics  are  computed  for  various 
dielectric  constants.  It  is  seen  that  the  increase  in  permittivity  increases  the  prop¬ 
agation  constant  and  decreases  the  characteristic  impedance,  regardless  of  which 
layer  the  strip  is  located.  Fig.  2.5  shows  the  variation  of  dispersion  characteris¬ 
tics  with  the  strip-width.  As  the  width  increases,  the  transverse  current  density 
also  increases  and  the  longitudinal  current  distribution  becomes  more  oscilla¬ 
tory.  Both  propagation  constant  and  characteristic  impedance  drop  significantly 
as  the  width  is  close  to  the  distance  between  two  sidewalls.  Fig.  2.6  shows  the 
dispersion  characteristics  as  a  function  of  the  strip-offset.  The  transverse  current 
distribution  increases  asymmetrically  with  the  increase  of  strip-offset  6.  Both 
propagation  constant  and  characteristic  impedance  decrease  with  the  increase  of 
strip-offset  6.  Fig.  2.7  shows  the  frequency  behavior  of  a  single  line.  The  be¬ 
havior  of  the  current  distributions  is  similar  to  the  dependence  on  strip-width 
because  the  strip-width  becomes  electrically  wider  as  the  frequency  increases. 
Both  propagation  constant  and  characteristic  impedance  increase  with  frequency 
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duracttrisiic  Impedaiice  Zc  Normalized  Propagation  Cbnstant  p 


4J 


easel — — O' 


1 — n - 1 - 1 — - 1 - r 


1234S6789  10 

Permittivity  e 
(a) 


(b) 

Fig.  2.4  Characteristics  of  a  single  line  for  the  configuration  of  Fig.  2.3(b). 

(a)  Normalized  propagation  constant  versus  permittivity  e. 

(b)  Characteristic  impedance  versus  permittivity  e. 

Case  1:  cj  —  2.2,  ci  =  €3  =*  c,  Case  2:  cj  =  10.2,  ci  =  C3  =  c. 

Case  3:  ci  =  £3  =  2.2,  cj  —  e.  Case  4:  ci  =s  £3  =s  10.2,  £3  =  c, 
hi  ss  ha  =s  *3  =s  O.O2A0,  =  1,  ^  -  10, 5  =  0  for  all  cases. 
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Fig.  2.5  Characteristics  of  a  single  line  for  the  configuration  of  Fig.  2.3(a). 
(a)  Normalized  propagation  constant  and  characteristic  impedance  versus 
strip-width,  (b)  Normalized  longitudinal  current  density  versus  strip  posi¬ 
tion.  (c)  Relative  transverse  current  density  versus  strip  position. 

Cr  =  10.2,  h  =  0;02A<„  f  =  0.0,  ^  =  2.0,  ^  =  10.0 


8-0.00 


Fig.  2.6  Characteristics  of  a  single  line  for  the  configuration  of  Fig.  2.3(a). 
(a)  Normalized  propagation  constant  and  characteristic  impedance  versus 
strip-offset,  (b)  Normalized  longitudinal  current  density  versus  strip  posi¬ 
tion.  (c)  Relative  transverse  current  density  versus  strip  position, 
e,  *  10.2,  ft  s  0.02Ao,  f  *  1.0,  «  2.0,  ^  *  10.0 


Nomutized  nropagation  Constant  P 


h/Xo  h/Xo 

(a)  (b) 


Fig.  2.7  Characteristics  of  a  single  line  for  the  configuration  of  Fig.  2.3(a). 

(a)  Normalized  propagation  constant  versus  normalized  frequency. 

(b)  Characteristic  impedance  versus  normalized  frequency. 

(c)  Normalized  longitudinal  current  density  versus  strip  position. 

(d)  Relative  transverse  current  density  versus  strip  position. 

- 10.2,  f  -  1.0,  jt »  0.0,  V  =  2.0,  f  =  10.0 
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in  this  case.  The  number  of  expansion  functions  used  to  calculate  propagation 
constant  is  much  less  sensitive  than  that  used  to  calculate  the  characteristic 
impedance. 

Figs.  2.8-2.12  show  the  dispersion  characteristics  for  coupled  lines.  In 
Fig.  2.8,  two  dominant  modes,  even-mode  and  odd-mode  of  symmetric  coupled 
lines  are  shown  as  a  function  of  the  strip-offset  8  and  are  compared  with  the  results 
of  a  single  line  in  some  special  cases.  The  even-mode  for  ^  =  0  corresponds  to 
the  dominant  mode  of  a  single  line  with  the  strip-width  twice  wide,  while  the 
odd-mode  corresponds  to  the  dominant  mode  of  a  single  line  with  an  electric  wall 
at  z  =  The  comparison  shows  good  agreement.  Fig.  2.9  shows  the  dispersion 
characteristics  of  asymmetric  coupled  lines  for  two  independent  modes.  Mode 
1,  which  is  similar  to  the  odd  mode  of  the  symmetric  coupled  lines,  has  slower 
phase  velocity  than  mode  2  and  it’s  effective  dielectric  constant  is  close  to  2.2. 
This  indicates  that  the  fields  are  confined  in  the  middle  layer.  On  the  other  hand, 
mode  2,  similar  to  the  even-mode  of  symmetric  coupled  lines,  has  larger  phase 
velocity.  This  implies  the  fields  around  lines  are  less  confined.  Both  modes  are 
not  strongly  dispersive.  Fig.  2.10  shows  the  current  distributions  on  each  strip 
for  both  modes. 

The  spectral-domain  analysis  provides  unique  features  in  calculating  the 
dispersion  characteristics  of  shielded  microstrip  lines.  However,  care  must  be 
taken  on  the  choice  of  the  basis  functions  for  numerical  convergence  and  efficiency. 
The  actual  unknown  current  distributions  can  be  approximated  reasonably  well 
if  a  proper  choice  of  expansion  functions  is  made.  Fig.  2.11  shows  the  comparison 
with  the  theoretical  results  published  by  Krage  and  Haddad  [13]  for  the  symme¬ 
tric  coupled  lines.  Good  agreement  is  found  if  the  same  number  of  the  expansion 
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6/Xo  5/Xo 

(a)  (b) 


Fig.  2.8  CharActeristics  of  symmetric  coupled  lines.  (Fig.  2.3(c)) 

(a)  Normalized  propagation  constant  versus  strip-ofbet. 

(b)  Characteristic  impedance  versus  strip>offiiet. 

(c)  Even-mode  longitudinal  current  density  versus  strip  position. 

(d)  Even-mode  transverse  current  density  versus  strip  position. 
€r  »  10.2,  h  *  0.02AO,  ^  *  2.0,  f  a  1.0,  f  a  20.0 
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(otuns)  Normalized  Propagation  Constant  ^ 


Fig.  2.9  Characteristics  of  aasnnmetric  coupled  lines.  (Fig.  2.3(d)) 

(a)  Normalized  propagation  constant  verstu  normalized  frequency. 

(b)  Modal  impedance  versus  normalized  frequency, 

e,  -  2.2,  ea  »  10.2,  e,  »  1.0,  -  1.0,  jjf  *  1.0, 

g-  -  1.0,  =  0.8,  jg.  =  1.0,  - 1.5,  ^  -  10.0 
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U2m/w) 


2x/w  2x/w 

(c)  (d) 

Fig.  2.10  Current  distributiona  of  uymmetric  coupled  lines.  (Fig.  2.3(d)) 

(a)  Normalized  longitudinal  current  density  versus  strip  position,  (mode  2) 

(b)  Relative  transverse  current  density  versus  strip  position,  (mode  2) 

(c)  Normalized  longitudinal  current  density  versus  strip  position,  (mode  1) 

(d)  Relative  transverse  current  density  versus  strip  position,  (mode  1) 

«i  *  2.2,  ea  *  10.2,  cj  *  1.0,  *  1.0,  Jf  =  1.0,  jg-  =  1.0,  *  0.5, 

*  1.0,  jf*  1.5,  jf  =  10.0 
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Fig.  2.11  Characteriftics  of  symmetric  coupled  lines.  (Fig.  2.3(c)) 

(a)  Effective  dielectric  constant  versus  normalized  frequency. 

(b)  Characteristic  impedance  versus  normalized  frequency. 
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Fig.  2.12  ChAracteristics  of  asymmetric  coupled  lines.  (Fig.  2.3(d)) 

(a)  Effective  dielectric  constant  versus  frequency. 

(b)  Modal  impedance  versus  frequency. 
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functions  (  first  four  Legendre  polynomials  for  expanding  the  longitudinal  cur¬ 
rent  density  and  one  sinusoidal  function  for  expanding  the  transverse  current 
density)  is  used.  By  increasing  the  number  to  8  for  the  longitudinal  current 
density  and  to  2  for  the  transverse  current  density,  the  curves  are  much  close  to 
those  using  the  expansion  functions  in  Eqs.  [2.53], [2.54]  with  only  three  expansion 
functions  for  the  longitudinal  current  density  and  two  for  the  transverse  current 
density.  The  comparison  between  this  analysis  and  the  quasi-static  results[12] 
at  low  frequencies  are  shown  in  Table  2.1  for  symmetric  coupled  lines  with  dif¬ 
ferent  strip-offsets.  The  comparison  shows  good  agreement.  Fig.  2.12  shows  the 
comparison  for  asymmetric  coupled  lines  with  the  results  published  recently  by 
Garin  and  Webb[15].  The  dyadic  Green’s  function  used  in  [15]  doesn’t  consider 
the  fact  that  the  tangential  electric  field  on  both  sidewalls  is  zero  and  is  treated 
as  a  good  approximation  only  when  the  distance  between  the  two  sidewalls  is 
large  compared  to  the  strip-width. 
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Geometrical  parameters:  €r  =  9.5« h  =  0.004A0, hi  =  gfc, w  —  h,6  =  0 

Sfh 

Vpk 

Vpk  Ref  [12] 

Zc 

ZcRef[12) 

0.05 

0.431 

0.434 

19.09 

23.50 

m 

MM 

31.37 

32.50 

0.419 

n 

41.68 

IIQI^ 

Bi 

0.412 

45.80 

Sfh 

Vph 

VpARef(12] 

Zc 

ZcRef[12] 

0.05 

0.384 

65.64 

65.60 

0.25 

0.381 

0.383 

61.89 

62.00 

1.00 

0.380 

0.382 

54.67 

55.00 

2.00 

0.384 

0.386 

50.99 

51.00 

Table.  2.1  Characteristics  of  two  layer  symmetric  coupled  lines  in 
comparison  with  the  quasi-static  results  Ref  [12]. 


Chapter  III 


Full  Wave  Analysis  of  Shielded  Microstrip 
Line-To*Line  Transitions 

3.1  Green’s  Function  Formulation  and  the  Method  of  Moments 

Several  types  of  shielded  line-to-line  transitions  as  shown  in  Figs.  3. 1*3.3 
are  investigated.  The  methodology  applied  here  is  in  analogy  to  that  reported 
in  [19], [20]  by  expanding  the  currents  in  the  coupled  line  section  with  a  com¬ 
bination  of  entire  domain  and  subdonuun  modes,  however,  the  spectral-domain 
Green’s  function  and  the  numerical  procedure  are  very  much  different.  Since 
the  strip-offset  and  strip-width  are  comparable  to  the  waveguide  dimensions,  the 
transverse  current  component  should  not  be  neglected  and  a  complete  dyadic 
Green’s  function  for  an  enclosed  multi-layer  structure  is  required  in  the  integral 
equation  formulation.  The  derivation  of  integral  equations  for  overlay-coupled 
lines  is  similar  to  that  of  Eq.  [2.49]  except  the  currents  on  strips  should  be 

+  (31) 

and 

+  +  (32) 

where 

Jii, "  (3.3) 


30 


y 


3 


incident  wave 


Fig.  3.3  Overlay  coupled-to-single  line  transition 


•'«..=  (3.4) 

nsl 

(3.3) 

j;,.  =  Ec/r.(»'),  (3.6) 

mal 

J5l.  -  (re-'^' )£'(»' +  «»/),  (3.7) 

•/».  -  i:^X(»')-  (3.8) 

n«l 

“  *^*1  (3.9) 

JSt.  =  Ew.(^').  (3.10) 

mxl 


t7(*)  = 


(3.11) 


and 

1  ;  if  X  >  0 
0  ;  otherwise. 

ovl  is  the  overlap  of  the  two  microstrip  lines.  Superscripts  e  and  s  represent 
the  entire  domain  modes  and  the  subdomain  modes  respectively.  The  transverse 
dependence  Jm  and  propagation  constants  are  parameters  with  respect  to 
the  entire  domain  modes  of  the  microstrip  t  and  can  be  obtained  through  the 
analysis  in  Chapter  II.  The  combination  of  a  piecewise*8inusoidal(P\VS)  function 
and  a  pulse  function  is  used  as  subdomain  mode  and  is  defined  from  the  end  of 
each  line(shown  in  Fig.  3.4).  These  modes  are 


0  ;  otherwise 

V  ' 


(3.12) 


and 


/S(*) 


J  ;  for  |s  -  z«|  <  } 


(3.13) 


(  0  ;  otherwise  , 

where  d  is  the  half-Ungth  of  the  PWS  function  and  also  the  width  of  pulse  func< 
tion.  Zm  is  the  center  associated  with  each  expansion  function.  The  parameter 
Jbm-  can  be  chosen  quite  flexibly  within  0  to  A  nice  feature  of  the  combination 
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microstrip  1 
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Fig.  3.4  Layout  of  expansion  modes  for 

(a)  longitudinal  currents  and 

(b)  transverse  currents 


35 


of  subdomain  and  entire  domain  modes  is  that  when  the  method  of  moments 
is  applied,  the  quantities  of  interest,  T  and  T,  are  directly  obtained  by  matrix 
inversion.  Besides,  the  dimension  of  the  impedance  matrix  is  relatively  small, 
typically  <  40.  The  integral  equation  using  the  above  currents  can  be  expressed 
as 


Eixix,y  =  hiyz) 

E2xix,y  =  hi  +  h2,z) 
E2xix,y  =  hi  +  h2,z) 


+  Jux) 
'JitxiJhx  +  ^hx) 
+  Hlx) 


0 

0 

0 

0 


(3.14) 


where 


Jitx  -  I  Jux{x)sinQnxdx\ 

Jo 

Jug  ”  j  ^cosOifX  dx ) 

4  =  °r  ju^y>-'dx\ 

00 

4  =  r 

«/-00 


and  ^  is  expressed  in  Eq.  [2.50]. 


(3.15) 

(3.16) 

(3.17) 

(3.18) 


Using  a  Galerkin’s  procedure  in  the  method  of  moments,  the  testing  func¬ 
tions  chosen  as  Jiti(x)fii{z)  are  applied  to  the  integral  equations  3.14.  The 
x-directed  currents  of  entire  domain  modes  are  neglected  in  comparison  with  the 
z-directed  currents  of  entire  domain  modes.  This  results  in  2M-('2N+2  linear 


36 


equations  with  2M+2N+2  unknowns.  These  linear  equations,  when  expressed  in 
matrix  form,  are 


[^m]  [2iu»\  l^t^] 

’[A.]  ' 
-r 

'  izir.,\ ' 

.[Z«..I  iZitL]  IZn,.]  [ZSU], 

[hx] 

T 

where  each  submatrix  [  Zjj,,  ],  due  to  the  presence  of  both  x  and  z  directed 
currents,  contains  4  submatrices  as 


(3.20) 


[Zij$$x*]  [^ijxxxx] 

and  submatrix  ( Ij,  ]  and  [  Zfj„  ]  (superscript  p  identifies  the  entire  domain  modes 
re/,  tra,  and  me),  contain  two  submatrices  as 


(3.21) 


and 


(3.22) 


Each  element  in  these  submatrices  represents  the  reaction  of  different  basis  func¬ 
tions.  For  instance,  the  elements  of  [Zi2sum]  the  reaction  between  x-directed 
currents  of  subdomain  mode  associated  with  microstrip  1  and  z-directed  currents 
of  subdomain  mode  associated  with  microstrip  2.  The  elements  of  [^um] 
the  reaction  between  x-directed  currents  of  subdomain  modes  associated  with 
microstrip  1  and  z-directed  currents  of  reflected  propagating  mode.  These  sub- 
domain  modes  are  either  PWS  functions  or  pulse  functions.  The  entire  domain 
modes  are  composed  of  the  reflected  mode,  transmitted  mode,  and  incident  mode 
which  are  distinguished  by  the  abbreviation  re /,  fra,  and  me,  respectively.  The 
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computation  of  each  element  requires  both  infinite  summation  and  integration  in 
the  spectral  domain  and  their  expressions  are  in  the  general  forms  of 

f)  (3-23) 

*“00 


4«i’=  E  r  (3.24) 

ns-oo-'"** 

where  //•  and  Jj/^  are  Fourier  transform  of  the  longitudinal  dependence  of  the  A;th 
subdomaon  mode  and  entire  domain  mode  p  respectively.  The  spectral-domain 
expressions  of  longitudinal  dependence  for  both  subdomain  and  entire  domain 
modes  can  be  found  in  Appendix  D. 


The  analysis  for  other  line-to-line  transitions  are  similar  to  that  for  overlay 
coupled  lines.  For  edge-coupled  lines  the  y  dependence  of  currents  associated  with 
the  parasitic  coupled  line  should  be  the  same  as  that  associated  with  the  feedline. 
For  coupled-to-single  lines  the  entire  domain  modes  in  the  feedlines  should  be 
discussed  separately  due  to  the  «cistence  of  two  different  modes(even-mode  and 
odd-mode). 


3.2  Numerical  Results  and  Discussions 


Numerical  results  for  line-to-line  transitions  including  reflection  and  trans¬ 
mission  coefficients  as  a  function  of  the  overlap  and  line-offset  are  shown  in 
Figs.  3.5-3.12.  It  is  found  that  the  bandwidth,  defined  as  the  frequency  range  be¬ 
tween  two  cut-off  frequencies,  is  approximately  ^  and  the  maximum  coupling  oc¬ 
curs  around  ovl  »  with  a  wide  frequency-insensitive  range.  This  implies 

that  the  transitions  are  broad-band  and  are  very  useful  in  many  MMIC  applica¬ 
tions.  These  phenomena  can  be  explained  empericaliy  by  a  quasi-static  analysis 
of  symmetric  coupled  lines.  For  two  coupled  lines,  the  impedance  matrix[21]  is 
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Reflociion  Cocflicicnt  i  F  | 


Fig.  3.5  |r|  versus  overlap  for  the  configuration  of  Fig.  3.1 
Fig.  3.6  |r(  versus  overlap  for  the  cpnfiguration  of  Fig.  3.1 
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Reflection  Coeflicient  |  F  | 


Reflection  COciTicicnt  |  T  | 


Fig.  3.9  |r|  versus  overlap  for  even-mode  excitation, 
(configuration  Fig.  3.3) 

Fig.  3.10  |T|  versus  overlap  for  even-mode  excitation, 
(configuration  Fig.  3.3) 
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Reflociion  Coeflident  |  F  | 
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- 1 - 1 
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Fig.  3.11  |r|  versus  overlap  for  odd-mode  excitation, 
(configuration  Fig.  3.3) 

Fig.  3.12  |T|  versus  overlap  for  odd-mode  excitation, 
(configuration  Fig.  3.3) 
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(3.25) 


and 


where 


and 


Z\i  —  Z22  —  “jf  ZcoCOtOo) 

Z\2  ~  Z21  —  ~3 ^.ZcgCscOf  ~  ZggCsc 


(3.26) 

(3.27) 


0^  =s  ^ovi 


(3.28) 


Zee  s^nd  Zco  represent  respectively  even-mode  and  odd-mode  characteristic  imped¬ 
ance  associated  with  the  coupled  line  section.  Tlie  scattering  matrix  can  be 
found  from  the  impedance  matri.x.  As  a  result  the  reflection  and  transmission 


coefficients  are  obtained: 


z?.-z?.-z; 

(Z„  +  Z.)»-Z?a’ 


(3.29) 

(3.30) 


'  ■  (Z,.  +  Z.)«-Z?,’ 

where  Ze  is  the  characteristic  impedance  of  feedine  and  the  parasitic  coupled  line. 
From  Eq.  [3.30],  the  cut-off  condition  occurs  at  Zu  s  0.  That  is 

(3.31) 

Zeo  Sin  00 

For  the  general  case  Zee  »  Zeo,  the  solution  is  approximated  to 


nsO,  1,2,3,  < 

m 


(3.32) 


Since  the  fields  of  the  odd-mode  associated  with  the  coupled  line  section  in  an 
overlay  configuration  are  nM)stly  confined  in  the  middle  layer,  the  effective  dielec¬ 
tric  constant  of  the  odd-mode  is  approximately  equal  to  the  dielectric  constant 
of  the  middle  layer.  Therefore,  the  overlap  for  cut-off  is  about: 


For  edge  coupled  lines,  the  fields  associated  with  odd-mode  are  spread  in  the  two 
dielectric  layers.  Therefore,  the  relative  guided  wavelength  is  approximately  equal 
to  the  mean  value  of  the  guided  wavelength  with  respect  to  each  dielectric  layer. 
Assuming  these  two  dielectric  layers  with  permittivity  e^,  and  ej  respectively,  the 
overlap  for  transmission  cut-off  is  about: 


Figs.  3.5-3.8  show  that  the  coupling  efficiency  is  better  in  overlay  line- 
to-line  transition  than  in  edge  coupled  line-to-line  transition.  This  indicates  the 
former  constitutes  a  promising  element  in  realization  of  millimeter  wave  high-pass 
filter.  Figs.  3.9-3.12  show  the  results  for  the  case  of  an  overlay  coupled-to-single 
microstrip  transition.  It  is  seen  that  the  even-mode  coupling  depends  less  on  the 
strip-offset  of  parasitic  coupled  line,  as  compared  with  single-line  coupling.  It  is 
also  noted  that  the  even-mode  of  coupled  lines  can  couple  energy  to  a  centered 
parasitic  microstrip  line  while  the  odd-mode  can  not.  This  may  find  applications 
in  a  phase  detector. 

The  frequency  response  of  an  overlay  coupled  microstrip  transition  is 
shown  in  Figs.  3.13-3.14.  The  geometrical  parameters  are  specially  designed 
in  the  coupled  line  section  where  the  line-width  is  much  larger  than  the  spac¬ 
ing  between  two  lines.  It  is  seen  that,  in  a  wide  frequency  range,  the  coupling 
coefficients  are  almost  independent  of  frequency.  In  addition,  it  is  possible  to 
couple  more  than  95%  of  the  total  power  through  the  discontinuity.  This  geom- 
etry(so  called  suspended  stripline)  may  be  very  useful  due  to  these  two  excellent 
characteristics. 
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Phase  (degree)  Magnitude 


Fig.3J4 

Fig.  3.13  Magnitude  of  F  and  T  versus  frequency. 

Fig.  3.14  Phase  of  F  and  T  versus  frequency. 

Cl  s  <3  s  1,  €2  SB  10.2,  h|  s  hs  s  70mil,  ha  s  lOmil,  wi  s  wa  s  SOmil, 
as  ^3  s:  0,ov/  as  185mt7,a  sa  500mil,  (configuration  Fig.  3.2) 
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In  the  present  computations  of  the  above  transition  problems,  entire  do¬ 
main  mode  s  of  3  guided  wavelengths  long  and  9-18  subdomain  modes  for  both 
pulse  and  piecewise  sinusodial  fimctions  are  used  in  each  microstrip  line.  The 
numerical  results  obtained  are  within  1%  convergence  accuracy.  The  power  con¬ 
servation  is  also  checked  in  each  case  and  will  be  discussed  in  next  section. 

3.3  Power  Conservation  Check 

For  shielded  microstrip  transitions,  the  convergence  of  the  moment  method 
solutions  is  very  sensitive  to  the  type  and  number  of  expansion  function  chosen. 
Power  conservation  provides  a  nice  way  of  checking  the  accuracy  of  the  solutions. 
According  to  the  configurations  in  Figs.  3.1-3.3,  the  incident  power  should  be 
equal  to  the  summation  of  reflected  power,  transmitted  power  and  some  losses 
coupled  to  the  higher  order  modes  and  the  multi-layered  waveguide  modes.  These 
higher  order  modes  and  multi-layered  waveguide  modes  correspond  respectively 
to  the  zeros  and  poles  of  the  characteristic  equation  2.64[22].  Moreover,  all  of  the 
higher  order  modes  and  multi-layer  waveguide  modes  have  cut-off  frequency  while 
the  dominant  mode  is  zero  cut-off[23].  That  is  to  say,  by  choosing  the  relatively 
smaller  waveguide  dimensions  or  lower  frequencies,  one  can  allow  only  dominant 
mode  propagating  in  the  shielded  structure.  In  this  analysis  to  characterize  the 
shielded  line-to-line  transitions,  the  waveguide  dimensions  are  restricted  to  allow 
only  dominant  mode  for  each  line.  Therefore,  the  relation  of  power  conservation 
for  line-to-line  transitions  can  be  simply  written  as 

Pine  ~  Pref  "I"  Ptrai  (3.35) 

where  ■^'nei  P ref  and  Ptra  represent  incident,  reflected  and  transmitted  power, 
respectively.  With  the  reflection  and  transmission  coefficients  obtained  by  the 
method  of  moments  in  Sec.  3.1  and  the  characteristic  impedance  obtained  by  the 
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Poynting  vector  analysis  in  Sec.  2.3.  The  reflected  power  associated  with  the 
feedline  and  transmitted  power  associated  with  the  parasitic  coupled  line  can  be 
normalize  respectively  by  the  product  of  the  characteristic  impedance  and  the 
square  of  the  absolute  value  of  the  scattering  coefficient  for  each  line.  Thus,  the 
relation  of  power  conservation  can  be  expressed  in  terms  of  scattering  coefficients 
and  characteristic  impedance  as  follows: 

|r|»  +  =  1,  (3.36) 

"el 

where  Zd  and  Zc2  are  the  characteristic  impedance  of  the  feedline  and  the  para¬ 
sitic  coupled  line,  respectiveh*.  For  coupled-to-single  configuration,  the  relation  is 
slightly  different  due  to  different  definition  of  characteristic  impedance  for  sym¬ 
metric  coupled  lines.  The  relation  is 

iri’  +  =  1,  (3.37) 

where  Zee  wd  Zeo  are  characteristic  impedance  of  even-mode  and  odd-mode, 
respectively. 

The  conservation  check  for  Figs.  3.5-3.12  are  shown  in  Table  3.1‘^3.4  re¬ 
spectively.  It  is  seen  that  the  error  in  each  case  is  within  1%.  This  confirms 
the  excellent  accuracy  of  the  numerical  results  of  both  the  frequency-dependent 
characteristic  impedance  and  the  moment  method  solutions. 
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Microstrip  1 

=  0.075a 
/3pi  =  2.634 
Zei  =  45.90 

Microstrip  2 

62  =  0.075a 
Pp2  =  2.634 
Zc2  ”  45.90 


ovl 

iri’  +  fem’ 

ovl 

|r|*  +  |«|Tp 

•0.05 

1.000 

0.13 

1.000 

-0.03 

1.000 

0.15 

1.000 

-0.01 

1.000 

0.17 

1.000 

0.01 

i.000 

0.19 

1.000 

0.03 

1.000 

0.21 

1.000 

0.05 

1.000 

0.23 

1.000 

0.07 

1.000 

0.25 

1.000 

0.09 

1.000 

0.27 

1.000 

1.000 

0.29 

1.000 

Table.  3.1  Power  conservation  check  for  the  configuration  of 
Fig.  3.1.  Both  |r|  and  |T|  are  the  same  as  those  in  Fig.  3.5 
and  Fig.  3.6. 
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Microstrip  1 

6i  =0.000 

/3pi  =  1.946 
Zci  =  55.90 

Microstrip  2 

^2  =0.000 
/3p2  =  1.907 
Zc2  =  59.00 


ovl 

in^+l^iTp 

ovl 

.0.05 

1.000 

0.13 

0.997 

.0.03 

1.000 

0.999 

.0.01 

1.001 

0.17 

1.004 

0.01 

1.003 

wm 

1.003 

0.03 

1.004 

0.21 

0.05 

1.003 

0.23 

1.001 

0.07 

1.001 

0.25 

1.000 

0.09 

1.001 

0.27 

1.000 

0.11 

0.999 

0.29 

0.999 

Table.  3.2  Power  conservation  check  for  the  confignrati  on  of 
Fig.  3.2.  Both  |r|  and  |T|  are  the  same  as  thos  e  in  Fig.  3.7 
and  Fig.  3.8. 
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Microstrip  la  and  lb 

=  0.100a 
=  1.950 
Zci  =  63.80 

Microstrip  2 

^2  =  0.000 
/3p2  =  2.072 
Zc2  “  53.60 


ovl 

|r|2  +  ^\T\^ 

ovl 

iriHfem* 

-0.05 

1.000 

0.13 

1.001 

-0.03 

1.000 

0.15 

1.002 

•0.01 

1.000 

0.17 

1.002 

0.01 

1.000 

0.19 

1.001 

0.03 

0.21 

1.002 

0.05 

1.002 

0.23 

1.005 

0.07 

1.003 

0.25 

0.09 

1.002 

0.27 

0.11 

1.001 

0.29 

1.003 

Table.  3.3  Power  conservation  check  for  the  configuration  of 
Fig.  3.3.  Both  |r|  and  |T|  are  the  same  as  those  in  Fig.  3.9 
and  Fig.  3.10. 
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Microstrip  la  and  lb 

51  =  O.lOOo 
jdpi  =  2.291 
Zei  =  45.30 

Microstrip  2 

52  =  0.100a 
I3p2  =  2.075 
Zc2  =  53.60 


ovl 

|r|=  +  ^|Tp 

ovl 

-0.05 

mm 

I.0QS 

-0.03 

1.000 

mm 

1.0QS 

• 

P 

0 

1.000 

0.17 

1.003 

0.01 

1.000 

mm 

1.001 

0.03 

1.000 

0.21 

0.05 

1.000 

0.23 

1.000 

0.07 

1.000 

0.25 

1.002 

0.09 

1.001 

0.27 

1.0QS 

0.11 

1.003 

0.29 

1.007 

Table.  3.4  Power  conservation  check  for  the  configuration  of 
Fig.  3.3.  Both  |r|  and  |r|  are  the  same  as  those  in  Fig.  3.11 
and  Fig.  3.12. 


Chapter  IV 


Conclusions 

An  efficient  method  has  been  presented  in  Chapter  II  for  obtaining  the 
dispersion  properties  of  three  layer  shielded  microstrip  lines.  This  method,  which 
is  based  on  a  hybrid  mode  analysis  followed  by  the  method  of  moments  has 
a  number  of  unique  features.  Numerical  results  given  include  the  propagation 
constant,  characteristic  impedance,  transverse  dependence  of  longitudinal  and 
transverse  current  distributions  for  both  single  and  coupled  microstrip  lines  with 
arbitrary  strip-width  and  strip-offset.  The  dispersion  characteristics  and  their 
dependence  on  various  geometrical  parameters  have  also  been  presented.  This 
method  is  general  and  can  be  easily  extended  to  the  multiconductor  case.  The 
choice  of  expansion  functions  is  carefully  made  for  improving  numerical  efficiency. 
Good  agreement  is  shown  in  comparison  with  other  available  data. 

In  Chapter  III,  a  full-wave  analysis  is  proposed  to  develop  a  generalized 
dynamic  model  for  several  types  of  shielded  microstrip  line-to-line  transitions. 
The  results  obtained  from  the  method  of  moments  are  checked  by  power  conser¬ 
vation  with  error  within  1%.  The  results  presented  also  show  excellent  properties 
in  some  transitions  and  may  find  promising  applications  in  MMIC  coupler  and 
filter  designs. 
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Appendix  A 


The  pertinent  coefficients  in  Eqs.  [2.18]~[2.29]  are  expressed  as 


^  _  -Jw/iosina„a:;  /e,  .  , 


(A.1) 


,  -  I)sinhqihismhq2h2  fil 

-jufio  sincK„Xi 

jfc3  Den2  ~  ^2^2  /21 

1)  cosh  91^1  /n 


(A.2) 


-jufip sinofnJ^l  [  cosh 92(^1  + /i2)cosh93A3 

kl  Deni  y_  sinh  93(^1  +  ^2)  sinh93h3 


(A.3) 


jufip  sin  Ofn^i  [  ^sinhqjhisinh  92(^1  +  ^2)  cosh  93^3 
k^  Deni  1  ^  cosh  92(^1  +  ^2)  8>nh  93^3 


(A.4) 


■"  ^3)  /31 

•  sinh  9iAt  cosh  92^1  cosh  93^3 
-  «i)  /ji 

ja>/io  sina^Xi  •cosh9iAi  sinh  93^1  cosh  93^3 

>cosh9i  sinh  93(^1  +  ^3)  cosh  93^3 
-i^S(c3  -  ei)  /n 

V 'C08h9ihi  cosh  93(^1  +  ^3)  sinh  93^3 


(A.5) 


5S 


c;= 


-lun,  3ina,i, 

kl  Dtnl  Dtn2 


-  ca)  /21 

•  sinhqihi  sinh92Ai  cosh^s^s 

~  ^0  /21 

•  cosh  qihi  cosh  gjhi  cosh  93^3 
+i^(€3  -  fn 

•  cosh^i  cosh  93(^1  +  A3)  cosh  93/13 

+i^^(c2  -  Cl)  /ll 

cosh^iAi  sinh  93(^1  +  h2)s\nhq^h3 


(A.6) 


D  =  sinongi  I  if  JJco8h93A3sinh93A3 
kl  Dtnl  y  ^.asinh93A3Cosh93A3 


(A.7) 


,  (i)®5Jc3(|f  “  l)«nh93A38inh93A3 /u 

•jujfio  sin  ori|X| 

kj  DtnlDen2  -  1)  c<»h93A3C08h93A3  fn 

~  1)  cosh 93^3  /31 


(A.8) 


-JUHoCOtQni 

kl  Dtnl 


(A.9) 


»  ( -OnJcaCJf- l)»inh9iAi8mh93A3 /31 

-jufio  cos  anX| 

DtnlDtn2  -“»*«*(«"  1) <=<*^91^1  co8h93A3  /21 
^+«n«i(if-l)co8h9iAi  fn 


(A.IO) 


B  =s  cosonxl  I -J-8inh9iAico8h93(Ai  +  hfjcoaliqzhs  \ 
kl  Dtnl  I  _3inh9iAi8inh93(Ai  +  A3)8inh93A3  ) 


C  =s  (  J«ot9iAi8inh93(Ai  +  h2)cothq3h3 

kl  Dtnl  I 


(A.12) 


5« 


“OtnfJ‘(C2  —  C3)  /2I 
•  sinh qihi  cosh  92^1  cosh 


cosan^l 

*  kl  Den\  Den2 


•coshgiAi  sinh  92^1  cosh  93^3 
-an(c3  -  fn 


•  coshfi  sinh  92(^1  +  ^2)  cosh  93^3 


\ 


-«ng(C2  -  «l)  fll 

cosh hi  cosh 92(^1  +  ^2)  sinh 93^3  j 


(A.13) 


c'  ss  cdaofngi 

*  h|  Deni  Den2 


^®a*|-(c2  -  C3)  /21 
*8inh9ihi  sinh  93^1  cosh  93^3 
+an(Sif^-ei)/2i 
•  cosh  9ihi  cosh  92^1  cosh  93^3 
+o«(<s  -  fu 


•  cosh  9i  cosh  92(^1  +  ^2)  cosh  93^3 


\ 


+tt»g(<3  ~  «i)  fn 

y •  cosh 9i hi  sinh 92(^1  +  h2)8inh93h3  ^ 


(A.14) 


-jufip  cos  a^x[  I  Jf"- cosh  93h3  sinh  92h2 
h|  Deni  ^  +  si  sinh  93h3  cosh  92h2 


(A.15) 


1'  3s  COSOnXi 

*  kl  Deni  Den2 


®nj«3(-f  “  I)sinh93h38inh92h2  /u 

+Oa«3(if  -  1)  cosh 93h3 cosh 92h2  fn 


\ 


\-On«3(U  -  1)  coshgshs  /21 


(A.16) 


where 
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Den\  = 


i)en2  = 


9i  cosh qihx  cosh  92^2  sinh  q^hz 
cosh  qx  hx  sinh  qzh^  cosh  93^3 

-2 

+93  sinh  hi  cosh  92^2  cosh  93^3 
^+92  sinh 9i hi  sinh  92^2  +  sinh93h3^ 

£391  sinh9ihi  C08h92h2  cosh93h3 
+C2®^  sinh  9ihi  sinh  92h2  sinh  93h3 
+^92  cosh  9i  hi  sinh  q2h2  cosh  qzhz 
^  +C193  cosh  9i hi  cosh  92h2  sinh 93h3  ^ 


(A.17) 


(A.18) 


fix  =  sinh9ihi  cosh  92  h3  sinh  93h3  +  —  sinh9ihi  sinh  92/(2  cosh  93h3  (A.19) 

92 

fn  —  ^  cosh  9ihi  sinh  92h3  sinh  q^hz  +  sinh  9ihi  cosh  93h2  sinh  q^hz  ( A.20) 

92 

/i3  =  sinh  9i  hi  sinh93h3  (A.21) 

/ai  *  sinh  9i hi  sinh93h3  (A.22) 
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Appendix  B 


The  spectral-domain- dyadic  Green’s  function  in  Eqs.  [2.36]  [2.43]  is 

,  _  ^  ^  rr  .  .  r  rs] 


kj  [  Deni 


G.i„ - ^ 

Gu«  “  ■*■  Deni  "Den2^^“  ' 


G..,.  -  ;ct/u  +  DmlDen2^^"  '  (®'^) 

111*  =  — G'xiir  (B.4) 

(®-«) 

_  W/iof-On^^  ,  an/5  ^  ^  ,xl  „x 

+  Dtnl  Dei.2^-^*' '  ’  •'»>] 


/7enl  Z?en2 


(/ii  •  /s  +  /21  •  /s) 


=  -G, 


Gi2,,  = 


fl2xr  = 


-jufio  \kl-^^ 


"'2  M  f  j _ H 

kl  Deni  Deni  Den2 


-ji^fip  [^2 -Qn 


kl  [  Deni 


—  “G^12j» 


G„„  - 

G22XX  =  •~^?22» 

^21,*  =  C7i2*» 

G’21*x  =  Gi2ax 
G2IX*  —  Gi2xx 
G2IXM  —  Gi2xt 


(B.8) 

Ga..  -  ^^J-f„  +  —J—(f„.f,  +  f^,.f^)  (B,9) 

G,a.  =  +  DcnlDcn2^^^ '  A  +  /ir  A)  (B.10) 

^  +  DtnlDm2^^” '  <®'“' 


(B.11) 

(B.12) 

(B.13) 

(B.14) 

(B.15) 

(B.16) 


where  Denl,Den2t  /n,  /12,  /21,  /22  are  listed  in  Appendix  A,  the  functions  /a,  /«,  /s, 
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and  /e  are  expressed  as  follows: 


/s  =  — 


^2 

•C39i( - 1)  sinh^i^i  cosh  92^2  cosh  93^3 

C3 

— ^2(- —  1)  cosh  qihi  sinh  92^2  cosh  03^3 

^2  C3 

=  ~«i93(—  —  1)  cosh^ihi  cosh92^2  sinh  93^3 
~  ~^92(—  —  1)  cosh  hi  sinh  92^2  cosh  03^3 

t2  61 


h  —  ~^i93( - 1)  cosh  hi  sinh  ^3^3 


/s  —  — C39i(—  —  l)sinh<jrihi  cosh93h3 
C3 


(B.17) 

(B.18) 

(B.19) 

(B.20) 
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Appendix  C 

The  average  power  in  each  region  is  expressed  as 

(C.1) 

nsB-oo 

~  ^2v*^2xe)f^*e 
+(-^2xc-fl2y,  ~  ^2yeH2gg)f2e» 

+(-^2xc^2yc  ~  ^2yc^2xe)f^ee 

^3  “  ^  ^  [-^ax^sy/a**  -  ^ay^ax/scc]  (C.3) 

™  ns—oo 

where 


in,  =  (t;  -  al)D,  +  +  d'.)  -  jfia,D, 

(C.4) 

ii,  =  (kl  +  ql){D',  +  D',)  -  a,q,D,  -  jqiffD, 

(C.5) 

-  uei0{D',  +  D’^) 

(C.6) 

w 

/fly  =  Uti^Dx  -  jwCiOn-D* 

(C.7) 

^2x*  =  (^2  ~  ®n)^»  +  «n92(C'*  +  C^x)  “  jl^<Xn^x 

(C.8) 

^2xc  =  +  o„g2(5l  +  Bl)  -  ;/9anC, 

(C.9) 

&2y,  =  (fcj  +  92)(-®*  +  B^)  —  OCn<l2Cx  ~  j^xPCx 

(C.IO) 

hyc  =  (A:|  +  9|)(c;  +  c;)  -  On925x  -  kxPBx 

(C.ll) 

^2x»  =  -  u;e2/3(B'  +  5^) 

(C.12) 

^2xe  —  J^f’iq^Bx  —  U}€2fi{Cx  +  C'x) 

(C.13) 

^2y*  =  WC2/S5*  -  jlj}t20lnBx 

(C.14) 
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k^yc  =  U}€i/3Cr  -  iwejOnC* 

(C.15) 

■^3x  =  (*3  -  «n)^*  -  OlnqaiK  +  K)  “ 

(C.16) 

£^3v  =  (^3  +  +  -^r)  +  +  jqz0At 

(C.17) 

^3*  =  -J^^zq^A,  -  ijJ€zP{X^  +  i4^) 

(C.18) 

w 

/fsy  =  ^tzpAx  —  juezOnAx 

(C.19) 

and 

/iM  =  -^^n;  +  “-sinh2?iAi 

1  2  Aqi 

(C.20) 

flee  -  -  +  7— smh2?iAi 

2  4gi 

(C.21) 

f2$s  -  Sgrii  -^^2  +  ~(sinh2?2(Ai  +  A2) -"Smh2g2/ii) 

2  4^2 

(C.22) 

/sm  =  —(cosh  292(^1  +  Aj)  “  cosh  292^1) 

(C.23) 

fzet  =  ■S'^n,-  j^(co8h292(hi  +  ^2)  ”CO»h292hi) 

(C.24) 

/jec  *  ih3  +  ~(8inh  292(^1 +  h2)-«nh  292^1) 

2  492 

(C.25) 

fzx,  «  S^n,-  -ihs  +  ri-sinh  293^3 

L  2  493 

(C.26) 

fzee  *  ;rh3  +  7~sinh293h3 

2  493 

(C.27) 

Sgrti  is  equal  to  1  as  9,-  is  real  and  equal  to  -1  as  9,-  is  imaginary. 
Ail  ~  and  D'^  can  be  found  in  Appendix  A. 

The  functions 
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Appendix  D 


The  Fourier  transform  of  longitudinal  dependence  of  both  subdomain 

modes  and  entire  domain  modes  are  as  follows: 

7m  2fce,(cOS  Kid  -  COS  ^d) 

zmk^d(lP-kl) 

(D.l) 

2  sin  ^  zg 

2 

(D.2) 

+  i)  /  sin ^pxz^'^*dz 
«/-00 

(D.3) 

+  j)  1  sinfip2Z6'’^^*dz 
^-00 

(D.4) 

^i»~W  =  -  j)  sin 

(D.5) 

where 

sin  0^zc^‘>-dz  =  +  ;>  (f(/3  ±  0^)  -S(0:f0^)]  (D.6) 


